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1.1 Introduction 

Axions are well motivated particles proposed in an extension of the Standard Model 
(SM) as a solution to the CP problem [1] in strong interactions. On the other 
hand, there is the category of axion-like particles (ALPs) which appear in diverse 
extensions of the SM and share the same phenomenology of the axion. Axions and 
ALPs are hypothetical neutral particles which interact weakly with matter, being 
candidates to solve the Dark Matter (DM) problem. 
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Several experiments have been developed for axion and ALPs detection, based 
on laboratory regeneration searches, galactic halo searches (haloscopes) and solar 
searches (helioscopes), which will be briefly discussed in this chapter. The helioscope 
technique, which is the main topic of this work, will be explained in detail. 

1.2 The strong CP problem in QCD 

The Quantum Chromodynamics (QCD) is a gauge field theory that describes the 
properties of the strong interactions between quarks and gluons, being the SU(3) 
component of the SM. A problem appears in QCD, the so-called strong CP problem. 
The theory predicts a CP-symmetry violation in strong interactions which have not 
been observed experimentally. The strong CP problem is closely related to the U(l)u 
problem of QCD and to its solution proposed by’t Hooft in 1976, some features of 
this problem will be briefly described is this section. 

1.2.1 The U(l)u problem and its solution 

The Lagrangian of QCD can be written as: 

Cqcd = - ran)Vn - \g%G^ ( 1 . 1 ) 

n 

where G a ^ u are the field strength tensors of the gluon fields with a color index 
a=l,...,8; are the quark fields with the quark flavor n and m n the corresponding 
quark masses. D^ is the covariant derivative defined as D^ — igT a G^ L where g 

is the coupling constant and T a are the generators of the group of rotations of the 
quark fields in the color space. 

The QCD Lagrangian has a large global symmetry in the limit of vanishing quark 
masses m n — > 0; U(N)y x U(N)u, thus it is invariant under global axial and vector 
transformations. In reality u and d quarks masses could be considered approximately 
zero because m u ,md <C A qcd, being A qcd the dynamical scale of the theory. So 
strong interactions are expected to be approximately U(2)y X U(2)u invariant. 
Indeed, U(2)y = SU(2)y x U(l)y is the vector symmetry corresponding to isospin 
times baryon number, being a good approximate symmetry in nature. 

However, the corresponding axial symmetries have not been observed. This phe¬ 
nomenon is understood because the SU(2)u symmetry is not preserved by the QCD 
vacuum and by spontaneous breaking we get the Nambu-Goldstone bosons (the pi- 
ons), whose mass vanishes in the limit m Ut d —> 0. But there is no pseudoscalar 
state with vanishing mass corresponding to the Nambu-Goldstone boson of a U(l)u 
symmetry, this is known as the U(l)u problem [2], 

The resolution of the U(l)u problem came through the realization by t’Hooft [3, 4] 
that the QCD vacuum has a more complicated structure. The problem was bypassed 
by introducing an anomalous breaking of the U(1 )a symmetry, adding an extra term 
Cg to the Lagrangian: 
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£o = o£r 2 G%G^ ( 1 . 2 ) 

here G° IV is the gluon field strength tensor and Ga v its dual, g the coupling constant 
and 9 is an arbitrary angle between 0 and 2ir, which arises from the complicated 
structure of the QCD vacuum. More precisely, this means that the ground state is 
a superposition of an infinite number of degenerate states | n). The ground state, 
also called 0-vacuum can be expressed as: 

OO 

m= £ e-"l»> (1.3) 

n=— oo 

Taking into account the electroweak interaction, the angle 6 in equation 1.2 has to 
be substituted by 9: 


9 = 9 + 9 weak = 9 + arg(detM) (1.4) 

where M is the quark mass matrix, thus the Lagrangian becomes: 

C ~e = ^£ 2 °%^ (1-5) 

and the QCD Lagrangian can be written as: 

GqCD = Gpert + Gg ( 1 . 6 ) 

This additional term in the Lagrangian solves the 11(1)^ problem but is not invariant 
under CP, which leads to the strong CP-problem. 


1.2.2 The strong CP problem 

For a non zero choice of 9 in equation 1.5 the CP symmetry is expected to be 
violated as it actually occurs in electro-weak interactions, but experimentally the 
CP violation is not observed in strong interactions. This is the case of the electric 
dipole moment of the neutron d n (nEDM) which has been estimated by the MIT 
(Massachusetts Institute of Technology) bag model [5]: 


d n = 32.7 x 1(T 


" m u m d + m u m s + m d m. 


(1.7) 


here R is the bag radius R ~ (140 MeV) -1 , the quark mass ratios are estimated 
to be md/ni u = 1.8, m s /m.d = 20 and the bag value m s ~ 300 MeV, obtaining 
d n = 8.2 x 1O~ 16 0 e cm. 

Experimentally the nEDM is bounded to \d n \ < 2.9 X 10~ 26 e cm (90 % C.L.) [6], 
it constraints 9 < 10 -10 . Such small value of 9 is allowed in the theory, however 
the problem is beyond the unexplained smallness of 9. It implies that the two 
contributions of equation 1.4, which are in principle unrelated, have to cancel each 
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other with a precision of 10 10 . The strong CP problem is why this 9 angle coming 
from the strong and weak interactions, is so small. 

There are mainly three solutions proposed in order to solve the CP-problem: the 
first and least likely implies zero quark masses, the second one is to set 9 = 0, 
imposing CP symmetry in the QCD Lagrangian and the third and more compelling 
solution is the Peccei-Quinn mechanism that will be described below. 

1.2.3 The Peccei-Quinn solution 

The Peccei-Quinn (PQ) solution was proposed in 1977 [7, 8]. In order to equal 9 
to zero it postulates a new global and chiral symmetry U(1)pq (the PQ symmetry) 
that is spontaneously broken at the energy scale of the symmetry f a . It implies the 
existence of a new held a which appears as the pseudo Nambu-Goldstone boson of 
the new symmetry, the axion. This new held yields an additional term C a in the 
QCD Lagrangian Cqcd 


f-'QGD — Cpert T Cg T C a 


( 1 . 8 ) 


where C a is given by 


C a = -^d^ad^a + C int + C a ° 


9 


s~ifJLV 

fa 32M^ a 


(1.9) 


The hrst term take account of the kinetic energy, the second term represent further 
interactions of the axions and the third term introduces the axion held a and its 
coupling to gluons. C a is a parameter dependent of the model, g is the strong 
coupling constant and f„ is the energy scale of the spontaneous breaking of the 
PQ symmetry. The third term in equation 1.9 provides an effective potential for 
the axion held 14 //, so the vacuum expectation can be obtained by calculating its 
minimum. 


d!4// \ _ _ r< g" /^a 
da / Ca 32vr 2 /a ' 


= 0 


(a) 


( 1 . 10 ) 


Thus, the value for the vacuum expectation of the axion held (a), can be written as: 

(a) = -^9 (1.11) 

^a 


It solves the strong CP problem by canceling the 9 term for any value of f a , providing 
a dynamical solution to the problem. Expanding 14// around its minimum the axion 
get mass: 


m 


2 

a 


/ 9 2 Veff \ 

\ d 2 a / 


C a 


9 2 d 

32ir 2 f a da 


KM") 


(a) 


( 1 . 12 ) 
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So far the PQ mechanism is the more elegant solution to the strong CP problem 
and could be experimentally proven by the discovery of a new particle: the axion. 

1.3 Axion properties. 

The PQ solution fixes some properties of the axion that only depend on the energy 
scale of the spontaneous breaking of the symmetry f a . These properties are the 
axion mass m a and the coupling constant g a i of the axion to other particles denoted 
by i, both of them are inversely proportional to f a . 

1 1 . 

m a oc , g a i oc — (1.13) 

Ja Ja 


1.3.1 Axion coupling to matter 

Although some of the other properties of the axions are model-dependent, generically 
the axion couples with gluons and photons, which is a consequence of the former. 
Axions could also interact with fermions. Some features of these interactions will 
be described below. 

1.3.1.1 Coupling to gluons 

The coupling of axion with gluons is given by the third term in equation 1.9, by the 
expression: 


£aG 


CX-S 

8-irfa 




(1.14) 


with the strong fine-structure constant a s . Due to its coupling with gluons, axions 
can acquire mass by mixing with pions 


m a = 


tUtt/tt I 5 10 9 GeV 

fa \J (1 +Z + w)(l + z) ~ me fa 


(1.15) 


here m n = 135 MeV is the pion mass and f n ~ 92 MeV the pion decay constant, 2 
and w are the quark masses ratios m u /md and m u /m s respectively. The coupling of 
gluons to axion is generic to all the models, thus the axion mass mass, given by 1.15, 
is inversely proportional to f a . 


1.3.1.2 Coupling to photons 

Axion couples to photons through the mixing of pions with axions. The contribution 
of the axion-photon interaction to the Lagrangian is given by the expression: 


£ 


aj — 


/ ^9a r y FiCL — g a ^EBa 


(1.16) 
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here g a7 is the coupling constant from axions to photons, F /IU is the electromagnetic 
field-strength tensor and its dual, E and B are the electric and magnetic fields 
and a is the axion field. A further contribution in the coupling can appear in models 
in which standard fermions carry PQ-charges in addition to the electric charges. The 
axion-photon coupling constant can be written as: 


a ( E 2(4 + z + w)\ 

gai 2irf a \N 3(1 + z + w)J 


(1.17) 


here a is the fine structure constant, E/N is the ratio between the electromagnetic 
E and color N anomalies, z and w are the quark mass ratios introduced in equa¬ 
tion 1.15. Since E/N is a model dependent factor, the terms in the parenthesis of 
equation 1.17 are usually presented in literature as the coefficient C 7 . 


9ay 


a 

Infa 


C y , 


, _ f E 2(4 + z + w)\ 
7 “ \N “ 3(1 + z + w)J 


^-1.92 ±0.08^ (1.18) 


The axion-photon coupling is generic to all the models and most of axion searches 
strategies are based on this interaction. 


1.3.1.3 Coupling to fermions 

Depending on the model, axions can also couple with fermions, this interaction is 
postulated in the second term of the Lagrangian from equation 1.9 1 , in which the 
interaction of an axion to a fermion / is given by: 


(vnVf) 


(1.19) 


where Tj is the fermion field, nif the fermion mass and g a f the axion-fermion 
coupling constant, that can be written as [9]: 


9af 


C f m f 

fa 


( 1 . 20 ) 


The dimensionless combination of g a f plays the role of a Yukawa coupling with an 
effective PQ-charge Cf and a fine-structure constant of the axion a a f = g'/j/A t: can 
be defined. 

The axion could also couple to electrons at tree level, this process is possible if 
electrons carry PQ-charge, thus the coefficient C e 0. This is the case of the DSFZ 
model (see section 1.3.2.2) in which the axion-electron coupling constant g/ffi 6 is 
given by 


9 


tree 

ae 


^' - = 0.85 x 10 w m a C e eV 1 


( 1 . 21 ) 


1 Note that the second term Ci„t in the Lagrangian of equation 1.9 is defined as C int = C al + C a f 
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Even if the effective PQ-charge C e = 0, axions could couple to electrons in a higher 
order. This is the so-called radiatively induced coupling at the one loop level, 
however the <f 0 f d coupling is weaker than g^f e . 

Finally, an effective axion-nucleon coupling could be considered since there are no 
free quarks below the QCD scale A qcd ~ 200 MeV. There are two contributions to 
the axion-nucleon coupling: the coupling to light quarks at tree level and the mixing 
with pions. The coupling constant is given by 

9aN = = 1.57 x 10 ~ 7 m a C N eV 1 (1.22) 

Ja 

here Cat is a model-dependent coefficient of the interaction in which the coupling to 
protons and neutrons have different contributions. 

1.3.1.4 Further processes 

There are processes involving axions quite relevant in the frame of astrophysics. 
The predominant emission process of axion in stars is the Primakoff effect 7 + Ze —> 
Ze + a, in which a photon can be converted into an axion in the presence of strong 
electromagnetic fields (e.g. the electric held of the charged particles inside the 
plasma). This process is quite relevant in axion searches because axions could be 
reconverted into photons (and detected) inside a strong magnetic held via inverse 
Primakoff effect. 

Further processes in axion models with tree-level coupling to electrons might dom¬ 
inate the axion emission in white dwarfs and red giants. These processes are the 
Compton 7 + e _ — > e~ +a and the Bremsstrahlung emission + Ze —> Ze + e~ + a. 
Also, the axion-nucleon Bremsstrahlung N + N —> N + N + a could be relevant in 
supernova explosions. 

These processes and its implications in astrophysics will be revisited in section 1.3.3 

1.3.2 Axion models 

Axion could couple to different particles as it was presented in the previous section. 
These interactions are strongly model dependent and can be divided in two different 
trends: the hrst postulated axion also referred to visible or PQWW (Peccei-Quinn- 
Weinberg-Wilczek) axion and the invisible axion model which assumes a lighter 
axion. 

1.3.2.1 Visible axion model 

In the original PQ model [10], the U(1)pq symmetry was broken at the electroweak 
scale f a = A wea k with A wea k — 250 GeV. It implies axion masses larger than 150 keV 
and relative strong coupled to baryonic matter and radiation. Such axions would 
have been produced and detected in reactor and accelerator experiments. The non 




8 


Chapter 1. Axions and Axion Like Particles 


experimental observation of the PQWW axion together with the astrophysical limits 
on the evolution of the red giants ruled out the existence of the visible axion. 

1.3.2.2 Invisible axion models 

After the dismiss of the visible axion, a new invisible axion model was postulated. 
Indeed, the breaking scale of the PQ symmetry f a is an arbitrary value that could 
be larger than the electroweak scale. There are mainly two different invisible ax¬ 
ion models, the KSVZ (Kim-Shifman-Vainshtein-Zakharov) and the DFSZ (Dine- 
Fischles-Srednicki-Zhitnitskii), in both of them the axion is light and interact weakly 
with matter. 



Figure 1.1: Axion-photon coupling constant as a function of the axion mass for invisible 
axion models. Three lines for the KSVZ model are shown: E/N = 6 (solid blue line); 
E/N = 0 (dotted green line) and E/N = 2 (dashed black line). Also the DFSZ model 
with E/N = 8/3 (dashed red line) has been drawn. 


KSVZ model 

In this model [11, 12] the interactions of axions with matter only occurs via the 
axion-gluon coupling and ordinary quarks and fermions do not interact with axions. 
A new heavy quark Q has to be introduced, being the unique particle which carries 
PQ charge in the KSVZ model. However, the axion would still interact with the 
light quarks due to the color anomaly. The axion-photon coupling is given by 
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equation 1.17 in which E/N is a model dependent parameter, in general E/N in 
the KSVZ model is given by: 


§ = 6 Qheavy 3 ) 

where Qheavy is the electric charge of the new heavy quark Q and can take values 

of Qheavy = 2/3, —1/3,1,0. The most common values used in the literature are 
E/N = 2, E/N = 6 and the standard KSVZ model with E/N = 0, that are shown 
in figure 1.1. 

DFSZ model 

In contrast with the KSVZ model, in the DSFZ model [13] the fundamental fermions 
carry PQ charge and no exotic quark is needed. The disadvantage of this model is 
that a fine tuning is necessary in order to obtain a breaking scale larger than the 
electroweak scale. On the other hand, the advantage of this model is that it can 
be easily introduced in Grand Unification Theories (GUT). In this case the model 
dependent parameter takes the value E/N = 8/3 (see figure 1.1). 

1.3.3 Summary of relevant axion constraints 

The axion mass is initially arbitrary since the strong CP problem is solved for any 
value of the breaking scale f a . Due to its implications in astrophysics, cosmology and 
particle physics, the axion has been thoroughly studied and the most relevant con¬ 
straints come from astrophysical considerations. Indeed, axions could be produced 
in hot and dense environments like stars, globular clusters and white dwarfs. 

1.3.3.1 Solar model 

Axions could be produced in the solar plasma via Primakoff effect and electron 
Bremsstrahlung processes, being an additional energy loss channel that could reduce 
the lifetime of the star. The most relevant constraint comes from the fact that 
the emission of axions from the Sun implies an increment of the nuclear burning 
and consequently, of the solar temperature distribution, that leads an increase in 
the 8 B neutrino flux [15]. With these considerations the axion-photon coupling 
is constrained to g ai < 7 X 10” 10 GeV . Also, the axion-electron coupling is 
constrained to g ae < 2.5 x 1CU 11 GeV -1 by the Bremsstrahlung collisions. 

1.3.3.2 Globular clusters 

A globular cluster is a spherical collection of ~ 10 6 of densely packed stars that were 
formed about the same time. Two different types of stars in globular clusters are 
especially interesting for the axion bounds, the horizontal branch (HB) stars and 
the red giant branch (RGB). The RGB stars have degenerated helium burning core 
and a hydrogen burning shell. When the RGB helium core becomes hot and dense 
enough it rapidly increases the rate of fusion and becomes a HB star. Thus the 
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axion production via Primakoff effect should be larger in HB stars in comparison 
with RGBs, adding a new energy loss channel which is negligible for RGB stars. 
So the population of RGB stars should be relatively larger than the population of 
HB stars inside the globular clusters. By computing the observed population of HB 
and RGB stars [9], the axion-photon coupling is constrained to g ai < 10 -10 GeV , 
being the most restrictive astrophysical limit on the coupling constant. 

Also the RGB stars are supposed to have larger brightness if the degenerate helium 
core losses too much energy by emission of axions or neutrinos. In this case the axion 
emission would be dominated by electron-Bremsstrahlung processes. Observational 
measurements in the globular cluster M5 [16] lead to a limit on the axion-electron 
coupling of g ae < 4.7 x 1CT 13 GeV -1 . 

1.3.3.3 White dwarf cooling 

White dwarfs (WD) are a remnant of initial low massive stars, made of a degener¬ 
ate carbon-oxygen core and a helium burning shell. A WD is very hot when it is 
formed, but since it has no energy source it gradually cools down, by first the emis¬ 
sion of neutrinos and later by photons from the surface. Also, WD could increase its 
cooling speed by the emission of axions generated mainly by axion-Bremsstrahlung 
processes. This axion emission could be constrained by the comparison of the ob¬ 
served cooling speed with WD models. One special case is the study of the cooling 
of the ZZ Ceti stars [17] in which the most restrictive limit on the axion-electron 
coupling g ae < 1.3 x 10 -13 GeV -1 has been obtained. 

1.3.3.4 Supernova 1987 A 

A supernova (SN) is a stellar explosion that expels a big quantity of stellar material 
with a great force. There are different types of SN, the most interesting one is the 
type II SN. When a star is massive enough to burning its carbon-oxygen core (in 
contrast with white dwarfs), it rapidly consumes the carbon core producing heavier 
nuclei and increasing the core size, at some point the core exceeds the Chandrasekhar 
limit 2 [18] and the star collapses. A cataclysmic implosion takes place within seconds 
generating a huge quantity of radiation and a burst of neutrinos. 

Axions could be emitted from the core of a type-II SN via axion-nucleon 
Bremsstrahlung shortening the neutrino burst. For a very small axion-nucleon cou¬ 
pling g a N , the duration of the burst is not affected, but by increasing the coupling, 
the burst duration becomes shorter until it reaches its minimum (free streaming 
regime in figure 1.2). Also, if g a N increases further, axions are not able to escape 
(trapped regime in figure 1.2) and the burst duration becomes unaffected. 

The SN 1987A gave the first chance to observe directly the neutrino emission from 

2 The Chandrasekhar limit is the maximum mass that can be supported by electron degeneracy 
pressure without suffering a gravitational collapse 
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Figure 1.2: Relative duration of a type II supernova neutrino burst as a function of the 
axion-nucleon coupling g a N- The free streaming region represent axions emitted from 
the core, while in the trapping regime axions arc not able to escape. The dotted line is a 
continuation to guide the eye. Plot taken from [17]. 

an extragalactic source, it was measured by the neutrino experiments IMB 3 4 and 
Kamiokande II 1 operating at this time. A total of 20 neutrinos u e were detected with 
a burst duration of ~ 10 s which fits the expectations without an extra axion cooling. 
These measurements lead to an excluded range on the axion-nucleon coupling of 
3 x 10 -10 < g a N ^3x 10- 7 [19]. However, the low statistics accumulated from the 
neutrino burst and some incompatibilities in the measurements of the two different 
experiments create a big uncertainty in this excluded range. 

1.3.4 Axions as a Dark Matter candidate 

There are several observational evidences of the presence of Dark Matter (DM) in 
the Universe, like the rotation curves of the spiral galaxies, gravitational lensing 
and fluctuations in the Cosmic Microwave Background (CMB). According with the 
latest results of the PLANCK experiment [20] the Universe is composed of a 4.2% of 
baryonic matter, a 27.2% of DM and a 68.6% of Dark Energy (DE); in other words, 
the most part of the Universe has an unknown nature. The DE takes account of the 
accelerated expansion of the Universe, but his nature is not well understood, it can 
be explained as a cosmological constant, vacuum energy or a quintessence held. 

The DM would be composed by electrically neutral particles which interact weakly 
with matter. Popular candidates to DM are the Weakly Interacting Massive Parti- 

3 Irvine-Michigan-Brookhaven detector 

4 Kamioka Nucleon Decay Experiment 
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cles (WIMPs) which appear in supersymmetry models and are searched in several 
underground experiments. The axion is also an attractive DM candidate, axions 
could be thermally produced in the primordial plasma, but being light particles, 
they would contribute to the hot DM as neutrinos do. This production mechanism 
is more efficient for larger masses, so cosmological observations constrain the amount 
of hot DM and fix an upper bound in the axion mass of m a < 0.9 eV [14]. Other 
two related mechanisms could have produced non-relativistic axions in an early Uni¬ 
verse: the vacuum realignment and the decay of the topological defects which make 
the axion an attractive cold DM candidate. 

In the early Universe when the temperature dropped below the PQ breaking scale 
/ a , the axion field appeared with random initial conditions. This means that the 
initial values of the axion field were different in causally disconnected regions. Later 
on, during the QCD phase transition the axion potential rises and the axion acquire 
mass, then the axion held oscillates towards its minimum solving the CP problem 
dynamically. These oscillations take account of the non-relativistic axions, this is the 
so-called vacuum realignment [21] mechanism in which the population of the axions 
depends on the unknown initial value of the held. Moreover, discrete domains with 
topological defects like strings and walls quickly decay producing a large amount of 
cold DM axions. This is the decay of topological defects [22] mechanism. The great 
uncertainly of the axion density depends on whether the inhation happens after or 
before the PQ symmetry breaking. 

In the case that inhation happens after the PQ symmetry breaking, the axion held is 
homogenized by inhation. In this scenario the topological defects are dissipated and 
only the vacuum realignment contributes to the cold DM axion density. Contrary to 
the thermal production, this mechanism is more efficient for smaller axion masses. 
Since the population of the axions depends on the unknown initial value of the held, 
by tuning the initial value of the held, there is a wide range of axions masses which 
gives the proper density of CDM, this is the so-called anthropic window. 

If the inhation occurs before the PQ symmetry breaking, the initial conditions of 
the axion held are randomly distributed in causally disconnected regions and an 
estimation of the vacuum realignment production can be done. In this scenario the 
contribution of the decay of the topological defects has to be taken into account, 
but its calculations are not clear, it could be of the same order of the contribution 
of the vacuum realignment or considerably larger. 

In summary, the preferred mass range for axion cold DM is 1CU 5 < m a < 10 -3 eV, 
also referred as classic axion window [22], Smaller axion masses are still possible in 
fine-tuned models and also axions with masses above the classic window can be a 
subdominant component of cold DM and could coexist with other exotic particles 
like WIMPs. Moreover, axion like particles in a wide range of g a -y — m a space can 
solve the DM problem. A summary of the astrophysical and cosmological constrains 
on the Peccei-Quinn axion in combination with experimental searches are shown in 
figure 1.3. 
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Figure 1.3: Astrophysical and cosmological bounds in combination with experimental 
searches for axions. Dark blue labels represent experimental search ranges of CAST and 
ADMX experiments, while red labels correspond to future prospects. Telescope and lab¬ 
oratory searches (light blue labels) are also drawn. Orange labels indicate strong model 
dependent limits while the green ones are model independent. The limit imposed for the 
SN 1987A burst has been drawn with a different style due to the uncertainties of the mea¬ 
surements. 

1.4 Axion like particles (ALPs) 

Although axions are the best motivated particles, there is also the category of Ax¬ 
ion Like Particles (ALPs) or more generically Weakly Interacting Slim Particles 
(WISPs). ALPs share the same phenomenology of the axion, being light particles 
that weakly couples to two photons. The interaction can be written as: 

£<h = ~ -^9<h < t >F ^F >lu (1-24) 

where 0 represents the ALP held and g its coupling constant. In contrast with the 
PQ axion, ALPs are not motivated by the strong CP problem and g^ and are 
independent parameters. ALPs may arise as pseudo Nambu-Goldstone bosons from 
extensions of the SM in which a new symmetry is broken at a high energy scale. 
ALPs also appear in string theory as the axion does. ALPs have an important role 
in the context of the low energy frontier of particle physics [23]. 

Beyond ALPs, in the category of WISP, there are some particles like hidden photons. 
They appear in extensions of the SM in which hidden sectors are included, these 
sectors do not directly interact via the gauge boson forces of the SM, being the 
interaction through the interchange of very heavy particles. Hidden photons can 
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oscillate with standard photons via a kinetic mixing and thus they show the same 
behavior of the axions and its coupling with photons, in this case a magnetic held 
is not necessary for the conversion. 

ALPs and hidden photons separately could provide all the amount of DM. The non 
thermal mechanisms of production are the vacuum realignment and the decay of 
the topological defects, the same mechanisms described for the axion, introduced in 
section 1.3.4. The contribution of the WISP to the total amount of DM has been 
recently studied [24] and a wide range of the parameter space (g ( p- r — rri^) can contain 
models with the right density of DM. 


1.5 Astrophysical hints for axions and ALPs 

Several astrophysical and astronomical observations could be interpreted as hints 
of axions and ALPs. Like the excessive transparency of the Universe to very high 
energy photons (VHE) and the anomalous cooling rate of white dwarfs (WD). These 
hints will be briefly described below. 

1.5.1 VHE transparency 

Very high energy photons (VHE) with energies of ~ 1 TeV traveling through the in- 
tergalactic medium have a non negligible probability to interact with the background 
photons permeating the Universe. VHEs could interact with the extragalactic back¬ 
ground light (EBL) via pair production e _ e + . Thus, the Universe is expected to be 
relatively opaque to distant VHE sources like active galactic nuclei (AGN). 

The EBL density is estimated by the measurements of the spectra from distant 
blazars by HESS [25] and Fermi [26], in good agreement with theoretical models. 
However, there are different observations that indicate a high transparency of the 
Universe to VHE photons, even in the most favored EBL models. Actually, experi¬ 
ments based on atmospheric Cherenkov telescopes like HESS [27] and MAGIC [28] 
have measured VHE photons with arrival directions clearly correlated with AGNs 
with an spectra which requires a low EBL density or an anomalous spectra in the 
origin. 

The high transparency of the Universe to the VHE photons could be explained 
by ALP-photon oscillations. VHE could be transformed into ALPs in the local 
magnetic fields at the origin, in the intergalactic magnetic fields and in the Milky 
Way. These ALPs could travel through the intergalactic medium without inter¬ 
acting. These oscillations lead to a more transparent Universe to the VHE pho¬ 
tons. This kind of ALPs [29] requires a low mass < 1CU' eV in order to pre¬ 
serve the coherence over large magnetic fields and relative large coupling constants 
g^ ~ 10 -12 — 10 _1 ° GeV -1 . Although these parameters are far away from 

standard axion models, ALP models in which g and are uncorrelated could 
take account of these oscillations. 
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1.5.2 The WD cooling 

The WD cooling was introduced in section 1.3.3.3 and a limit on the axion-electron 
coupling could be extracted from the study of the WD evolution. However, recent 
studies indicate an anomalous cooling speed which could be a hint of an axion 
emission. Indeed, the luminosity function 5 of the WDs is predicted with a great 
accuracy by theoretical models. An axion emission could be a source of extra cooling 
and may suppress the luminosity function of the WDs at some point. 

Recent works based on well studied WD models and the measurements of different 
WDs luminosity population point to the possibility that a small amount of energy 
loss coming from axions is feasible. Axions with coupling constants of g ae ~ 2 X 
ICR 13 [30] seem to fit the experimental measurements. Also the evolution of the 
pulsating period of certain WDs can provide a direct measurement of their cooling. 
The different pulsing of ZZ Ceti WDs, G117-B15A [31] and R548 [32] have been 
computed, obtaining a similar result of g ae ~ 5 X 1CU 13 . These results create some 
tension with another limits coming from astrophysical considerations. In any case, 
the WD models seem to improve by adding an extra cooling coming from axions 
with masses above m a ~ 3 — 20 meV and can be interpreted as a hint of axions. 
Also, ALPs models could fit to the observed cooling rate of WDs. 

1.6 Axion and ALPs searches 

As it was presented in previous sections, axions and ALPs could be generated in an 
early Universe and also can be produced in hot and dense environments like stars. 
So that the different search strategies can be divided depending on the source: halo- 
scopes and telescopes are looking for cold DM axions permeating the galactic halo, 
helioscopes that are searching Solar axions and photon regeneration experiments in 
which axions could be produced and detected in the laboratory. All these strategies 
are based on the Primakoff effect, in which axions could be converted into photons in 
strong electromagnetic fields. The most sensitive experiments for axions and ALPs 
searches are shown in figure 1.4. Also, the most relevant constrains and hints are 
drawn. 

1.6.1 Galactic halo searches 

Galactic halo searches could be divided in two different kind of experiments: halo- 
scopes which are based on resonant cavities in a magnetic held and telescopes that 
are looking for the decay of the axion to two photons. 

1.6.1.1 Haloscopes 

The haloscope technique was proposed by Sikivie [33] in 1983. As it was described 
in section 1.3.4, axions are attractive cold DM candidates. These relic axions could 

5 Number of stars per luminosity interval 
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m a (eV) 


Figure 1.4: Axion/ALP parameter space, the most sensitive experiments are shown: lab¬ 
oratory experiments (ALPS), helioscopes (CAST), haloscopes (ADMX) and telescopes. 
Also the most significant constrains from horizontal branch (HB) stars and hot dark matter 
(HDM) are drawn. The yellow band represents the most favored axion models and the 
green line corresponds to the KSVZ model ( E/N = 0). The orange band is the "classic 
axion window" with the most interesting models of cold dark matter axions. The region 
inside the blue dashed line indicates the preferred region from the WD cooling. For more 
generic ALPs, the parameter space below the dashed red line may contain cold dark matter 
ALPs. The dashed gray line on the left represents the most favored ALPs models in the 
context of the transparency of the Universe to VHE photons, while the maroon region has 
been excluded by the HESS experiment. Plot taken from [135]. 


be detected by resonant cavities inside strong magnetic fields, this is the so-called 
haloscope technique. Relic axions could be converted into photons via Primakoff 
effect in a magnetic field. Being these axions non relativistic, the converted photons 
would be monochromatic with energy around m a . The resonant process in the cavity 
is enhanced when its frequency matches m a . This technique only allows to scan a 











1.6. Axion and ALPs searches 


17 


narrow range in the axion mass. However, the resonant cavity is tunable and the 
data taking could be performed by scanning narrow ranges of m a . 
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Figure 1.5: Top: First ADMX results in contrast with RBF and UF searches, taken 
from [36]. Bottom: Latest ADMX results with a 90% of confidence level; plot taken 
from [37]. 


The first haloscopes were developed by Rochester-Brookhaven-Fermilab (RBF) and 
the University of Florida (UF) [35]. They demonstrate the technique in a wide range 
of axion masses 5 < m a < 16 peV but with a poor sensitivity (see figure 1.5 top). 
Later on, a second generation haloscope was developed, the Axion Dark Matter 
experiment (ADMX), being the first haloscope with enough sensitivity to reach the 
most favored axion models. In a first stage ADMX scanned an axion mass region 
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from 1.9 < m a < 3.3 ^eV [36], in a second phase the microwave receivers was 
substituted by SQUID (Superconducting QUantum Interference Device) in order to 
improve the sensitivity and scanning a mass range from 3.3 < m a < 3.53 /reV [37] 
(see figure 1.5 bottom). However, these results assume that axions are the main 
component of cold DM in the Universe. 

1.6.1.2 Telescopes 

Although the decay of the cold DM axions to two-photons is expected to be ex¬ 
tremely odd, it could be detected in thermally produced axions in the eV regime. 
The signal will be a closely monochromatic line emitted from galaxies and could be 
observed by telescopes used in astrophysics. Different works have excluded axion 
masses of 3 < m a < 8 eV [38, 39]. 

1.6.2 Laboratory searches 

In laboratory experiments, photons could be converted into axions or ALPs inside 
a strong magnetic field via Primakoff effect. The advantage of this technique is that 
it does not rely in astrophysical or cosmological assumptions about the origin of the 
axions or ALPs. On the other hand, the sensitivity of this kind of experiments is 
usually lower than other search strategies. There are mainly two different techniques 
in laboratory searches: the Light Shining through Wall (LSW) and the polarization 
experiments. 

1.6.2.1 LSW searches 

LSW or photon regeneration experiments [40] use high intensity laser beams inside 
a strong transverse magnetic field in order to produce axions or ALPS, in which an 
opaque wall is placed to block the photons. On the other side of the wall another 
magnetic field allows the axions or ALPs, which could easily pass the wall, to be 
transformed again into photons that can be detected. 

Several experiments have excluded significant regions in the ALP parameter space: 
the LIPSS 6 [41] experiment at the Jefferson Laboratory, the BMV 7 [42] collaboration, 
the GannneV 8 9 [43] experiment at Fermilab, the OSQAR°[44] experiment working at 
CERN and finally the most sensitive LSW so far, the ALPS 10 [45] experiment at 
DESY. The results of all these experiments are shown in figure 1.6. Although the 
experiments related before use laser beams, the same principle of detection could be 
adapted to the microwave regime. This is the case of the CROWS 11 experiment [46] 
at CERN. 

6 Light Pseudoscalar or Scalar Search 

1 Birefringence Magnetic du Vide 

8 Gamma to mili-eV particle search 

9 Optical Search for QED vacuum magnetic birefringence 

10 Axion-Like Particle Search 

11 CERN Resonant Weakly Interacting sub-eV Particle Search 
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Figure 1.6: Exclusion limits for axion like particles in different LSW and polarization 
experiments. The dashed and dotted lines correspond to the PVLAS experiment. Plot 
taken from [45]. 


1.6.2.2 Polarization experiments 

In contrast to LSW experiments, it is possible to obtain an axion signal in the laser 
beam itself. When a polarized light passes through a transverse magnetic held it 
could generate axions or ALPs and thus, it would induce a detectable small rotation 
and ellipticity in the laser beam. 

Experiments like the BFRT 12 [47] at the Brookhaven National Laboratory and the 
PVLAS 13 experiment at the the Legnaro National Laboratory exploited this tech¬ 
nique and excluded some regions of the parameter space (see figure 1.6). Moreover, 
the PVLAS experiment reported a signal, but it was ruled out later on since it was 
interpreted as an artifact of an instrument [49]. 


1.6.3 Solar axion searches 

The Sun is expected to be a powerful axion source as it was presented in sec¬ 
tion 1.3.3.1. Solar axions could be detected in the laboratory via inverse Primakoff 
effect. Two different techniques for the detection of solar axions have been devel¬ 
oped: helioscopes and Bragg scattering experiments. 

12 Brookhaven Fermilab Rutherford Trieste collaboration 
1J Polarizzazione del Vuoto con LASer 
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1.6.3.1 Helioscopes 

The helioscope technique was proposed by Sikivie [33] in 1983. The principle of 
detection is to reconvert the solar axions into photons inside an intense transversal 
magnetic held. Solar axions may have energies in the keV regime, thus the expected 
signal are X-rays in the detectors located at the edges of the magnet. Helioscopes 
can cover axion masses in a wide range of the parameter space (up to ~eV) because 
contrary to haloscopes, the signal is independent of the axion mass. The helioscope 
technique will be revisited in more detail in section 1.7 



m a (eV) 

Figure 1.7: Excluded region for different experiments based on solar axion searches at 
95% of C.L. Dashed lines correspond to Bragg scattering experiments; from top to bot¬ 
tom: SOLAX & COSME (black), CDMS (brown), EDELWEISS (blue) and DAMA (red). 
The solid lines represent helioscopes searches: Lazarus et al. (blue), Tokyo Axion Helio¬ 
scope (green) and CAST (red). Typical axion models are contained in the orange band, 
the solid gray lines correspond to the KSVZ model ( E/N = 0) and the DFSZ model 
(E/N = 8/3). 


Different experiments based on the helioscope technique have been developed, the 
first one was performed by Lazarus et al. [50] validating the technique. Later on, 
the Tokyo Axion Helioscope (SUMICO) [51] applied the same principle of detection 
increasing the sensitivity. Finally, the CERN Axion Solar Telescope (CAST) started 
in 2003 being the most sensitive helioscope until now. The excluded regions of these 
experiments together with the Bragg scattering searches are shown in figure 1.7. 
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The CAST experiment, which is the main topic of this work, will be detailed in 
chapter 2. 

1.6.3.2 Bragg scattering searches 

Crystalline detectors can also be used for solar axion searches. The axion-photon 
conversion could take place in the Coulomb held of the nuclei in the crystal lattice. 
The conversion occurs when the angle of the incident axions fulfills the Bragg con¬ 
dition with the plane of the crystal. Thus, the expected signal is a characteristic 
Bragg pattern in the detector. 

Even though there are not many dedicated experiments, a great number of under¬ 
ground experiments developed for WIMP searches have been looking for these pat¬ 
terns. Such as SOLAX 14 [52], COSME 15 [53], DAMA 16 [54], CDMS 17 [55] and more 
recently EDELWEISS 18 [56]. However, the sensitivity of this technique cannot com¬ 
pete with dedicated helioscope experiments, as is shown in figure 1.7. 


1.7 Solar axions and the helioscope technique 

Axions could be generated in the core of the Sun, being the most powerful axion 
source closest to the Earth. The expected production of axions in the well known 
Solar Model will be described in the next section. The detection of solar axions 
comes through its conversion into photons inside strong magnetic fields in which the 
probability of the axion-photon conversion has an important role and will be further 
detailed. 

1.7.1 The Solar Model and the axion production 

The solar axion flux can be extracted from the well established Standard Solar 
Model. The most relevant channel in the axion emission comes from the Primakoff 
conversion which is dominant in hadronic axion models such as the KSVZ model. 
In the case of non-hadronic axions (DFSZ model) with a tree level interaction with 
electrons, other processes have to be taken into account: Bremsstrahlung, Compton 
and axion recombination. The contribution of these channels could be considerably 
larger than the Primakoff emission. However, in the helioscope technique only the 
Primakoff contribution is usually taken into account because it is more reliable to 
suppose the same processes involving the generation and detection of the axions. 
Also the astrophysical constraints on g ae are more restrictive than the limits that in 
principle could be reached with helioscopes. 

14 SOLar AXion search in Argentina 

15 Germanium detector located in the Canfranc Underground Laboratory 
lc DArk MAtter searches, operating at Gran Sasso 
17 Cryogenic Dark Matter Search, located in Soudan 
18 Experience pour DEtecter Les Wimps En Site Souterrain, at Modane 
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1.7.1.1 Hadronic axion emission 

Axions could be generated in the core of the Sun inside the strong electric fields 
from the charged particles of the solar plasma by the process 7 + Ze —> Ze + a. 
Thus, the energies of the converted axions correspond to the blackbody radiation of 
the core of the Sun (keV regime). 

The Primakoff process in non-relativistic conditions is relevant when the mass of 
electrons and nuclei are larger than the energies of the ambient photons. Follow¬ 
ing [9], neglecting recoil effects, the differential cross section is given by: 


do- a -> 7 _ 9 2 ai Z 2 a |p 7 x p a 
dfl 87T q 4 


(1.25) 


here Z is the charge (Ze) of the involved particle, a is the fine structure constant, p 7 
and p a are the momentum of the photon and the axion respectively and q = p 7 — p a 
is the momentum transferred. 

The cut-off of the Coulomb interaction in vacuum for m a 7^ 0 is given by the 


minimum momentum transfer q m %n = m^/2Fn for m a E a and the total cross 
section can be written as: 



(1.26) 


In a plasma, the Coulomb potential is cut-off due to screening effects and thus the 
differential cross section is modified by an additional factor in equation 1.25. 


dcjg^ _ g 2 a 1 Z 2 oi |p 7 X p a \ 2 \q 2 \ 
dtt 87r q 4 k 2 +\q 2 


(1.27) 


In a non degenerate medium, the screening scale k s is given by the Debye-Hiickel 
formula 



(1.28) 


where T is the temperature of the solar plasma, ng is the baryon density and Y e 


and Yj are respectively the fraction of electrons and different nuclear species j per 
baryon. Using equations 1.27 and 1.28, the total scattering cross section can be 
calculated. Summing over all the targets, the transition rate r 7 _> a of a photon with 
energy into an axion of the same energy E a = E^ is obtained. 
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here, the effective photon mass and axion mass are negligible in comparison with 
E a . The differential axion flux at Earth can be calculated by the convolution of the 
transition rate and the blackbody emission of the Sun. 
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(1.30) 


being cZq = 1.5 X 10 s km the average distance from Earth to Sun and Rq = 
7.0 x 10 5 km the solar radius. Using the Standard Solar Model from [57] an analytical 
approximation of the axion differential spectrum can be derived [59]. 


= 6.02 x 10 10 ( - 9 -°E - (E a / keV) 2 ' 481 e --E»/i-205keV [ cm - 2 s -1 keV -:1 l 

V10“ 10 GeV -1 / v ' ’ V J 

(1.31) 



Figure 1.8: Differential solar axion flux at Earth for hadronic axion models. The calcu¬ 
lations have been made using a coupling constant of g ai = 10“ 10 GeV -1 . Two different 
lines are drawn, corresponding to the solar model from 1982 [58] (dotted black line) and 
an updated solar model from 2004 [57] (solid red line). 

The differential solar axion flux is proportional to the coupling constant with an 
average energy of < E a >= 4.2 keV and a maximum around 3 keV (see figure 1.8). 
Also the total axion luminosity of the Sun has been estimated: 

L a = 1.85 x 10 -3 ( -- T ^) L(7) (1.32) 

\10 -10 GeV -1 / ° v ' 


where Lq is the solar luminosity. 
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1.7.1.2 Non-hadronic axion emission 


In non hadronic axion models with an axion-electron coupling at tree level (DFSZ 
model), the solar axion emission is dominated by the axion-electron coupling. Fol¬ 
lowing the work of Redondo [60], the most important production comes from the 
interactions: electron-ion Bremsstrahlung (e + I — > e + I + a), electron-electron 
Bremsstrahlung (e + e— > e + e + a), Compton scattering (7 + e—>■ 7 + e + a), axio- 
recombination (e + / — > I + a) and axio-deexcitation (I* —> / + a). The Feynman 
diagram of these processes together with the Primakoff process are shown in figure 
1.9. 
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Figure 1.9: Feynman diagrams of the most relevant solar axion production channels. Plot 
taken from [60]. 


The Bremsstrahlung or free-free (ff) electron transitions on hydrogen and helium 
nuclei dominates solar flux together with the electron-electron (ee) Bremsstrahlung. 
The axio-recombination, also know as free-bound (fb) electron transitions, of metals 
(mainly O, Ne, Si, S and Fe) has a significant contribution. The axio-deexcitation, or 
bound-bound (bb) electron transitions, is dominated by Lyman transitions (mostly 
Ly-a) being significant in the case of iron (~ 6.4 keV). The Compton (C) production 
process has a poor contribution. 

The emission rate of solar axions for all the processes described before, can be 
expressed as: 


r a = r l f + r/ 6 + rf + rf + rf (1.33) 

where Fis the emission rate of the sum over all the processes. 

Electron-ion Bremsstrahlung, axio-recombination and axio-deexcitation 

In processes in which a photon is emitted when a electron makes an atomic transition 
ei — > ef. The spin-averaged matrix element of emitting an axion of energy E a is 
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proportional to the analogous matrix element of emitting a photon of the same 
energy. Thus in this case the emission rates of photons and axions are proportional 
and can be formulated as: 


Y^\E a ) _ 1 g\ e El 
Y^\E a ) 2 e 2 m 2 


i = ff, fb, bb 


(1.34) 


In the case of the electron-ion Bremsstrahlung (ff), the photon production rate is 
given by: 
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where Z is the charge of the involved nuclei and n z its number density, n e is the 
electron density, F(E a , k s ) take account of the screening effects and k s is the screen¬ 
ing scale introduced in equation 1.28. Therefore the transition rate for electron- 
Bremsstrahlung processes can be derived: 
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Electron-electron Bremsstrahlung 

In this case the expression 1.34 is not valid because the photon emission is highly 
suppressed in comparison with electron-ion processes. The production rate involving 
this process has been calculated in [9] and is given by: 


pP,EEO ( p \ _ 2 

' a \E a ) — CX 9, 


An 


n% 


~—E. 


3 \fTmJ 2 E a 


/ T F{E a ,V2k s /^/2m e T) (1.37) 


Compton scattering 

In the non-relativistic limit, the cross section of the axion production in Compton¬ 
like processes is given by [19]: 
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(1.38) 


Thus the transition rate can be expressed as: 


r a' C {E a ) 
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(1.39) 


here the ratio between the photon and the axion production rate is slightly different 
from equation 1.34. 
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r a P ’ C (E a ) _ g\ e El 
n y ,C (E a ) e 2 m 2 


(1.40) 


All the contributions to the axion emission presented before can be expressed in 
terms of the photon absorption coefficient with the exception of the electron-electron 
Bremsstrahlung. Thus, the total axion flux in earth can be calculated analogous to 
equation 1.30 
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where is the emission rate of all the processes presented in 1.33, this term depends 
on the local parameters of the solar plasma. By integrating equation 1.41 using the 
solar model from [61] and the different contributions presented before, the axion 
flux at Earth has been extracted [60] and is shown in figure 1.10. 



to [keV] 


Figure 1.10: Differential solar axion flux at Earth from non-hadronic processes, assuming 
a coupling constant of g ae = Kb 13 . Different contributions have been drawn: atomic re¬ 
combination plus atomic deexcitation (FB + BB), Bremsstrahlung process (FF) and Comp¬ 
ton. The solid blue line is the Primakoff emission with g ai = 10 -12 GeV -1 which has 
been multiplied by a factor 50 to make it visible. Plot taken from [60] 
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Contrary to the hadronic axion emission, the total axion flux cannot be approx¬ 
imated to an analytical formula, due to the narrow lines from the bound-bound 
and free-bound transitions. Also, the maximum of the spectrum is reached at lower 
energies (~1 keV), in contrast with the Primakoff emission. On the other hand, 
the Bremsstrahlung (B) and Compton (C) processes could be approximated to an 
analytical formula: 
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here, the Bremsstrahlung flux includes electron-ion and electron-electron processes. 

1.7.2 Probability of the axion-photon conversion 

In the helioscope technique, solar axions could be converted into photons inside 
strong magnetic fields by the inverse Primakoff effect. This conversion is only ef¬ 
fective when the magnetic field is transversal to the direction of propagation of the 
incoming axion. The wave equation of a particle propagating along the z-axis in a 
transverse magnetic field B is given in [62] and can be formulated as: 


id z 



(1.44) 


here An is the amplitude of the photon field parallel to the magnetic field, a is 
the amplitude of the axion field, T is the inverse absorption length for the photons 
and wiry is the effective photon mass, which is given by the plasma frequency of 
the medium. A first-order solution using a perturbative approach has been derived 
in [59] and is given by the expression: 
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In general B, T and m 7 are functions of 2 . By considering the case in which the 
magnetic field B and the density are uniform and T and m 7 are constant. The 
conversion probability P a ^ 7 of an axion into a photon inside a magnetic field with 
an effective length z = L, can be derived: 
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P a ^ 7 = | < A||(L)|a(0) > | 2 
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where q is the momentum transfer between the axions and the photons in the 
medium, which is given by: 


2 2 
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(1.47) 


Two different cases will be detailed below, the first one considering the conversion 
inside vacuum and a second case with a buffer gas inside the conversion volume. 


1.7.2.1 Probability of conversion in vacuum 

For the vacuum case, the probability of conversion can be simplified considering the 
absorption term negligible (T ~ 0). Applying this condition to equation 1.46, the 
probability of conversion is reduced to 19 : 
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(1.48) 


Also in this case the effective photon mass could be considered negligible m 7 ~ 0. 
Thus, the momentum transfer q from equation 1.47 can be simplified. 


q = 


ml 

2 E a 


(1.49) 


This yields a condition in which the coherence is maximized for qL < 27T due to the 
sinusoidal term, that can be expressed in terms of axion mass: 


m a < 



(1.50) 


Assuming a magnetic field of 9 T with an effective length about 10 m and an axion 
energy of E a ~ 4.2 keV, the coherence condition is fulfilled for axion masses of 
m a < 0.02 eV (see figure 1.11 top). Also, in the case of full coherence qL <C 2ir, the 
equation 1.48 can be simplified 20 . 

p a ^ 7 = for m,a < ( L51 ) 


19 Note that the relation sin 2 x = 
20 By using the approach sin a; 
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Figure 1.11: Axion to photon probability conversion as a function of the axion mass, the 
calculations are made for axions with an energy of E a = 4.2 keV, inside a transversal 
magnetic field of B = 8.8 T over a length of L = 9.26 m and with a coupling constant of 
g ai = 10“ 10 GeV” 1 . The black line represent the vacuum case, while the red line assumes 
3 He as buffer gas with a density of p = 5 x 10~ 4 g cm -3 , which corresponds to a photon 
mass of m 1 = 0.525 keV. The figure on the top is a general view, while the figure at the 
bottom is a zoom of the buffer gas case. 


1.7.2.2 Probability of conversion in a buffer gas 

In order to be sensitive to higher axion masses the conversion region can be filled 
with a buffer gas. In this case the effective photon mass from equation 1.47 has to 
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be taken into account. 

While in vacuum the photons are massless, an effective mass can be considered when 
the photons go through a medium in which the traveling speed is less than the speed 
of light c. This can be interpreted as an effect of mixing between the photon and 
the quantum excitations of the matter. The effective photon mass is given by the 
plasma frequency of the medium oo p : 

lAnarie 

= u p = \j^r 

here n e is the electron density of the buffer gas, which is related with the gas density 
by the expression: 


(1.52) 


n e = n A^P (1-53) 

Z and W are the corresponding atomic number and atomic weight of the gas, Na 
the Avogadro’s number and p the gas density. 

The coherence is restored when the momentum transfer tends to zero q = 0 (see 
equations 1.46 1.47), when the axion mass and the effective photon mass are equal 
m a = m~, (see figure 1.11 bottom). Also, the absorption T could not be considered 
negligible and has to be taken into account. This technique has been used in helio¬ 
scopes in order to scan higher axion masses by increasing the buffer gas density in 
small steps. 
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2.1 Introduction 

The CERN Axion Solar Telescope (CAST) experiment is located at building SR8 
at CERN and is looking for solar axions since 2003. CAST exploits the helioscope 
technique using a decommissioned LHC 1 dipole magnet [63] (see figure 2.1) that 
provides a magnetic held up to 9 T in which solar axions could be converted into 
photons. The magnet is mounted on a movable platform that allows tracking the Sun 
~1.5 hours two times per day, during sunset and sunrise. The magnet is composed 
by two magnetic bores with X-rays detectors placed at the bore ends. The signal 
of axions would be an excess of counts in the X-rays detectors while the magnet is 
pointing the Sun. 

The CAST research program can be divided in two phases, the first one with vac¬ 
uum inside the magnetic bores and a second phase with buffer gas in which 4 He and 
3 He were used separately in different periods. During Phase I, the X-ray detectors 
installed were a conventional Time Projection Chamber (TPC), a Micromegas (MI¬ 
CRO MEsh GAseous Structure) detector and a Charge Coupled Device (CCD) in 
the focal plane of an X-ray focusing device. During Phase II, the TPC was replaced 

1 Large Hadron Collider 
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Figure 2.1: Experimental set-up of the CAST experiment. The dipole magnet (blue) is 
mounted on a movable platform (green) that allows the horizontal and vertical movement. 


by two Micromegas detectors. This chapter describes the main features of the CAST 
experiment together with a review of the X-ray detectors that have been working 
at CAST. 


2.2 Technical description 


The CAST experiment is composed by several additional systems required for the 
data taking operation. The magnet operates at a nominal temperature of 1.8° K 
and a cryogenic cooling system is required. Also, a vacuum system is installed 
around the magnet bores in order to increase the transparency of the X-rays from 
the axion-photon conversion and to isolate the magnet from the environment. 

The magnet is mounted on a movable platform that allows a vertical movement of 
±8° and 80° in horizontal (see figure 2.2). The movement of the magnet is con¬ 
trolled by a tracking system program which allows to point the Sun automatically. 
Moreover, a gas system has been installed to fill the magnet bores with a buffer 
gas ( 4 He and 3 He) in small steps in order to restore the sensitivity to higher axion 
masses. The monitoring of all the systems is performed by the slow control. The 
different systems of the CAST experiment will be detailed below. 
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Figure 2.2: Schematic drawing of the CAST experiment. The different parts of the magnet 
and the movable platform are labeled. 

2.2.1 The CAST magnet and cryogenics 

The CAST magnet is one of the first prototypes designed for the LHC [64] (see 
figure 2.3). In contrast to the final bend design of the LHC, the CAST magnet 
has two straight beam pipes of 9.26 m that allow to exploit the total bore aperture 
of 4.30 cm of diameter (14.52 cm 2 cross-section). The magnet is made of a super¬ 
conducting Niobium-Titanium (NbTi) alloy working at a nominal temperature of 
1.8° K. During normal operation the magnet is loaded at 13 kA, providing a mag¬ 
netic field of 8.8 T. The magnetic field is perpendicular to the beam pipes as it is 
shown in figure 2.4, which allows the axion to photon conversion in the direction of 
the beam pipes. 

The cryogenic system is supplied with liquid helium [65], using the cryogenics from 
the former LEP 2 collider and its experiment DELPHI ’. The magnet is supplied by 
the Magnet Feed Box (MFB) placed on the top of the magnet, above the rotary pivot 
(see figure 2.2). The MFB provides liquid helium and high current to the magnet 
via flexible cables. The cooling of the magnet is performed in different phases, in 
the final stage the liquid helium becomes superfluid and circulates in a cycle over 
the magnet from the MFB to the Magnet Return Box (MRB). 

During magnet operation an abnormal termination can occur, a part of the su¬ 
perconducting coils can enter in a resistive state and raise the temperature of the 
surrounding region. This phenomenon is called quench and can be destructive for 
the magnet. In order to prevent a it, a quench protection system has been installed 
at CAST. When the protection system is triggered, the quench heaters are acti- 
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Figure 2.3: Standard cross section of a LHC dipole magnet. 



Figure 2.4: Orientation of the magnetic field in a dipole magnet. 

vated, generating a controlled increase of the magnet temperature that is uniformly 
distributed. The rise of the temperature provokes an overpressure in the system and 
the helium is purged in order to prevent damages in the magnet. 
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2.2.2 The tracking system 

The tracking system is crucial for the required precision of CAST while is point¬ 
ing the Sun. For this purpose, an accurate hardware and software were designed 
[66]. Also, two complementary alignment tests are performed regularly, the grid 
measurements and the Sun filming. All these features will be described below. 


2.2.2.1 Hardware 

The CAST magnet is mounted on a movable platform that allows a vertical move¬ 
ment of ±8° (polar angle) and from 46° to 133° in horizontal (azimuth angle). The 
rotary pivot is placed underneath the MFB in order to ensure the stability of the 
cryogenics while the magnet is moving. The MRB side of the magnet is supported 
by a trolley with four wheels mounted on rails, which allows the horizontal move¬ 
ment (see figure 2.2). The vertical movement is performed by two jacks spinning 
around two lifting screws mounted on the trolley. 



Figure 2.5: Left: Vertical encoder and lifting jacks. Right: Driving wheel and horizontal 
encoder. 


The horizontal motion is performed by a driving wheel connected to a motor, for 
the vertical movement another motor control the speed of the lifting jacks. Both 
horizontal and vertical motions are monitored by two different motor encoders (see 
figure 2.5). The values of the encoders can be translated into polar and azimuthal 
angles to check the orientation of the magnet. 
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2.2.2.2 Software 

In order to perform the Solar tracking a Labview based software was developed (see 
figure 2.6). The program calculates the position of the Sun using the NOVAS 4 [67] 
software. During the tracking the program communicates with the hardware. It 
calculates the position of the Sun in the next minute, then the motor velocities are 
modified using the position of the magnet from the motor encoders. The expected 
and the real tracking position are monitored and thus the precision of the tracking 
is checked every minute. 
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Figure 2.6: Labview based tracking program of the CAST experiment 


2.2.2.3 Tracking system alignment 

The precision of the solar tracking is essential in the helioscope technique. Two 
different methods are used in order to check the reference position of the magnet 
and the real one: the grid measurements and the Sun filming. 

The Sun filming consists in monitoring the visible light of the Sun through a window 
in the experimental area. These measurements only can be performed in March 
and September, when the Sun is reachable. For this purpose a photo camera is 
aligned with the magnet and takes photos of the Sun during morning trackings. 
Afterwards, a software is used to compare the deviation of the center of the Sun 
using the reference position of the camera. 

4 Naval Observatory Vector Astrometry Software 
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Comparison with September 2007 GRID (2011-2007) 




Figure 2.7: Left: Photo of the Sun during the Sun filming and analysis program that 
calculates the center of the solar disc. Right: 2011 grid results in comparison with the 
2007 grid for the two magnet bores V1 (blue points) and V2 (red points). The green circle 
corresponds to a deviation of 1 arcmin and the red circle represents a 10% of the solar 
radius. The units of the deviation are in mm per 10 m. 


The grid measurements are performed at least once every year and consist in mea¬ 
suring different reference positions of the magnet that are compared with previous 
years. In case of any deviation of the reference position the new grid can be imple¬ 
mented in the software program. 

2.2.3 The vacuum system 

The vacuum system at CAST can be divided in different subsystems depending 
on its purpose (see figure 2.8). The cryostat vacuum isolates the magnet from the 
environmental temperature and the general vacuum line separates the magnet from 
the detectors. At the end of the line four gate valves are installed (VT1, VT2, VT3 
and VT4), one per detector bore end. During normal operation the gate valves are 
open, but they can be closed separately in case of failure. Finally, in the detector 
side of the gate valves there are several vacuum subsystems for the different detector 
lines. 

2.2.4 The gas system 

The sensitivity to higher axion masses can be restored by the addition of a buffer 
gas inside the magnet bores. The CAST experiment has used 4 He and 3 He as buffer 
gases in different periods. In order to scan a wide range of axion masses in small 
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Figure 2.8: Scheme of the CAST magnet vacuum. The cold windows (green) separate the 
buffer gas volume (red) from the general vacuum system (yellow). The cryostat vacuum 
(blue) provides thermal insulation from the environment. 


steps, precise amounts of gas has to be inserted inside the cold bores and thus, an 
accurate gas filling system was installed. 
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The buffer gas has to be confined inside the magnetic field region, for this purpose 
four cold windows were installed at the end of the magnet bores (see figure 2.8). The 
requirements of the cold windows are: high X-ray transmission, a low leak rate of 
helium from the cold bore to the vacuum side and robustness against rapid increases 
of the pressure (e.g. during a quench). 



Figure 2.9: Left: Photo of the cold window placed on the strongback. Right: Drawing of 
the cold window and its different parts. 

The cold windows are made of 14 /mi polypropylene foil glued on a stainless steel grid 
structure also called strongback (see figure 2.9). The cold windows was manufactured 
and tested by the CERN Central Cryogenics Laboratory (Cryolab). 

The gas system was upgraded when the 4 He was replaced by 3 He. Since the 3 He 
is an extremely expensive gas, additional systems were installed to prevent any 
leak [68]. A scheme of the filling system is shown in figure 2.10 in which the dif¬ 
ferent elements are labeled: storage volume, purging system, metering volumes and 
expansion volume. 

The storage volume has a capacity of 963 1, here the gas remains below atmospheric 
pressure in order to preserve the gas in case of leak. The purging system is composed 
by two charcoal traps in order to ensure the purity of the gas. The metering volumes 
MV2 and MV10 have capacities of 1.63 1 and 8.58 1 respectively. The MV10 volume 
is used to fill big amounts of gas while the MV2 volume is used to insert more 
accurately small amounts of gas. For a precise calculation of the amount of gas, 
both metering volumes are inside a thermal bath with a constant temperature of 
36.0° C. 

The gas is sent from the storage volume to the respective metering volume through 
the purging system. Later on, the gas in the metering volume is inserted inside the 
conversion volume of the magnet. The gas system is controlled by a PLC 5 system 
that allows the transfer of the gas to the different volumes. The PLC is also used to 
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Figure 2.10: Scheme of the gas system. 
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perform the recovery process. The PLC is integrated with a SCADA 6 system that 
provides a GUI' which communicates with the pneumatic valves and pumps of the 
gas system. 

In case of quench, the pressure of the conversion volume inside the magnet can 
rise dramatically. In order to prevent a break of the cold windows, the PLC sends a 
signal to open two electrovalves that connect the conversion volume to the expansion 
volume where the buffer gas can flow. The expansion volume has a capacity of 450 1 
and is placed on the top of the magnet. After a quench the recovery process is 
started and the gas is transferred from the expansion volume to the storage volume. 

2.2.5 The slow control 

The slow control is a centralized system designed to monitor and control the main 
parameters of the different subsystems. The CAST slow control is a Labview based 
data acquisition and plotting system. It communicates with multiple NI S 8 data 
acquisition cards connected to different sensors of the CAST experiment. A great 
number of parameters are monitored such as: pressures (detectors, buffer gas and 
vacuum systems), temperatures, magnet movement, valve status and many others. 

All these parameters are measured in real time and can be displayed on the program 
(see figure 2.11). The data is stored every minute and periodically transferred to 
the AFS 9 storage system. Also, different alarms are triggered if some parameters 
are out of a given limit, in this case a SMS 10 is sent to the corresponding responsible 
of the system. 

2.3 The CAST research program 

The CAST experiment started in 2003 being the most sensitive helioscope so far. 
The data taking at CAST was divided in two different phases: a Phase I with 
vacuum and a Phase II with a buffer gas inside the magnet bores. Furthermore, two 
different periods could be distinguished during Phase II: the first one by using 4 He 
as buffer gas and the second with 3 He. 

The Phase I data taking period started in 2003 and finished in 2004. The detectors 
working at this time was a CCD on the focal plane of an X-ray telescope, a TPC 
and a Micromegas. In these conditions, CAST obtained an experimental limit on 
the coupling constant of g ai < 8.8 x 10 -11 GeV -1 at a 95% of C.L. for axion 
masses m a < 0.02 eV [69]. 

In order to restore the coherence to higher axion masses the magnet bores were filled 
with a buffer gas during Phase II. From 2005 to 2006 4 He was used as buffer gas, 
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Figure 2.11: Snapshot of the slow control program of the CAST experiment 

providing an experimental limit of g a7 < 2.17 X 10- 10 GeV” 1 at 95% of C.L. for 
axion masses 0.02 < m a < 0.39 eV [70]. Since the saturation pressure of 4 He at 
1.8° K is around 17 rnbar and the 4 He could liquefy, CAST did not go to higher 
pressures in order stay within safety limits. 

With the purpose of scanning higher axion masses the 4 He was replaced by 3 He. In 
this case the saturation pressure at 1.8° K is around 135 mbar and the axion mass 
range could be extended up to ~1.2 eV. During 2007 the CAST gas system was 
upgraded and also the TPC detector was replaced by two Micromegas detectors. 

The 3 He phase started in 2008 and ended up in 2011. During this period, axion 
masses from 0.39 < m a < 1.17 eV were scanned and different experimental limits of 
g a7 < 2.3 x 1CP 10 GeV -1 for 0.39 < m a < 0.64 eV [71] and g a7 < 3.3 x 1CP 10 GeV -1 
for 0.64 < m a < 1.17 eV at 95% of C.L. [72] were obtained. Thanks to these 
achievements, the KSVZ line (one of the most favored by theoretical models) was 
crossed for the first time. The excluded regions during the different periods are 
shown in figure 2.12. 

In parallel, CAST has also been looking for more exotic axions, like axions with an 
energy of 14.4 keV from the 57 Fe transitions [73], high energy axions from the 'Li 
and D(p, 7 ) 3 He nuclear decays with a gamma-ray calorimeter [74] and non-hadronic 
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m axion( eV ) 


Figure 2.12: Axion-photon coupling limit as a function of the axion mass at the CAST 
experiment. Three different data taking periods are shown: Phase I (CAST Vacuum) in 
black, the 1 He Phase II in blue and the 3 He of the Phase II in red. 

Solar axions [75] that were introduced in section 1.7.1.2. Moreover, during 2013 a 
SDD 11 was installed for more exotic particles searches [76], like chameleons. 

Although CAST finished its original research program in 2011, the data taking pe¬ 
riod has been extended. During 2012 the 4 He phase was revisited, improving the 
previous limit in a narrow mass range [77]. In 2013 CAST started a new data taking 
campaign revisiting the vacuum phase, motivated by the improvement of the back¬ 
ground levels of the detectors. Also, the rescanned vacuum phase continued during 
2014 when a dedicated X-ray focusing device in the focal plane of a Micromegas de¬ 
tector was installed and a considerable improvement of the sensitivity is expected. 
These features will be described in chapter 7. 
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2.4 X-ray detectors in the CAST experiment 

The X-ray detectors at CAST are installed at the magnet bore ends. During the 
sunrise, solar axions that enter through the MRB could be converted into photons 
inside the magnetic field and detected in the X-ray detectors installed on the other 
side, close to the MFB, that are referred as Sunrise detectors (see figure 2.2). During 
the sunset the opposite process occurs and thus, the detectors close to the MRB are 
called Sunset detectors. 

While the magnet is not pointing the Sun the detectors are taking background data. 
Later on, tracking and background levels are compared and an excess of counts 
during tracking might indicate an axion signal. So the sensitivity of CAST can be 
improved by lowering the background level of the X-ray detectors. 

Mainly three different kind of X-rays detectors have been working at CAST since the 
beginning of the experiment: a Charge Coupled Device (CCD) in the focal plane of 
an X-ray telescope, a Time Projection Chamber (TPC) covering two magnet bores 
and finally several MICRO MESh GAseous Structure (Micromegas) detectors, for 
which different technologies were used in different periods. The principle of operation 
of these detectors will be described below. 




Figure 2.13: Left: Photo of the front view of the ABRIXAS X-ray focusing device. 
The red circle represents the CAST magnet bore. Right: Efficiency of the telescope as a 
function of the energy in which two cases are displayed: the 2003 data taking period (blue 
line) and the following campaigns of data taking (green line). The loss of efficiency is 
due to the realignment of the telescope during 2004 in order to center the spot in the CCD 
detector. 


2.4.1 The X-ray telescope and the CCD system 

The use of an X-ray telescope is one of the innovations of CAST. The X-ray telescope 
is installed in the Sunrise side of the magnet and focuses the total magnet aperture 
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of 14.52 cm 2 to a spot of 9 mm 2 on a pn-CCD detector. The advantage of the use of 
an X-ray focusing device is that the expected signal area is much smaller and thus 
the signal to background ratio is increased, in this case by a factor 150. 



Figure 2.14: Left: Photo of the pn-CCD detector installed at CAST. Right: Projection of 
the focusing spot in the CCD, the image has been obtained during the alignment procedure 
with a laser. The white circle corresponds to the core of the Sun with a size of 0.2 solar 
radius. Plot taken from [80]. 

The X-ray telescope operating at CAST is a Wolter I prototype built for the 
ABRIXAS 12 mission [78]. It is made by 27 nested and gold coated parabolic and 
hyperbolic nickel mirror shells with a diameter from 76 mm to 163 mm and a focal 
length of 1600 mm. The telescope is divided into six sectors and only one is big 
enough to cover the entire magnet aperture of 43 mm of diameter. The X-ray tele¬ 
scope was fully characterized at the PANTER [79] facilities of the MPE 14 and the 
sector with the best performance was selected (see figure 2.13 left). The use of an 
X-ray telescope entails a loss of signal efficiency as it is shown on the right side of 
figure 2.13. 

The CCD detector [80] operating at CAST (see figure 2.14 left) is a prototype 
of the ESA 14 XMM-Newton 15 [81] mission. It is fully depleted and 208 pm thick 
with a sensitive area of 2.88 cm 2 distributed on 200x64 pixels with a size of about 

12 A BRoad Imaging X-ray All-sky Survey 

lj Max Planck Institute for Extraterrestrial Physics, Munich 

14 European Space Agency 

15 X-ray Multi-mirror Mission named Newton 
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150x150 //in 2 each. It is placed on the focal plane of the X-ray telescope where the 
focusing spot has a diameter of about 19 pixels (see figure 2.14 right). The pn-CCD 
detector operates at a nominal temperature of 143° K and under vacuum. 



Figure 2.15: The x-ray telescope and the CCD detector system in the CAST experiment. 


The advantages of this kind of solid state detector are a high quantum efficiency 
(close to one from 1-10 keV), an excellent energy resolution (~ 3% of FWHM at 
6 keV) and a good spatial resolution. It allows to discriminate background events 
from X-rays via pattern recognition. The nominal background level of the CCD 
detector at CAST is about 5 — 8 x 10 -5 c cm _2 s _1 keV _1 from 1-7 keV, which is 
worse than other detectors working at CAST and may be induced partially by the 
internal radioactivity of the materials surrounding the detector [82], 

The X-ray telescope and the CCD are kept under vacuum in order to avoid a con¬ 
tamination of the telescope which could produce a loss of efficiency. The system 
has additional gate valves that separate the different subsystems: magnet, telescope 
and CCD detector (see figure 2.15). The CCD is aligned with the X-ray telescope 
using a parallel laser that can be detected in the CCD. The position of the focusing 
spot is checked regularly with an X-ray finger source installed on the other side of 
the magnet. 
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2.4.2 The Time Projection Chamber 

The TPC was taking data at CAST from 2003 to 2007 covering the two magnet bores 
of the Sunset side. It is a gaseous ionization detector that combines the technology 
of the Multi Wire Proportional Chambers (MWPC) [83] and drift chambers. In the 
CAST TPC two different regions can be distinguished: the conversion region and 
the amplification region(see figure 2.16). In the conversion region the interacting 
particles ionize the gas and the electrons are drifted to the amplification region 
due to the electric held applied (~700 V cm -1 ). The amplification region is made 
between anode and cathode wires in a strong electric held of about 5 kV cm -1 , here 
an avalanche process take place amplifying the signal. The resulting electrons and 
ions generated in the avalanche are captured in the anode and cathode wires, which 
are transversely distributed, thus conferring spatial resolution to the detector. 

Anode plane Cathode planes 
of wires of wires 


for X-Rays rin 9 s 

Figure 2.16: Scheme of the TPC of the CAST experiment in which the different part are 
listed. Plot taken from [84] 

The CAST TPC [84] has a conversion volume of 10 X 15 x 30 cm 3 . The drift region 
is parallel to the magnet bores and has a length of 10 cm, while the amplification 
region has a gap of 3 mm. Its cross section (15 x 30 cm 2 ) allows to cover the 
two magnet bores by the use of 4 /rnr thick aluminized mylar windows glued to a 
metallic grid (also called strongback). The thin mylar windows separate the vacuum 
side from the detector and allow the transmission of the X-rays without a significant 
loss of efficiency. The TPC operates at atmospheric pressure and in a gas mixture 
of 95% Ar and a 5% of CH 4 . 

The two dimensional readout of the CAST TPC detector allows the offline discrimi¬ 
nation of the background events via pattern recognition with a high rejection factor. 
In addition, a shielding was installed during 2004 in order to reduce the background 
level. The shielding was designed according to the results of simulations of the envi¬ 
ronmental 7 ’s which were measured in the experimental area [85]. It was composed 
by different layers of materials with different purposes (see figure 2.17). In the in- 
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TO THE 



Figure 2.17: Left: Scheme of the CAST TPC shielding. Right: Photo of the experimental 
set-up. 



nermost part, the TPC was covered by a copper box with a thickness of 5 mm, 
used as a Faraday cage in order to reduce the electronic noise and also to block 
low energy photons. The second layer was made of 2.5 cm of lead that attenuates 
medium and high energy environmental 7 ’s. The external part of the shielding was 
composed by 22.5 cm of polyethylene with a cadmium sheet in the innermost part. 
Environmental neutrons could be thermalized through the polyethylene and then 
absorbed in the cadmium sheet. The background level was reduced by a factor 4.3 
after the shielding upgrade. 

The background level after the upgrade was about 4 X ICC 5 c cuD 2 s _ 1 keV _1 from 
2-10 keV. This background is worst than the levels that could be reached by the Mi- 
cromegas detectors that will be described in the following chapters. The CAST TPC 
was replaced by two Micromegas detectors during 2007, motivated by the improved 
background and the better performance of the novel Micromegas technology. 

2.4.3 The Micromegas detectors 

Different types of Micromegas detectors have been installed at CAST in different 
periods and places. At the beginning of the experiment only one Micromegas de¬ 
tector with a classical technology was installed in the Sunrise side of the magnet. 
Afterwards, novel technologies in the manufacturing techniques bulk and microbulk 
were developed and installed at CAST. In 2007, the TPC of the Sunset side was 
replaced by two Micromegas detectors, since then, three of the four X-ray detectors 
operating at CAST are of the Micromegas type. 

The Micromegas detectors are a technological evolution of the MWPC and drift 
chambers. Micromegas are in the frame of the novel Micro-Pattern Gaseous Detec¬ 
tors (MPGD) technology. In contrast to the classical TPC’s the wired readout has 
been replaced by a printed circuit board (PCB). The Micromegas detectors were 
developed by Giomataris in 1996 [ 86 ]. The main feature was the introduction of 
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a thin parallel electrode (also called mesh ) over the printed circuit with a narrow 
gap (50-100 /im), that makes the electric field constant in the amplification gap 
obtaining a better homogeneity of the gain. Another advantage is that the drift and 
the amplification regions are decoupled, that allows to optimize the different electric 
fields separately. The main technological challenge was the homogeneity of the gap 
between the thin mesh (3-5 pm) and the PCB. Therefore, different manufacturing 
techniques were developed from the classical Micromegas to the bulk and microbulk 
technologies. 

2.4.3.1 Classical Micromegas 

It was the technology of the first Micromegas installed at CAST on the Sunrise 
side in which the mesh and the readout planes were built separately [87]. The 
electroformed mesh and the anode were separated by periodic insulator spacers, 
also called pillars , made of kapton 16 that was deposited by standard lithographic 
methods at the anode. The mesh was stretched and glued to a frame and placed 
on the pillars. The electrostatic forces between the mesh and the anode ensure the 
homogeneity of the gap. However, this process was complicated and an expertise 
operator was required. Also, the performance of this type of Micromegas was limited 
by the manufacturing technique. 

2.4.3.2 Bulk Micromegas 

In the Micromegas detectors of the bulk type [88], the mesh is made by a commercial 
woven wired material (Au, Cu, Fe, Ni and Ti) of about 30 pm thick and the manu¬ 
facturing process is relatively easy. The woven mesh is stretched and encapsulated in 
a photoresistive material, named vacrel, and glued on the top of the anode readout 
(see figure 2.18 left). Then, the vacrel is etched using a photolithographic method 
creating the pillars (see figure 2.18 right). In order to guarantee the homogeneity of 
the gap the amplification region has a thickness of 128 pm. 

The Micromegas bulk detectors are robust and inexpensive, also it is possible to 
construct large areas. They have an acceptable energy resolution and a high gain. 
The advantages of this manufacturing technique motivate the replacement of the 
classical Micromegas of the Sunrise side at CAST by this novel bulk technology. On 
the other hand, the performance of the bulk detectors is not as good as the microbulk 
Micromegas. Also, the PCB readout is not radiopure which is not desirable for low 
background experiments. 

2.4.3.3 Microbulk Micromegas 

In the microbulk Micromegas [89] the anode and the mesh are produced together 
during the manufacturing process (see figure 2.19 left) by the use of a double sided 
copper coated kapton foil (copper-kapton-copper) of about 50 pm thick. The pattern 

16 Flexible polyamide foil 
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Base Material 


Lamination of Vacrel 


Positioning of Mesh 


Encapsulation of 
Mesh 



Figure 2.18: Left: Scheme of the bulk manufacturing process. Right: Photo of a bulk 
detector showing the pillars (green). Plots taken from [88]. 


anode readout is made by removing the copper with a photolithographic process. 
In addition, several single side coated kapton foils may be attached for the require¬ 
ments of the readout, these foils are etched and different vias are constructed in 
concordance with the readout pattern. Finally, the mesh grid is made by a photo¬ 
chemical process creating a hole pattern (see figure 2.19 right). Then, the kapton is 
etched and removed in order to build small pillars. 
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Figure 2.19: Left: Scheme of the microbulk manufacturing process. Right: Microscopic 
view of the mesh holes pattern in a microbulk detector. Plots taken from [89] . 


The manufacturing process of the microbulk Micromegas detectors shows the highest 
homogeneity in the amplification gap. This leads to the best energy resolution 
reached in a Micromegas detector with a high gain. Also, the materials used in the 
construction, mainly copper and kapton, are radiopure. The better performance 
of this technology motivates the replacement of the bulk Micromegas working at 
CAST. In contrast with bulk detectors, microbulk Micromegas are more expensive 
and large areas can not be easily constructed. 

The Micromegas detectors are very versatile and are used in different experiments 
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such as T2K 17 [90], nTOF 18 [91], COMPASS 19 [92] and they are projected in many 
other experiments like in MAMA 20 [93] for the ATLAS detector upgrade at the LHC 
and MIMAC 21 [94], These experiments include accelerator and rare event searches, 
in which different technologies as bulk and microbulk are used. Also, CAST is the 
best representative experiment for Micromegas detectors in which three of the four 
detectors currently installed are of the microbulk type. A more detailed description 
of the Micromegas detectors working at CAST together with the different detector 
systems will be introduced in next chapter. 


17 Tokai to Kamioka 
18 neutron time-of-flight 

19 Common Muon and Proton Apparatus for Structure and Spectroscopy 

20 Muon Atlas Micromegas Activity 

21 Micro TPC MAtrix of Chambers 
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3.1 Introduction 

The Micromegas detectors at CAST have shown an extraordinary evolution in terms 
of background level and detector performance since the beginning of the experiment. 
These improvements could be explained by the development of new manufacturing 
techniques as it was shown in section 2.4.3.3. For a better understanding of the 
processes involving this kind of detectors, a briefly description of the interaction of 
particles in gaseous detectors will be presented. Also, the working principle of the 
Micromegas detectors will be detailed. 

The understanding of the processes involving the detection of axions at CAST comes 
through a deep knowledge of the detector systems. For this reason, the main features 
of the different lines at the Sunset and Sunrise side will be described. Finally, the 
evolution of the background levels of the Micromegas detectors in different data 
taking campaigns will be presented. 

3.2 The Micromegas detectors 

Micromegas are gaseous ionization detectors with two different regions separated by 
a mesh: the conversion region and the amplification region. In the conversion region 
the interacting particles ionize the gas generating primary charges, these interactions 
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will be described in section 3.2.1. In the conversion region an electric field is applied 
and the electrons are drifted to the amplification region. In this region an avalanche 
of the primary electrons occurs, due to the strong electric field present in the gap, 
generating a readable signal. These processes which involve the signal generation in 
the Micromegas detectors will be detailed in section 3.2.2 

3.2.1 Interactions of particles in gaseous detectors 

In order to detect a particle, it must interact with the detector and transfer energy in 
a readable manner. This energy transfer can be translated into ionization, generating 
electron-ion pairs (in ionizing detectors) and into the excitation of the atoms of the 
detector. Also, the deexcitation of these atoms can be performed mainly via three 
different channels that can be measured: heat (or phonons) in bolometers, light in 
scintillators and electron-ion pairs in ionization detectors. Micromegas are gaseous 
ionization detectors and thus this section will focus on the related processes. 


3.2.1.1 Interaction of charged particles 

The energy transfer of charged particles through gaseous detectors is mainly done 
by electromagnetic interactions in which the most dominant process is the Coulomb 
scattering. This process could be explained by the interaction of a charged particle 
in the Coulomb field of the atoms of the detector. It yields an energy transfer from 
the incoming particle to the detector ionizing the gas. The energy loss of the charged 
particles per unit length is given by the Bethe-Bloch [95, 96] equation: 


(IE Z 47 rz 2 e 4 I" / 2 m e c 2 f3 2 \ 

_ = NAP Aphn e c 2 [ n V/(l-02) ) 



(3.1) 


here Na is the Avogadro’s number, Z and A are the atomic number and weight while 
p is the density of the target material; z is the charge and (3 the velocity relative 
to c of the incoming particle (v = /3c)] m e and e are the electron mass and charge 
respectively and I is the mean ionization potential in the target material. 

Although the Bethe-Bloch equation is valid for moderate relativistic charged heavy 
particles, additional corrections are necessary in equation 3.1 for the density effect 
at high energies and shell corrections at low energies. The shape of the Bethe-Bloch 
equation is shown in figure 3.1. The energy loss decreases as the energy increases 
until a minimum, followed by an almost constant value also called Fermi plateau. 
Thus, for soft relativistic particles (10 < /3q < 100 where 7 = (1 — /l 2 )” 1 / 2 ) the 
energy loss is independent of the energy of the particle and is relatively weak. These 
kind of particles are called minimum ionizing particles (MIPs) and generically the 
energy loss can be approximated to ~2 MeV g -1 cm -1 . 

The main task in the Bethe-Bloch equation is the determination of the mean ex¬ 
citation energies that have to be estimated by experimental stopping powers mea¬ 
surements [97] (see figure 3.2 left). Moreover, the energy loss is the average value 
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Energy (MeV) 

Figure 3.1: Bethe-Bloch parameterization of the stopping power in Ar for different par¬ 
ticles: electrons (red), muons (blue) an photons (black). Equation 3.1 has been used for 
protons and muons and equation 3.2 for electrons. 


of many collisions into the absorber and it is only valid for thick or dense enough 
materials. In the case of thin absorbers, the energy loss can be described by the 
Landau distribution [98], in which the most probable value is lower than the one 
extracted from the Bethe-Block equation. However, for very thin absorbers such as 
TPC’s, the distributions are wider than the original Landau function (see figure 3.2 
right). 

Furthermore, strong electronic collisions may transfer a big kinetic energy to single 
electrons, generating a secondary ionization. These electrons are known as <5-rays 
and its energy could be large enough to be distinguished as tracks. Moreover, 
electron collisions have a different behavior than heavier particles, because in this 
case the collisions are between particles with similar masses and thus the Bethe- 
Bloch equation from 3.1 needs some modifications. 
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(3.2) 


The energy loss of high energy electrons (above 10 MeV) is dominated by radiative 
interactions (see figure 3.3 top) called Bremsstrahlung. In this interaction an incom¬ 
ing electron is deflected in the electric field of a nucleus and a portion of the kinetic 
energy of the electron is converted into a photon. On the other hand, the energy 
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Figure 3.2: Left: Mean excitation energies (I/Z) for different materials, plot taken from 
[97]. Right: Cumulative energy loss (straggling function) for particles with /3 7 = 3.6 
in 1.2 cm of Ar (solid line), where A p is the most probable energy loss, (A) is the mean 
energy loss from the Bethe-Bloch equation and w is the full width half maximum (FWHM) 
of the distribution. The dashed curve represents the original Landau function. Plot taken 
from [98]. 

loss of low energy electrons (such as (5-rays) is dominated by collisions in which the 
electrons are continuously deflected and its trajectories are not straight, in contrast 
with heavy ions. In this case it is more practical the use of the range of electrons, 
defined as the average path length traveled for a charged particle with a given en¬ 
ergy. A parameterization of the range of low energy electrons (up to ~100 keV) can 
be empirically extracted: 


R = 0.71 E 1 ' 72 (3.3) 

here R is the range in g cm -2 and E is the energy of the electron in MeV. Although 
equation 3.3 reproduces the shape of the range of the electrons, more accurate 
parameterizations can be found in the literature, like the CSDA 1 range from the 
ESTAR 2 [99] database from NIST 3 . The differences between both parameterizations 
are shown on the bottom of figure 3.3. 

3.2.1.2 Interaction photons 

In contrast with charged particles, photons can travel some considerable distance 
without interacting. Also, the interactions of the photons lead to a partial or total 
energy transfer to the medium. So in this case it is more appropriate the use 
of the interaction probability of a photon and its cross-section, which is usually 
parameterized by the intensity I{x) of a beam of photons of a given energy after 
crossing a material with a thickness x, by the expression: 

1 Continuous Slowing Down Approximation 

2 Stopping power and range tables for electrons 

3 National Institute of Standards and Technology 
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Figure 3.3: Top: Stopping power of electrons in Ar. The dashed red line represents the 
energy loss from the electronic collisions and the red dotted line is the contribution of 
the radiative losses (Bremsstrahlung). The solid black line is the total stopping power. 
Bottom: Range ( R/p ) of electrons in Ar as a function of the energy. The black line is the 
parameterization from equation 3.3 while the red triangles have been extracted from the 
ESTAR database. 


I(x) = I 0 e~ ppx (3.4) 

where Iq is the intensity of the beam before crossing the material, p is the density of 
the material and p is the mass attenuation coefficient in cm 2 g -1 . In the literature 
there are also different definitions for this interaction with other interpretations 
related to p. such as: the cross section o = pp/N (where N is the density of atoms 
in the material) or the mean free path A = l/(pp). 
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Figure 3.4: Mass attenuation coefficient of photons in Ar as a function of the energy of 
the incident photon. Three different processes are drawn: photoelectric (red dashed line), 
Compton (blue dotted line) and pair production (green dashed and dotted line). The solid 
black line corresponds to the total mass attenuation, which is the addition of the different 
contributions. Data taken from the XCOM [100] database. 

Photons can interact with matter via mainly three processes: photoelectric effect 
which dominates at low energies, Compton scattering that becomes significant at 
mid energies and pair production at high energies (see figure 3.4). The total cross 
section is given by the sum of the different cross sections of these processes fJ-Total = 
h Photo T I^Compt on 4“ UPair 

In the photoelectric process a photon transfers its entire energy to a bound electron 
of an atom in the medium. This results in the emission of a photo-electron with the 
energy of the incoming photon minus the binding energy of the electron in the atom 
E pe = E. y — Eshell- The total cross section is given by the sum of the different shell 
contributions in which the K-shell dominates the cross section of the process. For 
instance, this effect can be observed in figure 3.4, where it causes a sharp increase 
of the mass attenuation around ~3 keV, that corresponds to the K-shell binding 
energy of the Argon. 

After the photoelectric interaction, the atom may return to its ground state via two 
processes: fluorescence or the emission of an Auger electron. In the fluorescence 
process the vacancy of an inner i th shell is filled with an electron from an upper ] th 
shell with lower binding energy and a photon with and energy E t — Ej is emitted. 
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On the other hand, the total energy available could be released by the emission of 
an Auger electron, however this process is more unlikely. 

The fluorescence processes in the detector may cause the occurrence of additional 
peaks in the spectra. If the photoelectric effect occurs inside the conversion region 
and the X-ray from the fluorescence is not absorbed, the detected energy would be 
^escape = A,, — E s heii, this is the so-called escape peak. If the photoelectric effect 
occurs outside the conversion region, only the X-ray from the fluorescence could 
be detected. These processes only represent a few percent of the events and are 
strongly dependent on the geometry of the detector. They can also be distinguished 
as different depositions of energy in a TPC detector. In general the photo-electron 
and the fluorescence photon are absorbed at the same point in the detector and thus 
the total amount of energy from the photon is detected. 

The Compton scattering becomes dominant at higher energies than the photoelectric 
does. In the Compton interaction a photon is scattered in presence of an electron, 
transferring a part of its energy to the electron. In a first approximation the electron 
can be considered free (not bounded to a nucleus) and the interaction could be 
parameterized by: 


K = - n 7 (3.5) 

7 -| , (1—cos 0)Esy v ' 

JL o 

m e c z 

here E[ f is the energy of the scattered photon and E 1 the energy of the original 
photon; m. e c 2 is the mass of the electron and 0 the dispersion angle of the scattered 
photon. This expression yields a maximum energy transfer Ej^max from the photon 
to the electron for 6 = it: 


Et 


max 


2 E 2 

m e c 2 + 2 E 1 


(3.6) 


So that, for a photon is not possible to transfer more energy to an electron in a 
single collision via Compton scattering. This causes a sharp cut-off at this energy 
which is commonly named Compton edge. Low energy transfers to electrons via 
Compton scattering are relatively common and can have an important contribution 
to the background at low energies. 

Finally, in the pair production process a high energy photon can generate an elemen¬ 
tary particle and its antiparticle in a presence of a nucleus, usually electron-positron 
pairs are generated. Because of the required energy of the process, the pair pro¬ 
duction is only possible at high energies, in principle far away from the y’s of the 
natural radioactivity and thus this process will not be treated. 


3.2.2 Working principles of Micromegas detectors 

As it was explained previously, Micromegas detectors could be divided in two dif¬ 
ferent regions, the conversion region and the amplification region, separated by a 
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transparent conducting mesh (see figure 3.5). In the conversion region ion-electrons 
pairs are generated by the interacting particles in the gas and the electrons are 
drifted to the mesh in the presence of an electric field between the cathode and the 
mesh. The mesh is provided by holes that allow the pass of the primary electrons 
from the conversion region to the amplification region (see figure 3.6). 



Figure 3.5: Scheme of the working principle of a Micromegas detector. Two different 
regions are separated by a mesh: the conversion region and the amplification region. 


The amplification region has a gap of about 50 /jrn thick between the mesh and the 
anode readout, here the electric field is big enough to generate the avalanche of the 
primary electrons. This avalanche process generates a multiplication of the primary 
electrons in which the resulting ions are collected into the mesh and the electrons 
into the anode readout. Finally, two different signals proportional to the primary 
ionization are generated: one in the mesh that provides time resolution and another 
in the anode readout which confers spatial resolution. 

Once reviewed the working principle of Micromegas detectors, this section will be 
focused on the different processes involving the signal generation: primary charge 
generation, drift and diffusion and avalanche multiplication. 
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Figure 3.6: Electric field in a Micromegas detector between the conversion region and the 
amplification region. Both regions have a constant electric field with a rapid transition in 
the mesh holes. An appropriate selection of the ratio between both electric fields makes 
that most of the primary electrons pass to the amplification region. 


3.2.2.1 Primary charge generation 

As it was described before, the ionization takes places when a particle transfers 
energy above the ionization potential to the medium. The electron-ion pair genera¬ 
tion is a probabilistic process and has some fluctuations. The total number of pairs 
produced can be parametrized by the expression: 


E 


np = 


W 


(3.7) 


where np is the number of ion-electron pairs generated, E is the deposited energy 
and W is the average energy required to produce an electron-ion pair. The W 
value depends on the gas of the chamber, for instance W = 26.3 eV in Argon. 
The generation of the primary charges seems to follow the Poisson distribution and 
implies a limitation on the energy resolution of the detectors. However, this process 
is not purely Poissonian and the variance has to be corrected by the Fano factor 
[101] with a -2 = npF (instead of a 2 = np) and thus, the best energy resolution that 
can be reached is given by the expression: 


FW 

R = 2.35 ] j — (3.8) 

here R is the FWHM 1 resolution and F the Fano factor. For instance the Fano 
factor in Argon is F = 0.17 and the resolution is limited to a R = 6% at 6 keV in an 

4 Full Width Half Maximum 
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Argon TPC based detector. Additional effects during the avalanche lead to a worst 
energy resolution. 

3.2.2.2 Drift and diffusion 

After the generation of the primary charges in the conversion region the electrons are 
drifted to the mesh and the ions to the cathode. The ions in the conversion region 
are not involved in the signal generation and will not be treated. The electrons 
are considered to move at a constant velocity during the drift, they experiment 
acceleration due to the electric field, but at the same time are slowed down due to 
the collisions with the gas atoms. The velocity of the electrons u during the drift 
can be parameterized by the expression: 

u =- Et (3-9) 

2 m e v ' 

here e is the electron charge, m. e the electron mass, E the value of the electric field 
and r the mean time between two collisions. The collision cross-section, and thus r, 
has a strong dependence with the electric field due to complex quantum-mechanics 
processes, this is the so-called Ramsauer [102] effect. It makes the calculations of 
the electron velocity really complex and usually are performed via Monte Carlo 
simulations in which the Magboltz [103] software is the most popular. 

The drift velocity can be improved by the addition of small quantities of another 
gas, called quencher , to a pure noble gas. This is due to the fact that in noble gases 
the cross-section is small, but the collisions are not inelastic. On the other hand, 
the quenchers have higher cross-sections but the collisions are inelastic. Then the 
optimum case is a gas with small cross-sections but with very inelastic collisions. 
For instance, the detector gas in the Micromegas at CAST is an Ar and a 2% of 
isobutane(iC 4 Hio) mixture and the drift velocities are around 5 cm [i s^ 1 . 

During the drift the electrons are deviated from its desirable trajectory due to 
collisions with the gas atoms. This phenomenon is called diffusion and may affect 
the topological information of the TPC readout. The transversal diffusion can be 
parameterized by a Gaussian distribution in which the standard deviation a is given 
by: 


a = y/2Dt = 



(3.10) 


here D is the diffusion coefficient, which is dependent on the electric field, t is the 
drift time that can be written as a function of the drift velocity u = l/t, where l 
is the drift distance. The diffusion can be improved by the addition of a quencher, 
because of the increment of the drift velocity. Also, the collisions are more inelastic 
and thus, the diffusion coefficient is minimized. 
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3.2.2.3 Avalanche multiplication 

Electrons in strong electric fields acquire enough energy to ionize more atoms and 
thus, generate additional electron-ion pairs, it results in a chain reaction called 
avalanche. The multiplication in the avalanche can be parameterized by the number 
of electrons N generated after a path dx, by the expression: 


dN 

dx 


= Na 


(3.11) 


here a is the Townsend coefficient, the inverse to the mean free path of the electrons 
in the avalanche. Integrating equation 3.11 the gain factor G can be obtained, 
defined as the number of ionizations created by a single electron in an avalanche: 


N(x) = N 0 e ax G = — = e ax (3.12) 

N 0 

However, the gain cannot be increased further, there is a limitation of G ~ 10 8 given 
by the Raether limit, due to the sparks in the detector. In practice the operation 
point is around four orders of magnitude less than the Raether limit, mainly due 
to imperfections in the mesh or in the anode. Also, the quencher has an important 
role in the gain because it absorbs secondary photons created in the avalanche that 
can provoke sparks in the detector. Although the intrinsic gain is reduced, it allows 
to reach higher voltages in the mesh. 

On the other hand, fluctuations in the avalanche deteriorate the energy resolution of 
the detector. This effect is caused due to the big uncertainties in the number of ion- 
electron pairs generated in an avalanche. These fluctuations can be parameterized by 
the Polya distribution in function of the gas gain. The fluctuations in the avalanche 
together with the uncertainty of the primary electrons described in section 3.2.2.1, 
imply a technical limitation in the resolution. For instance the best resolution 
reached by a Micromegas detector is R ~ 11% FWHM at 6 keV, in an Argon 
isobutane mixture. 


3.3 Micromegas detectors at CAST 

As it was presented in section 2.4.3.3, different types of Micromegas detectors have 
been working in the CAST experiment. This section will focus on one of the latest 
microbulk Micromegas designs, specifically constructed for CAST, as well as the 
detector systems installed during 2011. 

The microbulk detectors at CAST have experimented different design upgrades, 
motivated by the improvement on detector performance and the intrinsic radiopurity 
of the materials in the manufacturing process. In this design, the detector anode 
is made of interconnected square pads. The pads in one direction are connected 
at the level of the anode and the pads of other direction are connected in another 
layer of copper coated kapton foil (see figure 3.7 left). This leads to a 2-dimensional 
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readout of 106 x 106 strips with a pitch of 550 /jru and an effective area of about 
60 X 60 mm 2 (see figure 3.7 middle). The mesh holes are constructed on the top of 
the pads in order to ensure the transmission of the primary electrons. The cathode 
is made of aluminized mylar of about 5 /jrn thick, glued to a circular aluminum 
strongback 5 with a diameter of 110 mm (see figure 3.7 right), which is attached to 
the vacuum side of magnet. The body of the detector is made of Plexiglas in which 
the chamber and thus the conversion volume has 30 mm height. The Micromegas 
is glued to a Plexiglas base called rciqiLette, which is bolted to the chamber of the 
detector. The raquette has a circular shape in the detector area and a thin neck for 
the strips connections to the electronics. 



Figure 3.7: Left: Microscope view of the square pads of the Micromegas readout anode. 
Middle: Detector active area glued on the Plexiglas raquette. Right: Photo of the chamber 
mounted on the Micromegas detector in which the strongback is visible in the top of the 
active area. 

The chamber gas is an Ar + iC. 4 Hio mixture, the fraction of quencher is a 2% in the 
Sunset detectors and a 2.3% in Sunrise. The detectors are working at a pressure of 
1.4 bar, this value has been selected in order to maximize the quantum efficiency 
in the 1-10 keV range. Also, the gas in the chamber is continuously renewed with a 
flow of about 3 1/h. 

The detectors are placed at the magnet bore ends covering the total magnet aperture 
of 14.5 cm 2 . The cathode separates the vacuum side from the Micromegas chamber, 
however small amounts of the gas in the chamber can permeate through the mylar 
window. In order to avoid the deposition of the gas in the cold bores and either in 
the cold windows, a differential pumping system was installed in the detector lines. 
It consist in two regions separated by a 4 ^m mylar window that can be pumped 
separately. The side closest to the detector is referred as bad vacuum while the the 
magnet side is called good vacuum. 

Although the features related so far are common in all the Micromegas systems at 
CAST, there are significant differences between the detectors lines of the Sunset and 

5 Metallic grid 
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the Sunrise side and both systems will be described separately. 

3.3.1 Sunrise Micromegas system 

This design corresponds to the upgrade done during 2007 when it was planned the 
installation of a focusing device in the Sunrise Micromegas line [104], Unfortunately 
the optics was discarded, but the upgrades in the line remained, for this reason the 
detector was placed on the extreme of the magnet platform. 



Figure 3.8: Schematic drawing of the Sunrise Micromegas system in which the different 
parts are labeled. 

A sketch of the Sunrise Micromegas line is shown in figure 3.8. The main novelties 
were the installation of a shielding, upgrades on the gas system, a new vacuum 
system and an automatic calibrator. The shielding was composed of different layers 
of various materials, inspired by the TPC shielding that has been described in 
section 2.4.2. However, in this case the inner shielding has a cylindrical shape 
(see figure 3.9 left) and it is more compact. It consist of 5 mm of copper in the 
innermost part, 25 mm of archaeological lead 6 and 2 mm of cadmium foil at the 
end. The shielding is covered by a Plexiglas cylinder and also, the outermost part 
allows the addition of polyethylene blocks up to 250 mm. 

The gas system is designed to maintain a constant pressure in the Micromegas 
detector. It consist of a manometer in the input and a mass flow controller in the 
output of the gas line (see figure 3.9 middle). The manometer set the pressure in the 
chamber and the mass flow controller regulates the gas flow, with this system the 
pressure remains constant at a working point of 1.4 bar. The pressure and the flow 

°Lead with a negligible concentration of 210 Pb, one of the natural radiative isotopes of lead 
with a half life of T t / 2 = 22.3 y. In order to ensure the decay of the 21t) Pb it had to be melted in 
the antiquity. The Roman lead is one of the most popular. 
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in the system are monitored in the control box and in the slow control. Additionally, 
two electrovalves are installed for safety reasons. 



Figure 3.9: Left: Photo of the Sunrise Micromegas shielding. Middle: Photo of the gas 
panel in which the different parts are labeled. Right: A sketch of the Sunrise manipulator 
installed on the line during a calibration. For background runs the source is retracted and 
hidden inside the flange. 

During normal operation, CAST Micromegas detectors are calibrated at least once 
per day with a ' 55 Fe source (main peak at ~6 keV). Calibrations are essential in 
order to apply a discrimination criterion to the background events (see chapter 4 
for further details). The calibrations are performed by a pneumatic manipulator 
attached to the vacuum line (see figure 3.9 right). During background runs the 
source is hidden inside the flange, while during calibrations the manipulator moves 
the source to the middle of the pipe and thus the X-rays can be detected. The 
manipulator is controlled remotely by the DAQ 1 . 

The vacuum in the line is performed by a differential pumping system. The two 
different regions are connected by a bypass (see figure 3.10). When the pumping 
starts, the bypass has to be opened in order to ensure the integrity of the differential 
window. Only is able to close the bypass when both regions reach a good vacuum 
level and then the different systems are pumped separately. Just in this case is 
possible to open the gate valve VT3 directly connected to the magnet vacuum. 
Moreover, the pressures in the different volumes are monitored and a protection 
system is implemented. In case of failure, the bypass is automatically opened to 
prevent the break of the differential window, also the gate valve is closed to avoid 
the contamination of the magnet bores. 

Furthermore, an X-ray focusing device was installed in the Sunrise Micromegas line 
during 2014. The line has a new detector design and a shielding upgrade, these 
features will be described in chapter 7. 

1 Data AcQuisition system 
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Figure 3.10: Scheme of the vacuum system of the Sunrise Micromegas line. 


3.3.2 Sunset Micromegas system 

The related system was installed during the 2007 upgrade, when the TPC was 
replaced by two Micromegas detectors. The main features are the same as described 
for the Sunrise Micromegas line. However, in this case both detectors share the same 
shielding, vacuum and gas system. These peculiarities will be described below. 



Figure 3.11: Left: Drawing of the Sunset system. Right: Lead shielding covering the two 
detectors. 

The shielding was inherited from the former TPC. A drawing of the Sunset system 
and a photo of the shielding are shown in figure 3.11. In the innermost part, a 
copper box with a thickness of 5 mm covers the two detectors. A second layer made 
of archaeological lead covers the copper box, with a thickness of 50 mm on the top, 
bottom and on the back of the detectors, in the laterals it is only 25 mm thick and 
the side of the magnet is not shielded. The lead is surrounded by a cadmium sheet 
2 mm thick and polyethylene blocks in the outermost part. A new shielding design 
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was installed during 2012, improving the background level in a factor ~4.5. The 
advantages of the new shielding will be described in chapter 7. 
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Figure 3.12: Scheme of the vacuum system and gas system for Sunset Micromegas detec¬ 
tors. 

The gas system is composed by a manometer in the input and a flowmeter in the 
output. Also, a gas reservoir with a volume of 5 1 is installed in the line in order 
to avoid big fluctuations in the pressure. Both detectors are connected in series, 
working with the same gas. In order to perform the vacuum a differential pumping 
system is implemented, it has the same characteristics explained before for the 
Sunrise line. However, in the Sunset side the two pipes are connected and pumped 
at the same time. A scheme of the gas and vacuum system of the Sunset Micromegas 
line is shown in figure 3.12. 

The calibrations are performed with two pneumatic manipulators of the same type 
as described for the Sunrise system. However, in the Sunset line the sources are 
closer to the detectors (see figure 3.13). 

3.4 Micromegas background history at CAST. 

The Micromegas detectors at CAST have experimented a background reduction 
of about two orders of magnitude since the beginning of the experiment. These 
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Figure 3.13: Left: Photo of the Sunset Micromegas in which the calibrators are visible. 
Right: Photo of the sources inside the pipes when the manipulators are in the calibration 
position. 


improvements are mainly due to the better performance of the novel Micromegas 
technology and its radiopurity, a better discrimination of X-ray like events in the 
analysis and different shielding upgrades. 



Date 


Figure 3.14: Micromegas background evolution in the CAST experiment. 
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The evolution of the background level of the Micromegas detectors at CAST is shown 
in figure 3.14. The blue dots correspond to the first Micromegas detectors installed 
at CAST of the classical technology and the background reduction is due to the im¬ 
provement of the discrimination algorithms. During 2007 the Micromegas detectors 
were shielded, also novel technologies as bulk and microbulk were introduced. These 
improvements result in a background reduction of a factor 4.5 represented by the 
green squares in figure 3.14. Finally, the red triangles correspond to the shielding 
upgrade of the Sunset detectors during 2012, also new electronics for the readout 
were installed during 2013. These features together with the current understanding 
of the background will be described in chapter 7. 
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4.1 Introduction 

Once described the physical processes involving the signal generation in the Mi¬ 
cromegas detectors, this chapter will focus on the signal processing. For this pur¬ 
pose a detailed description of the electronics implemented for the mesh and strips 
readout will be presented. Also, the processing of the data from the DAQ to the 
definition of the different observables will be described. Finally, this chapter will 
focus on the discrimination method developed in order to select X-ray like events in 
the background runs. 

4.2 Micromegas readout and electronics 

As it was described in section 3.2.2, there are two different signals in the Micromegas 
detectors: the mesh signal, generated by the ions in the avalanche process and the 
strips signals, induced by the electrons in the anode readout. 

The mesh signal gives information of the deposited energy in the detector and also its 
shape is related with the distribution of the primary charges along the perpendicular 
direction to the readout. Moreover, it is used to generate the main trigger of the 
acquisition. The mesh signal is preamplihed in a first stage using a Canberra 2004 
preamplifier. Afterwards, the signal is shaped and amplified in an ORTEC 471 NIM 1 
timing filtering amplifier, which allows the tuning of different timing parameters. 

1 Nuclear Instrumentation Module 
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The amplified signal is duplicated in a linear Fan In-Fan Out module, one signal is 
sent to a quad discriminator in order to generate the main trigger and other signal 
is sent to a digitizer. The digitizing of the mesh signal is performed by a VME 2 
Matacq 3 [105] board, in which the signal is recorded in a 2500 ns window with a 
sampling rate of 1 GHz and a dynamic range of 12 bits. 

The strips signals give spatial resolution and energy information of the event. The 
strips are read by the front end Gassiplex cards [106], that allows the processing of 
96 channels each and thus four Gassiplex cards are used to read the 106 x 106 strips 
readout. The Gassiplex cards integrate the charge of the strips that are stored in 
different analog memories. In order to perform the acquisition of the strips the 
Gassiplex are controlled by three digital inputs: track & hold, clock and clear signals 
(see figure 4.1). When the DAQ is triggered a delayed track & hold signal is sent 
to the Gassiplex and the strips charges are integrated. Later on, when the first 
clock signal becomes low the analog memories are frozen. Conforming the 96 clocks 
signal arrive (one per channel) the analog values of the memories are sent through 
the output of the Gassiplex and multiplexed. When the clocks ends, the track & hold 
signal stops and a clear signal is sent in order to reset the analog memories. The 
input signals are controlled by a VME sequencer (CAEN V551B) and the output 
signals are digitized and stored in the CRAMS 4 (CAEN V550) modules. 
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Figure 4.1: Left: Scheme of the Gassiplex input and output signals. Right: Photo of a 
Gassiplex card. 

The different delays in the input signals (T1 and T2 in figure 4.1 left) can be 
optimized in order to get the maximum peaking in the output signal. In this set-up 
one finds experimentally that the output is maximized when the trigger and the 
track Sz hold are delayed ~900 ns and the delay between the track & hold and the 
clock is about ~100 ns. 

2 Versa Module Europa bus 

3 MATrix for ACQuisition 

4 CAEN Readout for Analog Multiplexed Signals 
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During the acquisition of an event there is a period in which the DAQ cannot record 
more events because it is busy, this period is called dead time. A busy signal is 
implemented in the DAQ and the incoming triggers are rejected in this period. 
The dead time is given mainly by the CRAMS and is about 10 ms, limiting the 
acquisition rate up to 100 Hz. 

Although the features described before are common for Sunset and Sunrise DAQs, 
there are significant differences between them and have to be described separately. 

4.2.1 Sunrise acquisition and software 

In this case the interface between the computer and the VME is performed by a 
NI 5 VME-MXI2 controller with a GPIB 6 connection. The DAQ software controls 
and communicates with the different VME modules in order to perform the acqui¬ 
sition. Apart of the Matacq, sequencer and the CRAMS described above, there are 
additional modules like a CAEN V560 scaler to count the total number of triggers 
and a CAEN V262 input-output register in order to control the manipulator for the 
calibrators. Also, several NIM modules are used to process the different signals. A 
scheme of the electronic chain is shown in figure 4.2. 



COMPUTER 


Figure 4.2: Sunrise Micromegas electronic diagram. 
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A Labview based software with a GUI' was developed in order to perform the 
acquisition, also different tools for on-line visualization were implemented (see figure 
4.3). There are predefined different types of runs such as pedestals, calibration and 
background. During background and calibration runs the DAQ performs a normal 
acquisition, the only difference is that the source is moved to the proper position. 
For the pedestal runs an external trigger is generated and only the strips values are 
stored. These values are used to estimate the noise level in the readout, defining a 
pedestal value for every strip. This method will be explained in section 4.3. 




Figure 4.3: Snapshot of the Labview based GUI in the Sunrise DAQ. Different windows 
are shown: RunControlAll.vi (top left) to start and stop the acquisition, RunControlMoni- 
tor.vi (top right) to set the different parameters in the DAQ and RunEv.vi (bottom) for the 
on-line visualization of the events. 

The data is stored in a binary file of 4 bytes array words. It starts with a file 
header that contains information about the run (start time, configuration, type, etc). 
Moreover, every single event is provided with an event header and an event footer, 

‘ Graphical User Interface 
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between them, mesh and strips info is written. The events are stored sequentially 
until the footer of the file which marks the end. 

4.2.2 Sunset acquisition and software 

The main difference with Sunrise is that in the Sunset DAQ both detectors are 
acquired together at the same time. The different mesh signals are amplified sep¬ 
arately and the resulting logic signals after the discriminators are passed through 
a logic OR gate that gives the general trigger to the DAQ. Also in this case four 
CRAMS are used to acquire the strips from both detectors. Another difference is 
that a CAEN V718 bridge with and optical link is used to communicate with the 
computer, instead of the GPIB connection used in Sunrise. A detailed scheme of 
the electronic chain is shown in figure 4.4. 

The acquisition software is a modified version of the Sunrise DAQ, which was ex¬ 
tended in order to acquire two detectors (see figure 4.5). In the Sunset DAQ two 
different calibration types, one per detector, are defined. Also, an autopilot mode is 
implemented with which the DAQ acquires in a predefined schedule and the differ¬ 
ent run types are acquired automatically. When the detector is in autopilot mode, 
calibrations, pedestals and background runs of both detectors are written in the 
same binary file. This creates a more complex file, however the different run types 
and detectors are labeled and in the first stage of the data processing the file is split 
and the related data are separated. 


4.3 Data analysis 

In a first stage of the analysis the information from the DAQ is processed and 
different observables are defined, the analysis software is written using standard 
C++ and ROOT [109] libraries. Also, the gain is calculated with the information of 
the observables related with the deposited energy in the calibrations files. This gain 
factor is applied later on to the corresponding background runs. Finally, a ROOT 
hie is generated, in which the different observables that define different parameters 
of the event are stored. 

For the mesh signal a pulse shape analysis (PSA) is implemented, in which the pulse 
is defined by different parameters. For the strips, a cluster analysis is performed 
after the pedestal subtraction, defining additional observables related with the strip 
readout. Both analysis methods will be described below. 

4.3.1 Pulse shape analysis 

The mesh pulse is parameterized by different observables that describe its shape (see 
figure 4.6). The energy of the event is extracted from the pulse height (amplitude) 
and its integral. The topological information of the pulse is mainly parameterized 
by the risetime and the width. 
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Figure 4.4: Electronic diagram of the Sunset Micromegas acquisition. Plot taken 
from [107], 
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Figure 4.5: Snapshot of the Labview based GUI in the Sunset DAQ. The different windows 
of the DAQ are shown: RunControlAll.vi (top left) to start and stop the acquisition and set 
the autopilot, RunControlMonitor.vi (top right) to set the different parameters in the DAQ 
and RunEv.vi (bottom) for on-line visualization. 


The pulse parameters can have different values depending on the definition used. 
The characteristics of the analysis method implemented are summarized below: 

• Baseline: Voltage offset of the pulse, calculated as the mean value of the first 
500 points. 

• Amplitude: The maximum value of the pulse after the baseline subtraction. 

• Integral: The integral of the pulse since its rise from 15% of the amplitude, 
until it drops to the 15%. 

• Risetime: Time difference between the rise of the pulse at a 15% of the am¬ 
plitude, until it reaches the 85%. 

• Width: Time difference between two points at a 50% of the pulse amplitude. 
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Figure 4.6: Sample of a mesh pulse recorded in the DAQ. The definitions of different 
observables are shown. 

Before the definition of the different pulse parameters, the pulse is processed by a 
FFT 8 analysis, in which high frequencies are subtracted. It reduces the fluctuations 
induced by noise and its shape is smoothed. 

4.3.2 Cluster analysis 

The analysis of the strips readout is performed by a cluster analysis in which the 
different charge depositions are parameterized. A cluster is defined as a consecutive 
number of triggered strips. In order to distinguish if a strip has been triggered or 
not, a pedestal value is defined. 

Every strip has a characteristic noise level that could be parameterized by its 
pedestal value. The pedestal calculation is made by computing the mean rrij and 
the standard deviation a t of the noise level for large number of events in every single 
strip. A strip is considered triggered when its value s* is above the mean plus three 
times its standard deviation: 


Si > rm + 3cjj (4.1) 

For the pedestal calculation a set of 15000 events is used, also it can be extracted 


Fast Fourier Transform 














4.3. Data analysis 


79 


from pedestal, background and calibration runs. For the two last cases a method to 
reject triggered strips was implemented. 




Figure 4.7: Single cluster event in the Micromegas readout after the pedestal subtraction. 
Two different plots are shown: one for the strips of the x-axis (left) and another for the 
y-axis (right). The red bars represent the triggered strips, while the blue bars are the non 
triggered strips. 

After the pedestal subtraction a cluster analysis is performed. In a first stage the 
two different readout planes are scanned using a cluster finder algorithm. A cluster 
is defined as two or more strips triggered consecutively. However, in some detectors 
there are problematic strips that have to be removed or are not visible. To avoid 
this problem, clusters up to two consecutive non triggered strips are allowed. The 
third consecutive non triggered strip marks the end of the cluster. After the one 
dimensional cluster finder different parameters are defined: 

• Cluster charge : Addition of the measured charge of the strips inside the 
cluster. 


• Cluster position: Mean position of the cluster weighted by the measured charge 
in each strip. 



I 


• Cluster size/sigma: Standard deviation of the charge in a cluster. It gives 
information about the shape of the cluster. 

a 2 = 1 c) 

c 

i 


Cluster multiplicity: Number of triggered strips inside a cluster. 
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An X-ray event in the detector generates a point-like deposition of the charge, 
producing single clusters events in the readout (see figure 4.7). So in principle 
events with more than one cluster might be rejected. However, noise in the strips 
or crosstalk could trigger strips, producing non physical clusters. For this reason a 
two dimensional cluster analysis is performed, by introducing a new parameter, the 
cluster balance. 

For events that have more than one cluster, all the possible cluster combinations 
are computed and the cluster with the higher charge is defined as main cluster. If 
the main cluster contains more than the 70% of all the deposited charge, the event 
is stored only taking account the charge inside the main cluster , if not the event 
is rejected. It allows the definition of the cluster balance as the deposited charge 
in the main cluster divided by all the deposited charge of the event. Usually more 
than the 95% of the calibration events are accepted in the analysis, being the cluster 
balance above a 90%. However, these values can be worse in detectors with a bad 
performance or with a low intrinsic gain. 

4.4 Background discrimination method 

As it was presented in section 1.7.1 the axion signal would be an excess of X- 
ray events while the magnet is pointing the Sun. For this reason a discrimination 
algorithm has been developed in order to select X-ray like events, in the CAST Range 
of Interest (Rol) from [2 — 7] keV. Therefore, the most representative observables 
from the daily 55 Fe calibration are selected. 

The 55 Fe decays to 55 Mn by electron capture, emitting characteristics X-rays of 
5.9 keV with a probability of a 25.4% and 6.5 keV with a probability of a 3%. 
These two different peaks can not be distinguished in the Micromegas because of 
the resolution of the detector. Nevertheless, the 5.9 keV peak shows a bump in 
the right side, which can be fitted. Also, the escape peak of the Ar explained in 
section 3.2.1.2, with an energy of about ~3 keV, is observed during calibrations. A 
typical spectrum of a Micromegas detector with a 55 Fe source is shown in figure 4.8. 
The 5.9 keV and the 3 keV peaks are used in order to define a selection criteria for 
the background runs. 

The most significant observables described in previous sections or combinations 
of them are used to define a selection criterion of the X-ray events. This set of 
observables is summarized in the following list: 

• Pulse risetime: Point-like charge depositions as X-ray events give small rise¬ 
time values. 

• Pulse width: It is a measure of the topology of the event in the z-direction. 
X-ray events usually have small values. 

• Energy ratio mesh: Ratio between the mesh pulse amplitude and the mesh 
integral with an expected value of 1. 
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Energy (keV) 


Figure 4.8: 55 Fe spectrum in a Micromegas detector in which the three different peaks 
have been fitted: the escape peak (green line), the 5.9 keV peak (red line) and the 6.5 keV 
line (blue line). The solid black line is the total fit, computed by the addition of the three 
Gaussians. 


• Cluster size/sigma X: Standard deviation of the charge in a cluster in the 
X-axis. X-ray events show a narrow distribution. 

• Cluster size/sigma X : Standard deviation of the charge in a cluster in the 
Y-axis. 

• Cluster multiplicity X: Number of triggered strips in the X-axis. 

• Cluster multiplicity Y : Number of triggered strips in the Y-axis. 

• Cluster balance: Charge of the main cluster in the strips readout divided by 
the overall deposited charge of the event. 

• Energy ratio mesh/strips : Ratio between the mesh pulse amplitude and the 
strips charge with an expected value of 1. 


For the observables listed before, the probability distribution is computed for the 
5.9 keV calibration events. These observables are allowed in a wide range and divided 
in bins, finally the distribution is normalized dividing by the number of events (see 
figure 4.9). The ratio of probabilities that an event has a certain observable value x 
is given by the odds: 


O(x) 


Pi?) 

1 - P(x) 


(4.2) 
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Figure 4.9: Normalized probability for several observables used in the analysis. The 
different observables have been labeled. 
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where P(x) is the probability of a given observable to have a value x and 0(x ) are 
the odds. The ratio of probabilities of an event to have different observables values 
Xi is given by the multiplication of the odds: 

n = J] O t (x t ) (4.3) 

i 

here II is the multiplication of the odds of the different observables denoted by i. 
Computing the logarithm of the odds, log-odds [110, 111], a value is obtained for 
every single event that gives information about how likely [112] is the event to an 
X-ray. 


- log{n) = - ^log(Oi(x;)) = - ^log(T^(xj)) + ^log(l - P*(xj)) (4.4) 

i i i 



Figure 4.10: Distribution of the log-odds for the 5.9 keV peak (blue), 3 keV peak (green) 
and background events (red). The blue filled area corresponds to a software efficiency of 
a 75% (cut-value = 12.8) at 5.9 keV and the green area to a 50% (cut-value = 13.9) of 
efficiency at 3 keV. 

Finally, the log-odds value is computed for calibration and background events. For 
the calibration events two different cumulative distributions are obtained, one for 
the 5.9 keV peak and another for the 3 keV escape peak. It allows to define a 
cut-value in the log-odds distribution below which a certain number of events are 
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accepted. For background events with energies above 4 keV the cut-value from the 
5.9 keV distribution is used, while for events with energies below 4 keV the cut- 
value from the 3 keV distribution is applied. This discrimination method is shown 
in figure 4.10, in which the distribution of the X-ray events during the calibrations 
is peaked to the leftmost part of the figure. On the other hand, the background 
events are uniformly distributed with a bump on the right part. This results in a 
rejection factor of about the 99% of the background events. 

The cut-values are determined by requiring a given software efficiency, defined as 
the number of accepted events divided by the total number of events at 3 and 5.9 keV 
during the calibrations. The software efficiency and its implications in the analysis 
will be revisited in section 5.3.2. 
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5.1 Introduction 

This chapter will focus on the results of the analysis of the Micromegas detectors 
during the 2011 data taking campaign at CAST, using the background selection 
method presented in section 4.4. Moreover, the capabilities of the background re¬ 
jection have dependence with the detector performance and will be studied. Also, 
the analysis leads to the definition of a software efficiency that has to be optimized 
in order to maximize the sensitivity. 

Finally, the determination of the tracking times and a study of the compatibility of 
the tracking and background events will be presented. 

5.2 CAST Micromegas detectors performance 

In this section the performance of the microbulk Micromegas detectors taking data 
during 2011 at CAST will be presented. Three different detectors were working at 
this time: the Mil in the Sunrise side, the M14 in the Sunsetl side and the M9 in 
the Sunset2 side. The gain evolution, energy resolution and the spatial resolution 
of these detectors will be presented. 
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5.2.1 Gain stability 

The gain of the background runs is defined using the closest calibration and also 
the background discrimination is performed with the related calibration run. The 
stability of the gain is a good indicator of the detector performance, so it is desirable 
a stable gain during the data taking period. The evolution of the gain during the 
2011 data taking campaign is shown in figure 5.1. 

The gain of the three detectors shows a good stability along time. However, there 
are some fluctuations whose effect is related with understood events: the reduction 
of the gain in Sunsetl on the 29 th of May is related with the decrease of the mesh 
voltage, while the reduction of the strips gain in Sunrise from the 6 th to the 8 th of 
June is due to an intervention in which the Gassiplex digital inputs were modified, 
also the period without data in Sunrise from the 10 th to the 15 th of June corresponds 
to a power cut in the line. 

5.2.2 Energy resolution 

The energy resolution of the detectors is a good indicator of the detector perfor¬ 
mance. Indeed, the detector on the Sunrise side shows the best energy resolution, 
while Sunset detectors have a modest resolution. The normalized spectra for a set 
of calibrations in the different detectors are shown in figure 5.2. 

The energy resolution is also measured for every single calibration run. A Gaussian 
fit is performed to the 5.9 keV peak and the FWHM (Full Width Half Maximum) is 
extracted, given by the difference between the two values at which the distribution 
is equal to the half of it maximum value. It is related with the Gaussian distribution 
by the expression: 


FWHM(%) = 100 2v/21 ° g2 - ~ 235- (5.1) 

/i fj, 

here o is the standard deviation and [i the mean value of the Gaussian. The evolution 
of the FWHM during the 2011 data taking campaign is shown in figure 5.3. 

The energy resolution shows a good stability along time for all the detectors. How¬ 
ever the Sunset2 detector shows a improvement of the resolution on the 1 st of May, 
related with an intervention in which the noise was reduced. The Sunrise detector 
has a very good energy resolution of about a 18%, being one of the best detec¬ 
tors working at CAST. On the other hand, Sunset detectors show a poor energy 
resolution, which also could be translated in a worst discrimination capabilities of 
the background events. Due to its modest performance, the Sunset detectors were 
replaced during 2012. 

5.2.3 Spatial resolution 

The spatial resolution of the strips readout is also a good indicator of the detector 
performance. Some detectors show dead areas which correspond to defects in the 
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Figure 5.1: Gain as a function of time during the 2011 data taking campaign. The three 
Micromegas detectors are shown: Sunsetl (top), Sunset2 (middle) and Sunrise (bottom). 
The blue circles correspond to the ADC value of 5.9 keV peak in the mesh during the 
calibrations and the red squares show the gain of the strips after the cluster analysis. 


strip readout. These defects are usually related with disconnected strips because of 
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Figure 5.2: Normalized spectra for a set of calibrations. Different detectors are shown: 
Sunset 1 on the top, Sunset2 on the middle and Sunrise in the bottom. The plots on the left 
show the spectra of the mesh (black bars) while the plots on the right side correspond to 
the spectra in the strips readout. 


a problem during the manufacturing process or strips in shortcut with the mesh that 
have to be removed. The spatial resolution of the different Micromegas detectors is 
shown in figure 5.4, in which the cumulative distribution of the mean cluster position 
for a set of calibration, also called hit-map, has been computed. 

The hit-map has a characteristic pattern that corresponds to the metallic grid cath¬ 
ode of the detector called strongback, described in section 3.3. In the Sunrise de¬ 
tector the strongback is clearly visible, due to its high spatial resolution, however 
some strips in the y-axis are missing inside the coldbore area. Sunset detectors have 
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Figure 5.3: FWHM as a function of the time during 2011 data taking period. The three 
Micromegas detectors working at this time are shown: Sunsetl (top left), Sunset2 (top 
right) and Sunrise (bottom). The blue circles correspond to the mesh and the red squares 
to the strips FWHM at the 5.9 keV peak. 


a poor spatial resolution and the strongback is barely visible, also some strips are 
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Figure 5.4: Hit-map distribution for a set of calibrations in the Micromegas detectors at 
CAST. The three different detectors are shown: Sunset 1 (top left), Sunset2 (top right) and 
Sunrise (bottom). The black ring corresponds to the colbore area, where the axion signal is 
expected. Note that the z-axis is in logarithmic scale. 


missing, but outside of the coldbore area. Moreover, this effect can be caused be¬ 
cause the calibration sources are relatively closer to Sunset than in Sunrise (20 cm 
and 1.5 m respectively). 


5.3 Micromegas efficiency 

The detection efficiency of X-rays events from the magnet bores has to be taken into 
account in order to discriminate an axion signal, these features will be described 
in chapter 6. Two different efficiencies can be distinguished in the Micromegas 
detectors: the quantum efficiency in which the loss of efficiency inside the detector 
chamber and through the detector line is computed and the software efficiency, 
introduced in section 4.4, that takes into account the efficiency of the discrimination 
method at different energies. 
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While the quantum efficiency is fixed, the software efficiency can be optimized in 
order to maximize the sensitivity in axion searches. The quantum and the software 
efficiency and its implications in the axion searches will be described below. 

5.3.1 Quantum efficiency 

The Micromegas detector efficiency was measured in the PANTER facilities of the 
MPE 1 . These results showed a good agreement with the simulations [113, 114] per¬ 
formed with the Geantl [115] toolkit. In the simulations the different elements 
present in the line, described in previous chapters, have been implemented: the 
aluminum cathode strongback with its characteristic square grid and the 4 /mi of 
aluminized polypropylene window with a coating of 50 nm of deposited aluminum; 
the detector chamber with a drift distance of 3 cm filled with an Ar+iC^io mixture 
working at a pressure of 1.4 bar; the differential window made of polypropylene with 
a thickness of 4 /mi and the cold window strongback with a 14 /mi polypropylene 
layer. Finally, the detector response was implemented by introducing the energy 
resolution. In order to do it, the efficiency has been convoluted with a Gaussian 
function [116]. The results of the simulation with the efficiency loss in the different 
elements are shown in figure 5.5. 




Figure 5.5: Quantum efficiency of the of the Micromegas detectors in the CAST experi¬ 
ment. Left: Results from the Geant4 simulation with the different efficiency loss: chamber 
efficiency (black line), strongback cathode (green line) and differential and cold windows 
(blue line). Right: Quantum efficiency convoluted with a Gaussian function. 


5.3.2 Software efficiency 

The background selection leads to the definition of a software efficiency, which take 
account of the calibration X-ray events that are not accepted by the analysis. Also, 
this efficiency loss has an impact on the sensitivity of the experiment. There is a 
compromise between the background level and the software efficiency that can be 

1 Max Planck Institute for Extraterrestial Physics, Munich 
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quantified in terms of sensitivity by the definition of the figure of merit [117] of the 
detectors F.O.M.det , given by the expression: 

F.O.M. det = £ = £ s£q (5-2) 


here £ is the efficiency of the detector, given by the multiplication of the software 
efficiency e s and the quantum efficiency e q , while b is the background level of the 
detector. Indeed, the sensitivity is maximized by increasing the software efficiency 
with a low background level in the detector. For this reason a systematic study 
for different software efficiencies has been performed. However, the F.O.M.det is 
strongly dependent on the detector performance and each detector has to be com¬ 
puted separately. 

Moreover, the software efficiency at different energies has to be taken into account. 
From the daily calibrations just the efficiencies at 5.9 and 3 keV can be extracted. In 
order to estimate the software efficiency at different energies, the detector response 
has been simulated for different spectra lines in the [1 — 10] keV regime [114]. For 
a given value of the efficiency at 5.9 and 3 keV, the software efficiency at different 
energies are extracted. Theses results have been confirmed by the measurements 
performed in a X-ray line beam in the CAST detector lab [118]. The different 
points obtained from the simulated and measured data are presented in figure 5.6. 



Figure 5.6: Simulated software efficiency (blue dots) in contrast with the measured data 
(red squares). In this case a software efficiency of a 80% at 5.9 keV and a 50% in 3 keV 
peak have been selected. 


The dependency of the software efficiency with the energy leads to the definition of 
a new figure of merit of the detector. And thus, equation 5.2 can be written as: 
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F.O.M. de t = 


n 


y 


(5.3) 


here £j is the efficiency for the energy bin i, bi is the background level and <hj is 
the normalized solar axion spectra, in which the rates at different energies of the 
X-rays from the axion to photon conversion are included. Following this criterion, a 
systematic study of the background level for different software efficiencies has been 
performed. The F.O.M. from equation 5.3 has been computed using a binning of 
0.5 keV in the Range of Interest (Rol), from 2 — 7 keV. The results for the different 
detectors are shown in figure 5.7 and summarized in table 5.1. 


Detector 

e s (5.9 keV) 
% 

e s (3 keV) 
% 

FOM det 

c -1 / 2 keV 1 / 2 cm s 1 / 2 

Sunset 1 

75 

40 

106 

Sunset 2 

70 

35 

89 

Sunrise 

80 

50 

138 


Table 5.1: Software efficiencies used for the Micromegas data analysis during the 2011 
campaign. The selected software efficiencies at 3 and 5.9 keV are the ones that maximize 
the F.O.M.. 


From the results presented in table 5.1, one finds that the Sunrise detector has the 
higher F.O.M. and thus the better capabilities for axion discovery. On the other 
hand, Sunset detectors have a smaller F.O.M.. This behavior could be explained 
by their modest detector performance, presented in section 5.2. Also, the region 
which optimizes the F.O.M. seems shifted to lower software efficiencies because the 
background level rises for higher efficiencies. The software efficiencies have been 
selected taking into account that for the same F.O.M. the background level and 
thus the tracking counts are minimized, also the background spectra have to be 
coherent within its expected shape. 


5.4 2011 Micromegas background and tracking results 

Following the selection criterion described in section 4.4 and using the software effi¬ 
ciency presented in section 5.3.2, the background data of the Micromegas detectors 
working at CAST during the 2011 data taking campaign have been analyzed. 

Background and tracking data are usually stored in the same file. In order to 
distinguish its origin, the event time of the analyzed data is used. The tracking 
times can be extracted from the tracking variables and using the criteria described 
below. 
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Figure 5.7: Figure of Merit ( F.O.M .) of the different Micromegas detectors during the 
2011 data taking campaign. From top to bottom: Sunsetl, Sunset2 and Sunrise. The x- 
axis represents the software efficiency at 5.9 keV and the y-axis the software efficiency at 
3 keV. The different colors show the value of the F.O.M. computed using equation 5.3. 
The selected efficiencies for the final analysis are marked with a black ellipse. 

5.4.1 Tracking definition and data taking overview 

The axion signal would be an excess of counts while the magnet is pointing the 
Sun. So the discrimination of the background and tracking events is mandatory. 
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The CAST magnet is considered in tracking when the conditions described below 
are fulfilled: 

• The magnetic field is present in the magnet. 

• The corresponding gate valves (VT1, VT2 and VT3) are open. 

• The tracking program is in tracking mode and the Sun is reachable. 

• The precision of the vertical and horizontal angles is lower than 0.01°. 


Applying these conditions to the information generated by the tracking software 
and the slow control, the tracking times are extracted and background and tracking 
events can be separated. 


Detector 

Number of 
trackings 

Tracking 
time (h) 

Background 
time (h) 

Sunrise 

46 

72.8 

2149.7 

Sunset 1 

45 

72.0 

2045.1 

Sunset2 

45 

72.0 

2105.6 


Table 5.2: Summary of tracking and background times in the 2011 data taking campaign 
for the Micromegas detectors at CAST. 

The 2011 data taking campaign started on the 11 th of May and finished on the 22 nd 
of July. The magnet bores were filled with 3 He with a pressure from 82 mbar to 
108 mbar at 1.8 K, covering 172 pressure steps. However, there were some technical 
stops during the data taking due to mechanical problems which translate into a data 
taking efficiency of a 65% (93 trackings of 144 possible). The main incidents during 
2011 are described below: 

• 30.05.2011 - 06.06.2011 There was a crash noise in the movement of the magnet 
during the tracking. After some investigations it seems that a wheel was 
damaged. Because it was not a serious problem the data taking was restarted 
and the wheel was replaced at the end of the year. 

• 22.06.2011 - 02.07.2011 A power glitch in the experimental area caused a 
quench in the magnet, the vacuum and the detector system were also stopped. 

• 10.07.2011 - 16.07.2011 The data taking was stopped due to a problem in the 
cryogenics, which produced a rise of the magnet temperature. 


It is remarkable that there was not tracking data lost due to a detector problem or 
malfunction. 
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5.4.2 2011 Micromegas background and tracking data 

As it was described in the previous section tracking and background data are com¬ 
puted separately. In a first stage of the analysis tracking and background levels are 
extracted. A summary of the background and tracking data is shown in tables 5.3 
and 5.2. 


Detector 

Background level 
c cm~ 1 s 1 keV 1 

Tracking level 
c cm 1 s 1 keV 1 

Sunrise 

(6.09 ± 0.10)xl0 -6 

(5.71 ± 0.55)xl0 -6 

Sunset 1 

(5.96 ± 0.10)xl0 -6 

(6.14 ± 0.57)xl0~ 6 

Sunset2 

(6.83 ± 0.11)xl0 -6 

(7.58 ± 0.63) xl0“ 6 


Table 5.3: Tracking and background levels for the three Micromegas detectors during the 
2011 data taking campaign. The data have been computed from 2 — 7 keV and inside the 
colbore area (14.52 cm 2 ). 




Figure 5.8: Comparison between background (blue dots) and tracking (red dots) levels 
during the 2011 data taking campaign. The tracking level is computed for single trackings 
while 240 h have been computed in order to extract the background level. The error bars 
represent 1 a. 












































































5.4. 2011 Micromegas background and tracking results 


97 


From the results presented in table 5.3, it seems that background and tracking levels 
are compatible for all the detectors. In this case background and tracking levels have 
been compared using the mean value in the entire data taking period. However, it is 
more suitable to define a background level related to the corresponding tracking. For 
this purpose the background events have been computed ~120 h after and ~120 h 
before a given tracking. These results are shown in figure 5.8 in which background 
and tracking levels are compatible within its error. 



Energy (keV) 


Figure 5.9: Comparison between background (blue bars) and tracking (red bars) spectra 
during the 2011 data taking campaign. A binning of 0.5 keV has been selected in order to 
compute both distributions. 
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Figure 5.10: Spatial distribution of the background (blue dots) and tracking (red dots) 
events during the 2011 data taking campaign. Only the events inside the coldbore area and 
inside the Rol ([2 — 7] keV) are shown. 


Also, background and tracking spectra are compared in the 1 — 10 keV range (see 
figure 5.9), in which three different peaks can be distinguished: the 8 keV peak, 
caused by the fluorescence of the external particles on the copper anode readout; 
the 6.4 keV peak due to fluorescence in the stainless steel pipe and the 3 keV peak, 
which correspond to the fluorescence of the Argon in the detector. A more detailed 
study of these events will be described in chapter 7. 

Finally, the spatial distribution of the background and tracking events in the anode 
readout are represented in figure 5.10. In a close look to figure 5.10 there are small 
regions without events. It is related to some dead strips described in section 5.2.3. 
The Micromegas on the Sunrise side has the larger dead area. However, this re¬ 
gion is negligible in comparison with the coldbore area. Apart from these issues, 
background and tracking events are homogeneously distributed. 

As shown, background and tracking data are compatible in terms of rate, spectrum 
and spatial distribution. However, additional statistical processing of the data has 


















5.4. 2011 Micromegas background and tracking results 


99 



Figure 5.11: Time binning schema used to study the statistics compatibility of tracking 
(red) and background (blue) data. 


been performed in order to check the compatibility of the tracking and the back¬ 
ground. Being events with a low probability, its occurrence should follow the Poisson 
distribution. In order to verify it, background and tracking events are divided in 
time bins (see figure 5.11) and the distribution of the events inside the time bins is 
extracted. 


Sunset 1 Sunset 2 



Sunrise 



Figure 5.12: Distribution of the number of counts using a time binning of 15 minutes. The 
distribution has been obtained following the scheme from figure 5.11. 


The distribution of the background and tracking events in the time bins are shown 
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in figure 5.12, in which a binning of 15 minutes has been used. The black line is 
the expected Poissonian, that has been calculated using the mean value of the back¬ 
ground level, the blue dots correspond to background events and the red squares are 
the distribution of the tracking events with an errorbar of 1 a. For instance, tracking 
and background distributions seems to fit properly to the expected Poissonian. 
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6.1 Introduction 

In order to distinguish the presence of an axion signal, the data have been analyzed 
using an unbinned likelihood method. In the case of absence of signal, a limit on 
the axion to photon coupling is extracted. All these features will be described in 
this chapter. 

Moreover, the gas dynamics of the 3 He inside the magnet bores is not trivial and 
could affect the calculations of the axion mass and will be described. Finally, the 
constraints on the coupling limit for the 2011 CAST Micromegas data together with 
the limit of the 3 He phase will be presented. 

6.2 Unbinned likelihood method. 

In order to take in account the variations of the density inside the coldbore during the 
tracking, an unbinned likelihood method was developed for the statistical analysis 
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of the data of the 3 He phase. Before the explanation of the unbinned likelihood 
method, the likelihood function will be introduced. 


6.2.1 The likelihood function 

The likelihood function is given by the probability density function (p.d.f.) evaluated 
with the data y, as a function of parameters f(y;6). For a set of independent jp 
measurements, the joint p.d.f. factorizes and the likelihood function is given by: 


m = l[nvi\o) (6-i) 

i=l 


here L(9) is the likelihood function. For the CAST analysis the p.d.f is given by the 
Poisson distribution and thus, the likelihood function can be written as: 


1 






=n 


( 6 . 2 ) 


2—1 


2—1 


where Lo is the normalization factor, the index i is referred to the energy bin, m 
is the number of counts measured in the i bin and pi is the expected number of 
counts in the i bin, given by the addition of the background and signal counts, by 
the expression: 


Hi = bi + Si (6.3) 

here bi is the expected number of background counts and s t is the expected number of 
X-rays from the axion-photon conversion. The expected signal Sj has a dependence 
on the gas density p, axion mass m a and coupling constant g a7 and is given by: 

fE+AE ifo rE+AE j 

Si(g aj , m a , p) = J (iT P a ^~/ £ iAAt(lE = g A arj j -^AtdE (6.4) 

where is the solar axion flux from equation 1.31. P a is the conversion probabil¬ 
ity of the axion into a photon in a strong magnetic held, presented in equation 1.46. 
Ei is the detector efficiency for the bin i, in the case of Micromegas detectors quan¬ 
tum and software efficiencies have to be taken into account. A is the magnet bore 
area and At the time exposure. The integral is performed between E and E + A E, 
referred to the corresponding energy bin i. In the second identity of equation 6.4 the 
formula has been reduced. Indeed, the expected number of counts are proportional 
to and the expected number of signal counts per time and energy unit. 

Instead of using the natural likelihood, it is easier to work with the logarithm of the 
likelihood InL, also called log-likelihood. Maximizing the likelihood is equivalent to 
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minimize the quantity —21nL. According to Wilks’ theorem [119], if certain gen¬ 
eral conditions are satisfied, the minimum of —2 In L approaches a y 2 distribution. 
Computing the logarithm of equation 6.2, the following expression can be derived 1 : 

2 n 

- y = In L = y] nj - m + nj log — (6.5) 

i=1 1 

This equation is valid in order to calculate the log-likelihood with a fixed density 
inside the magnet bores and for a given axion mass. For instance, it was used for 
the analysis of the vacuum and the 4 He phase. However, for the 3 He phase an 
unbinned likelihood method was developed, in order to compute the variations of 
the gas density inside the magnet bores during the tracking. 

6.2.2 The unbinned likelihood 

In this case the likelihood is given by the product over the events, not the bins. The 
unbinned likelihood is calculated using infinitesimal time bins in the limit of zero 
counts contribution and one count contribution, by the expression: 

ln£= ’Y^L k {n i = 0) + y ^L k (nj = 1) (6.6) 

k k 

S -V-' S -V-' 

Zero counts cont. One count cont. 

here the index k is referred to time bins, which are short enough to have either 0 or 
1 count (m = 0 or m = 1). Indeed, equation 6.5 can be written as: 

ln£= ~ 22 +* +'22 (l-Ah) + ln Ah (6.7) 

k m ,7ii=0 k\rii =1 

s -v-' s -v-" 

Zero counts cont. One count cont. 

Using equation 6.3, the expression introduced before can be rewritten: 

In C = n c - 22 bi + 22 Si + 22 ln (Pi + Si) (6.8) 

k-,rii=0 k\rii =0 k\rii = l 

k\rii=l k\rii =1 

here n c is the number of tracking counts. The equation from above can be expanded 
applying equation 6.4: 


In £ = n c — 22 b i+9, 

k\rii =0 
k\rii= 1 


% 


""f dn, 


'try 

dE 


A tdE+ 22 

k-,rii =1 


dbi 

dt 


At + 5a 7 Ai 


rEi~\~AE 

nr <dE 

e, dE 


(6.9) 


1 Using the approach Inn! ~ nlnn — n 
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here At is the corresponding time bin and E, and Ef are the initial an final energy of 
the Rol used in the analysis ([2 — 7] keV in the Micromegas detectors). Some terms 
are not considered in the calculations, since they do not depend on the coupling 
constant. However, in order to perform a likelihood fit they have to be taken into 
account, because the resulting % 2 has to be compared with the number of degrees 
of freedom. Nevertheless, this fact does not affect the calculations performed in this 
work and those terms have been considered as a constant value, obtaining: 


f dria-y . 

— tdE+ 
dE 

( 6 . 10 ) 

— ^2 in At 

k;rii=l 

C = n c - E *- E in At (6.11) 

k]rii= 0 k]7ii =1 

k-rii = 1 



Note that the fourth term in equation 6.10 is the sum over the tracking counts n c . 
Finally, the unbinned likelihood can be reduced to the expression: 


X i r 4 

- y = ln £ = - 9 


0,7 


f Ef dn a ry f dbi 

/ —A7 A tdE + > in —- + a 

J E dE ^ \dt 

J . i =0 v 


4 

a7 


rEi+AE 


dn, 


a 7 , 


dE 


dE +C 


Expected axion events 


Counts over the background 


( 6 . 12 ) 


here the first term is proportional to the expected number of X-rays from the axion- 
photon conversion, while the second term is the sum of the number of expected 
background and signal counts over the detected tracking counts. In order to perform 
the analysis, the axion mass is fixed and a y 22 is obtained for several coupling 
constant values. Finally the presence or absence of signal can be evaluated in the 
X 2 -9a-y parameter space and in the case of absence of signal a coupling limit can be 
derived. These features will be described in section 6.4. 


The method described above is explained for the case of one detector. However, 
the analysis is performed by taking into account the data of all the detectors. The 
contribution of the different detectors can be added in the log-likelihood, by the 
expression: 


Atot — | Ed in C-tot — ^2 ln ^d 


(6.13) 


2 Although it is not a pure y 2 distribution, for simplicity this notation will be used in the 
document. 











6.3. He as buffer gas 


105 


here d is referred to the corresponding detector. Using this method the detectors 
can be treated separately until the last stage of the analysis. 

6.3 3 He as buffer gas 

In this section the dynamics of the 3 He inside the coldbore will be treated. There 
are mainly two different issues concerning the buffer gas: the gas density and the 
effective coherence length inside the magnet bores. Both of them are related with 
the density profile in the magnet bores during the tracking and will be introduced. 
On the other hand, the absorption coefficient of the X-rays inside the 3 He is related 
with the density inside the magnet bores and will be described. 

6.3.1 Gas density 

As it was presented in section 1.7.2.2 the effective photon mass and consequently 
the coherence condition is determined by the density of the buffer gas inside the 
magnet bores, given by equations 1.52 and 1.53. 

The gas density inside the magnet bores is calculated using the cold bore pressure, 
measured at the MRB side and the magnet temperature, by computing the equation 
of state (EoS) of the 3 He gas [72], During trackings the pressure changes continu¬ 
ously due to the variations of the hydrostatic pressure inside the magnet bores and 
convection effects close to the cold windows. These effects are related with the 3 He 
dynamics inside the magnet bores. Moreover, the magnet temperature has small 
changes of about 15 mK due to the tilt of the magnet. 

In order to take into account these effects, the density inside the magnet bores is 
calculated continuously. The pressure at the center of the magnet is determined 
from the measured pressure and the hydrostatic pressure difference. The density 
in the center is calculated from the central pressure and magnet temperature using 
the EoS of 3 He. In this way, the dynamic effects of the 3 He inside the magnet 
bores are taken into account by using the measured pressure. The variations of the 
magnet temperature, pressure and the calculated density during a single tracking 
are represented in figure6.1. 

6.3.2 Effective coherence length 

The axion-photon conversion probability is computed using the magnet length (see 
equation 1.46). However, the dynamics of the 3 He inside the magnet bores may cause 
a reduction of the conversion region. This problem was solved by the introduction 
of an effective coherence length and the axion to photon conversion can be properly 
parameterized. 

The effective coherence length was introduced in order to take account the variations 
of the density inside the magnet bores. The density profile is not a constant value 
along the cold-bore. This effect has been included in the analysis by introducing 
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Time 


Figure 6.1: Evolution of different observables during a tracking: pressure inside the cold 
bore (black line), temperature of the magnet (blue line) and calculated density (red line). 
The data correspond to the morning tracking on the \2 th of May of 2011. 


the effective coherence length, that could be defined as the region inside the magnet 
bores in which the density has constant value. Moreover, the problem is more com¬ 
plex when the tilting of the magnet is introduced. Since the density profile cannot 
be measured directly, it has to be calculated by Computational Fluid Dynamics 
(CFD) simulations. 




Figure 6.2: Left: Density profile inside the magnet bores at a pressure of 83 mbar (red 
Line). The dashed blue lines represent the allowed region for the density while the dashed 
black lines mark the effective length. Right: Effective coherence length as a function of 
the cold bore pressure. The black dots are the result of the CFD simulations and the blue 
line is a linear fit. 
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The CFD simulations take account of the geometry of the magnet including the 
cold windows and link volumes together with physical phenomena such as gravity, 
convection and the EoS of the buffer gas. A detailed study of the density profile 
using the CFD simulations at different pressures and tilting was performed and 
the effective length of the magnet has been parameterized [120]. In this case the 
coherence criterion is that the variation of density along magnet bore is less than 
A p < 0.001 kg m -3 [72, 121]. 

The density profiles and thus the effective coherence length are extracted from the 
CFD simulation results with the magnet in a horizontal position. For this purpose 
several CFD simulations at different pressures were computed. The effective coher¬ 
ence length is defined as the region in which the deviations from the central density 
are below A p < 0.001 kg m -3 (see figure 6.2 left). Finally, the results of the differ¬ 
ent simulations were fitted to a line (see figure 6.2 right). So the effective coherence 
length is parameterized using the coldbore pressure, by the expression: 


L eff = 8.611 - 0.021Pcb (6.14) 

here T e // is the effective coherence length in meters and Pcb is the cold bore 
pressure in mbar. The dependency of the effective length with the cold bore pressure 
has been included in the statistical analysis of the data. 

6.3.3 X-rays absorption in ’He 

The absorption coefficient T of the X-rays in the buffer gas is included in the conver¬ 
sion probability, given by equation 1.46. This term takes account of the X-rays from 
the axion to photon conversion that are absorbed in the buffer gas. The absorption 
of X-rays in a medium is given by: 

l = e“ rL ; r = HP (6.15) 

to 

here ^ is the fraction of X-rays that are absorbed in the medium and L is the 
magnet length. The absorption coefficient F is parameterized by the attenuation 
coefficient p, and the gas density p. The attenuation coefficients at different energies 
have been extracted from the NIST 3 X-COM database and can be parameterized 
by the expression: 


fi(E) = exp (-1.5832 + 5.9195e“ a353808£: + 4.03598e“°' 970557i? ) (6.16) 

The agreement between the NIST X-COM data and the fit is shown in figure 6.3. 
For the statistical analysis of the data, the absorption coefficient F is calculated 
using equation 6.16 and the gas density is estimated with the method presented in 
section 6.3.1. 

3 National Institute of Standards and Technology 
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Figure 6.3: Attenuation coefficient as a function of the energy for He. The red circles 
show the data extracted from the NIST X-COM database, while the black line is the fit 
result from equation 6.16. 

6.4 Coupling limit for the 2011 Micromegas data 

In this section a limit on the axion to photon coupling for the 2011 Micromegas data 
will be derived. In a first stage, the axion mass coverage will be presented. Then the 
method to distinguish the presence or absence of an axion signal will be described. 
Finally, an upper limit on the axion to photon coupling will be extracted. 

6.4.1 Axion mass coverage 

The axion mass coverage can be estimated dividing the exposure time in density 
steps. However, it is more accurate the use of the expected number of signal counts 
ria-y- It takes into account the coherence process in the axion to photon coupling 
and the efficiencies of the detectors. The n a7 is calculated using the first term of 
equation 6.12, derived from equation 6.4 and given by: 

r E f r E f (i n 

n ai (gay,rn a ,p)= —^P a ^ 1 e(E)AAt(p)dE = g^j — t(p)dE (6.17) 

J E-i J E^ 

The expected number of counts is calculated for a given axion mass using the density 
exposure of the detectors. Figure 6.4 shows the number of expected counts during 
2011 data taking campaign for the different Micromegas detectors and the addition 
of them, together with the density exposure used for the calculations. 
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m a (eV) 


Figure 6.4: Top: Density exposure of the Micromegas detectors during 2011. Different 
detectors are plotted: Sunrise (blue line), Sunset (red line) and the addition of all of them 
(black line). Note that only one Sunset detector has been drawn because the exposure is 
the same in both Sunset. Bottom: Number of expected counts as a function of the axion 
mass, equation 6.17 has been used in the calculations, assuming a coupling constant of 
g ai = 10’ 10 GeV -1 . The three detectors have been drawn: Sunrise (blue line), Sunset 1 
(red line) and Sunset2 (green line). The black line is the number of signal counts calculated 
by the combination of all the detectors. 

6.4.2 Search for an axion signature 

As it was described in section 6.2 the minimum of —21nL follows a x 2 distribution. 
This feature can be used to estimate the most probable value of g a7 that is located 
at the minimum of the x 2 — Sa-y distribution. Moreover, the change in x 2 by one 
unit corresponds to one standard deviation (er) shift in the parameter estimation. 

In order to study the presence or absence of an axion signal, the unbinned likelihood 
is performed for several axion masses in the scanned region. For a given axion mass, 
the x 2 — di'y distribution is extracted and the minimum is located. Later on, the 
region surrounding the minimum is scanned and the values of g ^ at one x 2 unit 
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above the minimum are extracted. These values correspond to one a deviation from 
the most probable value. This analysis method is shown in figure 6.5. 
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Figure 6.5: \' 2 as a function of gf n for an axion mass of m a = 1.061 eV (black line). The 
dashed blue line marks the distance at one x 2 unit to the minimum. The red lines represent 
the positions of the points at one a. Note that the y-axis is not a pure y 2 distribution. 
However this fact does not affect the calculations. 


Using this method the standard deviation of the most probable value of g\ y has 
been calculated for several axion masses in the scanned region. Although the x 2 
is supposed to be a parabola in the minimum, the standard deviation from the 
X 2 extracted from the likelihood seems to be asymmetric. And thus, two different 
standard deviations are taken into account, the one on the left 07 and on the right 
a r from the minimum. As is shown in figure 6.5, some points lies in negative values 
of which does not have physical meaning, however they are significant from the 
statistical point of view. 

In absence of an axion signal, the minimum of \' 2 has to be compatible with g ay = 0. 
In order to verify it, the deviation of the minimum in er’s units has been calculated 
for several axion masses and are shown in figure 6 . 6 . 

From the results presented in figure 6.6 we conclude that there is no axion signa¬ 
ture, because the deviation is below 2 cr’s. In order to see the significance of these 
deviations, a set of data has been simulated. Obtaining a deviation of up to 2 <r’s 
only due to statistical fluctuations of the tracking counts [122], On the other hand, 
deviations above 3 it’s would require a dedicated analysis. 
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Figure 6 . 6 : Deviation in a units of the minimum of \ 2 from zero. The distributions with 
a minimum below zero are not represented and thus the sigma on the left 07 has been used 
in the calculations. 


6.4.3 Limit on the axion-photon coupling 

After discarding a possible axion signal, a limit on the axion to photon coupling can 
be extracted. An upper limit on the coupling constant for a given axion mass is 
calculated by integrating the Bayesian posterior probability from zero up to a 95% 
with a flat prior in g* , by the expression: 


/cf * 7 e +c dg ai _ £? Jo ai e "2 dg ai 
Jo°° e~^ +c dg ai e~^dg ai 


(6.18) 


In the practice e - *"^ 2 is integrated in a wide range starting from zero and then the 
which correspond to a 95% of the integral is extracted. This analysis method is 
shown in figure 6.7. 

Using the method presented before, a limit on the axion to photon coupling for 
the Micromegas detectors during the 2011 data taking campaign has been obtained. 
The different coupling limits for all the detectors separately and combined are shown 
in figure 6 . 8 . Also, for the scanned region during 2011, the one where the coherence 
condition is fulfilled, an average value of the coupling constant has been extracted: 


g ai < 3.90 x 10 - 10 GeV -1 for 1 < m a < 1.17eV (6.19) 


at 95% of confidence level (C.L.). This limit could be improved including the data 
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g 4 , (x 10 40 GeV 4 ) 

Figure 6.7: Calculation of the coupling limit, the black line represent the e - * 2 / 2 — gf n 
distribution for an axion mass of m a = 1.061 eV. The blue filled area is the integral from 
zero up to a 95%. The end of the tilled area marks the value of the coupling limit. 



Figure 6.8: Axion to photon coupling upper limit for the different detectors during the 
2011 data taking campaign. The solid blue line corresponds to the Sunrise detector, the 
solid red and green lines show the limit for Sunset 1 and Sunset2 detectors respectively. 
The solid black line is the combination of all the detectors and the dashed brown line 
represents the KSVZ model with E/N = 0. 


of the CCD detector, working at CAST during 2011. 
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6.4.3.1 Estimation of the systematic error on the 2011 coupling limit 

The determination of the coherence length is considered of the main source of error 
for the calculations of the systematics on the coupling limit for the 2011 data. 
Indeed, the profile density along the magnet cannot be measured and the CFD 
simulations have some uncertainties. Moreover, the effective coherence length could 
be affected by the tilting of the magnet and be reduced. On the other hand a more 
optimistic scenario could be considered, neglecting these effects and using the entire 
magnet length in the calculations. 

In order to estimate the systematic uncertainties on the coupling limit presented 
before, two extreme cases have been taken into account: the first one by considering 
all the magnet length of T e // = L, mag = 9.26 m as the coherence region and the 
second one by using a conservative scenario, substituting the formula for the effective 
length (from equation 6.14) by L e ff = 7.611 — 0.021 Pcb- These two cases have been 
implemented in the analysis and a coupling limit has been extracted, the results are 
shown in figure 6.9. 



Figure 6.9: Calculation of the systematic error for the 2011 Micromegas data. Two 
different coupling limits have been derived: an optimistic case by considering the total 
length of L e ff = 9.26 m (red line) and a pessimistic case using an effective length given 
by L e ff = 7.611 — 0.021 Pcb (blue line). The black line is the coupling limit presented 
before. 


The averages values of the two different cases have been calculated in the axion 
mass range from 1 < m a <1.17 eV. Obtaining systematic uncertainties below 7% 
in contrast to the previous result: 


g ai < 4.04 x 10 -10 GeV -1 for L eff = 7.611 - 0.0211^5 m 

9a-y < 3.65 x 10~ 10 GeV _1 for L e ff = 9.26 m 


( 6 . 20 ) 
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6.5 Coupling limit for the 3 He phase 

The data of the 3 He phase of the CAST experiment have been computed and a 
coupling limit has been extracted. These results correspond to the data taking 
campaigns from 2008 to 2011, which have been already published in [71] and more 
recently in [72], Although the data of all the operative Micromegas detectors have 
been taken into account in order to obtain the limit, the data of the CCD-telescope 
system, from 2009 to 2011, have not been included. The results of the 3 He phase 
are shown in figures 6.10 and 6.11. 



m. 


/eV) 


Figure 6.10: Parameter space m^—ga-y excluded by the CAST experiment for the vacuum 
and 4 He phase (black line) and the latest results of the 3 He period (red line). The yellow 
band is the most favored region for the axion models, while the green solid line corresponds 
to the KSVZ model ( E/N = 0). Plot taken from [72]. 

During the 3 He phase CAST extended its previous limit towards higher axion 
masses, obtaining an average value on the coupling constant of: 

g ai < 2.94 x 10 _10 GeV _1 for 0.37 < m a < 1.17eV (6.21) 

On the other hand, the limit shown in figure 6.11 has a characteristic shape. The 
high-frequency structure is because the different axion masses have a different ex- 
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Figure 6.11: Coupling limit for the CAST 3 He phase at a 95% of C.L., corresponding to 
an axion mass range of 0.37 < m a < 1.17 eV (solid red line). The green line is the KSVZ 
benchmark model ( E/N = 0). 

posure times, also statistical fluctuations of the data may cause this behavior. The 
worsening of the coupling limit for higher axion masses can be explained by different 
reasons: the time exposure was reduced to the half for axions masses above 1 eV; 
the continuous decrease of the effective coherence length and the increase of the 
absorption coefficient T with the density, also the data from the CCD-telescope for 
axion masses above 0.64 eV have not been not included in the analysis. 

CAST finished the 3 He phase extending the axion search up to m a = 1-17 eV, 
being the first time that an helioscope crosses the KSVZ ( E/N = 0) model, one 
of the most favored ones. Although there was no axion signal, CAST currently is 
rescanning the vacuum phase with an improved sensitivity. However, the goal is to 
scan more favored regions for axions and ALPs, this achievement cannot be reached 
with the existing CAST apparatus. It will require a dedicated magnet, optics and 
improved detectors, as it is proposed for IAXO, that will be described in chapter 8. 
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7.1 Introduction 

Low background X-ray detectors are mandatory for axion searches or more gener- 
ic.ally for rare event searches experiments. In this chapter the low background 
techniques developed in order to reduce the background level of the Micromegas 
detectors will be introduced. 

Moreover, the different set-ups for the Micromegas detectors together with the sim¬ 
ulations, developed in order to study the different contributions of the background, 
will be presented. These special set-ups were crucial in order to reduce the back¬ 
ground level in the Micromegas detectors in the CAST experiment. Finally, the 
different upgrades of the CAST Micromegas detectors and the background reduc¬ 
tion will be described. 
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7.2 Low background techniques 

CAST microbulk Micromegas exploit different strategies developed for the back¬ 
ground reduction: the intrinsic radiopurity of the detectors; the detector perfor¬ 
mance, which is related with the improvements on the manufacturing process; the 
event discrimination of the background events in which the front end electronics 
has an important role and finally the shielding strategy. All these techniques have 
been developed in the context of the TREX project by the University of Zaragoza 
group [123] and will be described in this section. 

7.2.1 Radiopurity 

The intrinsic radioactivity of the materials used in the construction of the detector 
may be a source of background. As it was presented in section 2.4.3.3, the microbulk 
Micromegas readout is made of kapton and copper, two radiopure materials. The 
intrinsic radiopurity of the microbulk Micromegas readout has been measured in 
the Canfranc Underground Laboratory (LSC) with a high purity Ge detector [124], 
These results have been confirmed recently by the measurements performed by the 
BiPo [126] detector at the LSC. 

On the other hand, the materials used in the construction of the Micromegas cham¬ 
ber may have an impact on the background. The chamber design was presented 
in section 3.3, the body of the chamber is made of Plexiglas which is intrinsically 
radiopure. However, the strongback is made of aluminum which has a non negligible 
contribution of 238 U and 232 Th [125]. The impact of the aluminum strongback in 
the background was measured at the LSC in an special set up, this contribution 
will be discussed in section 7.3.1. Moreover, an on-going program is measuring the 
contribution to the background of the different materials in the set-up (screws, gas 
gaskets, connectors, electronics, etc.). 

The thorough study of the intrinsic radiopurity of the detectors and the chamber 
lead to a new detector and chamber design which was installed at CAST during 2014 
in the Sunrise side. The features of this set-up will be described in section 7.4.2. 

7.2.2 Manufacturing technology 

The performance of the Micromegas detectors may have an impact on the back¬ 
ground and consequently in the sensitivity of the detectors. Indeed, as it was 
presented in section 5.2, during 2011 the Sunrise detector shows clearly a better 
performance than the Sunset ones. So the capabilities of the Sunrise Micromegas 
for axion discovery were higher. 

It is remarkable the progressive improvement of the performance of the Micromegas 
detectors that is closely linked to the development of the Micromegas technology. 
CAST has been one of the most demanding test-bench for Micromegas detectors 
from the classical Micromegas until the novel microbulk technology. Also, a sys- 
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tematic study on the microbulk Micromegas mesh and anode readout has been 
performed. This has led to the optimization of the different parts of the structure of 
the detectors, such as: the anode layout, the mesh pattern, the size and shape of the 
mesh holes and the pitch of the strips anode. This continuous improvement in the 
manufacturing process of the microbulk technology is reflected in the latest detec¬ 
tor design for the Sunrise Micromegas, the one with the best detector performance 
working at CAST, which will be presented in section 7.4.2. 

7.2.3 Event discrimination 

The algorithms implemented in order to discriminate X-ray like events have been 
improved since the beginning of the CAST experiment. The power of this discrim¬ 
ination is highly coupled with the quality of the readout, so improvements in the 
readout design or in the manufacturing process lead to improvements in discrimina¬ 
tion power. 

On the other hand, the rejection capabilities can be improved with the upgrade of the 
front-end electronics and the acquisition system. For this reason the novel front-end 
AFTER [108] electronics was installed at CAST during 2013 for all the Micromegas 
detectors. The main advantage of the new electronics is that every strip pulse is 
digitized and stored, in contrast with the Gassiplex-based electronics in which only 
the integrated charge was stored. The event information in the AFTER electronics 
is displayed in figure 7.1. 

The upgrade of the front end electronics in the strip readout leads to the modifi¬ 
cation of the electronic chain of the Micromegas detectors described in chapter 3. 
Also, a new acquisition software has been developed, based on C++ and ROOT 
libraries. Moreover, a new software analysis has been written, in this case the pulse 
shape analysis has been extended to every single strip and yields the definition of 
more parameters related to the strips pulses. Additionally, the cluster analysis has 
been extended to the z direction, which allows the definition of new observables: 
cluster size Z and sigma Z. These new parameters can be easily included in the 
discrimination method introduced in section 4.4. The upgrade of the front end elec¬ 
tronics in the strip readout supposed a reduction of the background level of about 
a ~25% in the Micromegas detectors at CAST. 

7.2.4 Shielding 

The shielding design of the Micromegas detectors during 2011 at CAST has been 
described in sections 3.3.1 and 3.3.2. The shielding was made of an innermost layer 
of copper 5 mm thick, 25 mm of lead shielding and a 2 mm thick cadmium sheet. 
After a thorough study of the different background sources from simulations (see 
section 7.3.3) and some measurements in special set-ups (see section 7.3.1 and 7.3.2), 
a shielding upgrade was proposed. 

The shielding strategy is to increase the thickness of the different copper and lead 
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PositionX (strips) 


PositionY (strips) 


Figure 7.1: Top: Pulses from a 5.9 keV calibration event digitized in the AFTER elec¬ 
tronics. Bottom: 2D view of a background event acquired with the AFTER electronics in 
which the two different projections XZ (left) and YZ (right) are shown. 

layers in the set-up. However, in the CAST experiment there are space and weight 
constraints on the magnet movable platform. Taking into account these constraints 
a new shielding design was projected with 10 mm of copper and 100 mm of external 
lead. Also, the cadmium sheet was removed because the contribution of neutrons is 
negligible. Moreover, the stainless steel pipe to the magnet bore has been replaced 
by a copper one due to the low intrinsic radioactivity of the copper and to avoid the 
stainless steel fluorescence from 5 — 7 keV, that is inside the Rol. Also, a PTFE 1 
cylinder is installed inside the pipe in order to attenuate the copper fluorescence. 

Finally, an active muon veto was installed in order to distinguish events related with 
cosmic muons. Although muons that interact directly with the detector are easily 

1 Polytetrafluoroethylene 
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rejected by the analysis, cosmic muons may provoke fluorescence in the surrounding 
materials, which can contribute to the background. 


7.3 Test benches and simulations 

The CAST Micromegas set-up described in sections 3.3.1 and 3.3.2 has been repli¬ 
cated and some measurements were performed underground at the LSC 2 and at 
surface level in the laboratories of the University of Zaragoza. The main purpose 
of these test benches is to determine the different contributions to the background. 
The measurements performed in special set-ups together with the simulations were 
crucial for the upgrade of the shielding of the Micromegas detectors at CAST and 
will be detailed. 


7.3.1 Underground measurements 

A CAST detector replica was installed at the LSC in order to understand the ori¬ 
gin of the background in the Micromegas detectors. The Canfranc Underground 
Laboratory (LSC) is situated at Canfranc (Huesca) in the Spanish Pyrenees with a 
depth of 2500 m.w.e. 3 . Here the muon flux is reduced by a factor 10 4 relative to the 
surface level [127]. So in this set-up we may consider negligible the contribution of 
the muons to the final background level. However, in contrast to surface operation, 
there are other sources of background that have to be taken into account, such as: 
the concentration of 222 Rn, that is several times higher than in surface and the 
environmental 7 and neutron flux that could be considerably larger. 



Figure 7.2: Left: CAST Micromegas replica at the LSC, in which the different parts 
inside the Faraday cage are labeled. Right: Set-up with 100 mm of external lead, in which 
the electronics and the dewar used to flush nitrogen are shown. 

Different measurements have been performed in the set-up at the LSC in which 
nitrogen is flushed to reduce the environmental 222 Rn level inside the shield- 

2 Laboratorio Subterraneo de Canfranc 

3 meter water equivalent 
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mg. In a first stage the CAST Micromegas shielding (5 mm of copper and 
25 mm of lead) was installed (see figure 7.2 left), obtaining a background level 
of ~ 6 X ICC 6 c cm~ 2 keV -1 s -1 , the same level reached at surface level with 
the same set-up. Later on, the lead shielding was increased to 200 mm and 
100 mm (see figure 7.2 right). These upgrades diminished the background to 
~ 2 x 10- 7 c cm -2 keV” 1 s” 1 [128], the lowest level reached with a Micromegas 
detector. 

Moreover, different contributions to the background have been measured at the LSC, 
such as the aluminum strongback and the 222 Rn (see figure 7.3). In order to see the 
impact of the aluminum cathode in the background, two different measurements 
were performed, one with a radiopure copper cathode and another with the alu¬ 
minum strongback, obtaining a difference of (5.2 ± 1.2) x 10” 7 c cm” 2 keV” 1 s” 1 
between them. In addition the nitrogen flow was stopped in order to see 
the effect of the environmental 222 Rn in our set-up, measuring a value of 
(3.0 ± 0.8) X 10” 8 c cm” 2 keV” 1 s” 1 per Bq/m 3 of air-borne 222 Rn in the 
surrounding atmosphere. This level has been quantified by measuring the 222 Rn 
concentration inside the Faraday cage with an alphaGUARD detector [130]. 



Figure 7.3: Background level in the LSC during 2011. An ultra-low background level (red 
aiTows) was interrupted due to the measurements of the aluminum cathode (blue arrow) and 
the 222 Rn intrusion (green arrow). Plot taken from [123]. 


The set-up at the LSC continues taking data, currently is measuring the contribu¬ 
tions of the non-radiopure materials (viton gaskets, brass gas connectors,...) close to 
the detector, which have been replaced by radiopure ones. These upgrades have led 
to a low limit closer to 10” 7 c cm” 2 keV” 1 s” 1 [123]. The set-up at Canfranc is an 
essential test bench for the Micromegas detectors in order to measure the intrinsic 
radiopurity of the materials surrounding the detector. 
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7.3.2 Surface measurements 

The measurements performed in the Zaragoza laboratory were a crucial key in order 
to verify the impact of the muons in the background. Moreover, this set-up was an 
important test bench to develop the DAQ for the new AFTER front end electronics. 

The most significant measurements performed at Zaragoza were the ones made with 
a replica of the shielding upgrade of Sunset Micromegas proposed for CAST, that will 
be described in section 7.4.1. From the results obtained at the LSC, an important 
background reduction was expected. After the first measurements in this set-up the 
background level was ~ 2 X ICC 6 c cm” 2 keV” 1 s -1 , one order of magnitude higher 
than the level reached at the LSC in the same conditions. Indeed, it indicates that 
the muons, that are highly suppressed underground, have a contribution in surface 
operation. 



Figure 7.4: Left: Set-up in the Zaragoza lab, both plastic scintillators for the muon 
rejection are shown. Right: Time difference between the delayed Micromegas trigger and 
the scintillator. The black dashed lines delimit the region for the rejection of events induced 
by muons in the off-line analysis. The shape of the distribution is due to the drift velocity 
in the Micromegas detectors. 

In order to see the effect of the muons in the set-up, one plastic scintillator was 
installed on the top and another in one side of the Micromegas shielding (see figure 
7.4 left). The signals from the scintillators are sent to a quad discriminator and the 
time difference between the delayed trigger from the Micromegas and the trigger 
from the scintillators is stored in a VME scaler CAEN V560. Computing this time 
difference, events in coincidence with a muon can be easily rejected (see figure 7.4 
right). 

After the subtraction of the events related with muons, the background level dimin¬ 
ished to ~ 1 x 10” 6 c cm” 2 keV” 1 s” 1 [129]. The resulting background spectrum is 
shown in figure 7.5, is remarkable the reduction of the 8 keV peak from the copper 
fluorescence. It may indicate fluorescences produced by the muons or muon showers, 
in a region surrounding the detector. The excellent results obtained in this set-up 
motivate the installation of a dedicated muon veto for the Micromegas detectors at 
CAST. 
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Figure 7.5: Background spectra from the measurements performed in the Zaragoza lab 
using the muon scintillator. The blue line corresponds to the level before the muon sub¬ 
traction, while the red line is the final level after the veto cut. 

The contribution of the muons to the background is not fully understood because 
the resulting background at surface level after the veto subtraction is significantly 
higher than the one measured underground. From the set-up shown in the left part 
of figure 7.4 the estimated veto coverage is about a ~ 95%, which is the fraction of 
muons that crosses the scintillators divided by the fraction of muons which crosses 
the chamber. The resulting background level might be induced by muons crossing 
the shielding or the pipe, or either far away from the detector. This effect could be 
minimized by extending the surface area of the scintillator. 

7.3.3 Simulations 

The motivation of the simulations is to understand the different contributions to 
the background in the Micromegas detectors. For this purpose the complete CAST 
set-up geometry has been implemented in the Geant4 toolkit. After simulating 
the physical events using the RESTSoft package [114], a code developed at the 
University of Zaragoza, the drift and the diffusion in the chamber are simulated 
using the parameters extracted from Magboltz [103]. Finally, the electronic response 
of the detector is implemented and the simulated data are transformed into the 
experimental data format. 

The main purpose of the simulations was to see the effect of the external 7 ’s in 
our set-up. And thus the environmental 7 flux was measured by a Nal detector 
in the CAST area. However, after computing the simulation, the background level 
obtained was significantly lower than the experimental background level, although 
the energy spectrum matches qualitatively with the one obtained experimentally. 
Moreover simulations show that most of the background events were generated by 
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gammas passing through the shielding outlets instead of penetrating through the 
shielding. Indeed, many of the final events are accumulated in the 5 — 7 keV region 
due to the fluorescence in the stainless steel pipe (see figure 7.6 left). While the 
8 keV peak is generated by the copper fluorescence in the anode readout, produced 
by the events crossing the detector. 



Energy (keV) 


Figure 7.6: Left: Background spectrum extracted from the simulations of the environ¬ 
mental 7 flux in the experimental CAST area. The contribution of the events coming from 
the cathode window or pipe is superposed. The spectrum has the same shape as the ones 
shown in figure 5.2. Right: Comparison between experimental background levels as a 
function of the lead shielding thickness in special set-ups and simulations. The black line 
corresponds to a fit for the simulated data (black circles) of environmental y’s. The red 
squares are the measurements performed underground while the blue squares arc the mea¬ 
surements in surface with a muon veto. Plots taken from [123]. 


As is shown in the right part of figure 7.6 the background level obtained from sim¬ 
ulations does not match with the surface measurements, while it fits properly with 
the measurements performed underground. It suggests that muons have a negative 
effect in our set-up, hypothesis that was confirmed by the surface measurements 
described in section 7.3.2. 

Due to the results of the simulations together with the measurements performed in 
the special set-ups, different upgrades were proposed for the Micromegas detectors 
at CAST: 


• The lead shielding has to be extended along the pipe to the magnet in order 
to avoid Compton processes in the pipe. 

• The stainless steel pipe has to be replaced by a copper one, in order to avoid 
the stainless steel fluorescence which is inside the Rol. 

• The use of an active muon veto is mandatory in order to minimize the negative 
effect of the muons in our set-up. 
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Following this roadmap, Sunset and Sunrise Micromegas detectors have been up¬ 
graded at CAST, resulting in a reduction of the background level of the detectors 
in a factor ~ 6. All these features will be described in the following section. 


7.4 CAST Micromegas upgrades: State of art 

Although CAST finished its research program in 2011 scanning axion masses up 
to 1.17 eV, a new 4 He and vacuum phase were proposed, partially motivated by 
the R&D in the background reduction on the Micromegas detectors. During 2012 
the Sunset Micromegas were upgraded with a novel shielding design and newly 
manufactured detectors. In these conditions the 4 He phase was rescanned in a 
narrow axion mass range m a ~ 0.4 eV [77]. In 2013 a new vacuum phase started, 
in which an improved sensitivity is expected due to the reduction of the background 
level of the detectors. During 2014 a dedicated X-ray optic was manufactured with 
a novel Micromegas detector in its focal plane at the Sunrise side. 

The different upgrades on the Micromegas detectors at CAST are the result of 
the low background techniques presented before. And thus the peculiarities of the 
Sunset and Sunrise upgrades will be described in this section. 

7.4.1 Sunset Micromegas upgrade 

Although the Sunset Micromegas shielding was upgraded during 2012, the set-up 
that will be described in this section is related with the 2013 and 2014 data taking 
campaigns. Even if during 2012 a veto was installed, its coverage was poor due to 
geometrical constraints in the experimental area, for this reason during 2013 two 
plastic scintillators was specifically constructed for this purpose. 



Figure 7.7: Left: Photo of the Sunset shielding partially opened in which the inner copper 
layer and the external lead are visible. Right: Scheme of the Sunset Micromegas shielding 
design in which the position of the plastic scintillators can be observed. 


The upgrade is focused on reducing the contribution of the external j’s and in 
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particular the steel fluorescence induced in the pipes. The lead shielding thickness 
was increased from 25 mm to 100 mm and the design is more compact, improving 
the shielding around the pipes to the magnet (see figure 7.7 left). The inner copper 
shielding has been increased from 5 mm to 10 mm in order to attenuate the Pb 
fluorescence and also the 45.6 keV 7 line from the 210 Pb. 



Figure 7.8: Micromegas background spectrum during the 2013 data taking campaign in 
the Sunset 1 (top) and Sunset2 (bottom) detectors. The blue line corresponds to the raw 
spectrum before the veto cut, the red line is the final background after the veto subtraction, 
while the green line is the contribution of the muons. Plots taken from [131]. 

Moreover, the connection to the magnet bores is done by a 10 mm thick copper pipe, 
which has an inner PTFE coating with a thickness of 2.5 mm in order to attenuate 
the 8 keV copper fluorescence peak. Additionally, the aluminum strongback has 
been replaced by a more radiopure copper one and all the components close to 
the detector have been carefully selected and cleaned. In order to discriminate 
events induced by muons, two plastic scintillators have been manufactured. Due 
to geometrical restrictions in the experimental area one scintillator is placed on the 
top of the shielding covering the detectors projection, while a larger one is placed 
on the back of the shielding (see figure 7.7 right). The coverage to cosmic muons is 
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estimated to be a ~ 95% [129]. 

In addition, the Sunset Micromegas detectors were replaced by two new ones spe¬ 
cially manufactured for this purpose. Also, during 2013 the Gassiplex cards were 
replaced by the novel AFTER front-end electronics for the strips readout. This 
upgrade led to a reduction of the background level in a factor 4 [129], 

After the implementation of all the novelties described before, the background level 
in the Sunset Micromegas detectors diminished to ~ 1 X 10 -6 c cm -2 keV -1 s -1 
[131]. It is remarkable the reduction due to the muon veto that can take account 
of a 50% of the background events after the analysis. The background spectra of 
both Sunset Micromegas detectors are shown in figure 7.8. In contrast with the 
background of the 2011 data taking campaign (see figure 5.9), the steel fluorescence 
([5 — 7] keV) has disappeared. However, the final background is dominated by the 
copper fluorescence at 8 keV and its escape peak at 5 keV. 

7.4.2 The new Sunrise Micromegas + XRT system 

During the 2014 CAST data taking campaign, a new X-ray focusing device was 
installed in the Sunrise side with a Micromegas detector in its focal plane. Being 
the first time that an X-ray optic is specifically built for axion research. Moreover, 
the Sunrise Micromegas detector has a novel design which collects all the R&D in 
low background techniques for the Micromegas detectors. 

The X-ray telescope (XRT) has been designed and built by the groups of LLNL 1 , 
DTU° and the University of Columbia. It is made of segmented glass substrates with 
13 nested layers covering an area of 14.52 cm 2 and has been manufactured using 
the same techniques developed for the NASA’s NuSTAR 4 * 6 [132] satellite mission. 
The XRT might improve the effective background of the Micromegas detector in 
a factor ~50, while reducing the efficiency in a factor ~2. The new XRT and the 
Micromegas detector were installed and aligned with the CAST magnet between the 
25 th August and the 4 ift September of 2014, the final set-up is shown in figure 7.9. 

A new Micromegas detector has been designed for the 2013 data taking campaign. It 
is the prime example of the current state of art in low background techniques for the 
Micromegas detectors. In contrast with the previous design described in section 3.3, 
the body and the chamber of the detector is made of 20 mm thick radiopure copper 
and all the gaskets are made of PTFE. Also, a new held shaper has been designed, 
printed on a kapton circuit and integrated in the chamber. It makes more uniform 
the drift held and reduces the border effects, also is covered by a 2 mm thick PTFE 
coating in order to avoid the copper fluorescence. The high voltage connections 
were implemented in the detector printed board, which allows an easy extraction of 
signals and voltages from the shielding. The design of the chamber and a photo of 

4 Lawrence Livermore National Laboratory 

’Technical University of Denmark 

6 The Nuclear Spectroscopic Telescope Array 
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Figure 7.9: Photo of the new Sunrise Micromegas + XRT system in the CAST experi¬ 
ment. The different parts of the set-up are labeled. 


the new Micromegas detector are shown in figure 7.10. 

Following the Sunset design, a copper pipe interface with a PTFE coating has been 
installed. However, in this case the aperture of the pipe has been reduced from 43 to 
20 mm of diameter. In this case the expected signal area is considerably smaller due 
to the focusing of the XRT, which allows a better coverage around the pipe. Since 
the body and the chamber are made of radiopure cooper, the inner shielding is not 
necessary, also it shields among the background events coming from the electronics 
through the raquette, a weak point in the Sunset design (see figure 7.11). The 
external shielding is made of 100 mm of lead. However, in some places the lead 
thickness is smaller (~70 mm) due to geometrical constraints in the experimental 
area. Finally, a plastic scintillator is installed on the top of the shielding. Although 
only one scintillator is installed due to the spatial limitations, the muon veto has 
been extended through the pipe in order to discriminate events induced by muons 
far away of the detector (see figure 7.9). 

New Micromegas detectors of the microbulk type have been manufactured for the 
line. In contrast with previous designs the readout has been modified and the strips 
pattern has a smaller pitch (500 pm instead of 550 pm) by keeping the same active 
area (60 X 60 mm 2 ) and thus the number of strips has been increased to 120 per 
axis. This new design is the result of the studies done on the Micromegas detectors 
in order to enhance its performance and the improvements on the manufacturing 
technique. So far is the detector with the better performance working at CAST 
with a 13% of FWHM in the 5.9 keV peak (see figure 7.12 left). Also, it shows 
an excellent spatial resolution and homogeneity of the gain in the active area (see 
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Figure 7.10: Top: Step file of the novel Sunrise Micromegas design in which the different 
parts or the chamber are labeled. Bottom: Photo of the new Sunrise Micromegas detector 
glued on the copper raquette. The active area is on the left while the strips signals are 
extracted from the central connector on the right. The HV is fed by the lines situated on 
the right. 

figure 7.12 right). 

In the new design the quantum efficiency of the Micromegas chamber has been 
enhanced, because the expected region of the axion signal is minimized after being 
focused by the XRT. Consequently the cathode pattern has been modified, now 
it has a spider web design with a central hole of 8.5 mm of diameter, big enough 
to contain the expected axion signal. Indeed, the X-rays focused by the XRT go 
through the 4 pm aluminized polypropylene window avoiding the grid structure, 
which was responsible of a ~10% of efficiency loss in previous set-ups. The new 
strongback together with its projection during calibrations in the Micromegas are 
shown in figure 7.13. The expected focusing spot region for the X-rays has been 
measured using an Amptek COOL-X [133] X-ray generator, placed on the other 
extreme of the magnet (~ 14 m far away). The result of these measurements are 
shown in the right part of figure 7.13. 

Also the DAQ is equipped with the AFTER front-end electronics for the strips 
readout and a new acquisition software has been developed, based on C++ and 
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Figure 7.11: Scheme of the XRT+Micromegas line in which the different parts of the line 
have been labeled. 



Figure 7.12: Left: 55 Fe calibration spectrum of Sunrise Micromegas. Right: Gain unifor¬ 
mity in Sunrise Micromegas. The dead areas (in purple) show lower values than the unity 
(in green) and lie outside the axion-sensitive area. 


ROOT. Although the mesh pulse is acquired by the Matacq board as the previous 
acquisition, the electronic chain has been modified in order to implement the stand¬ 
alone electronics for the strips. 

After the implementation of the upgrades described before, the background level 
of the Sunrise Micromegas detector drops to (0.8 ± 0.2) x ICR 6 c cm~ 2 keV -1 s _1 , 
in this case only a ~25% of the background events are rejected by the muon veto. 
This issue could be due to the non optimum coverage of the plastic scintillator in the 
set-up. Because of geometrical constraints in the experimental area, only a plastic 
scintillator has been installed. The background spectra is shown in figure 7.14, it 
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StripsX StripsX 

Figure 7.13: Left: Photo of the new copper strongback with the spider web design. Cen¬ 
ter: Projection of the cathode in the Micromegas detector during the calibrations. Right: 
Spot generated during the X-ray finger run in the Micromegas detector, the blue circle 
represents the central hole in the cathode with a diameter of 8.5 mm. 



Figure 7.14: Background spectrum of the Sunrise Micromegas detector during the 2014 
data taking campaign. The blue line corresponds to the spectrum generated before the 
muons subtraction while the red line is the final spectrum after the veto cut. 


is dominated by the 8 keV copper fluorescence, which can be a hint of the effect 
of the muons in the set-up. Nevertheless, this is the best background level reached 
in the CAST experiment with a Micromegas detector, the first time below the 
10~ 6 c cm” 2 keV^ 1 s -1 level in stable conditions. 

7.5 Future prospects 

The new XRT + Micromegas line in the CAST experiment has set a milestone in 
axion research. For the first time an X-ray focusing device has been constructed 
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specifically designed for axion searches, in its focal plane a low background Mi- 
cromegas detector has been placed, showing an excellent background level and de¬ 
tector performance. 



Figure 7.15: Expected sensitivity of the new vacuum phase in the CAST experiment (red 
line), in contrast with the current CAST limit (blue line). ALP hints and theoretical limits 
introduced in chapter 1 are also drawn. 

CAST will finish the rescanned vacuum phase at the end of 2015 improving its pre¬ 
vious limit, lowering the coupling limit to an expected value of ~ 6 x 10 -11 GeV -1 . 
This number has been obtained using a realistic Monte Carlo simulation imple¬ 
menting the background levels of the detectors described before and the unbinned 
likelihood method presented in section 6.2. The results of the simulation, together 
with the current CAST limit are shown in figure 7.15. 

Although CAST will improve its sensitivity for axions by reducing the background 
level of the detectors, the discovery potential of CAST is limited by the size of 
the magnet. In order to scan a wider region in the m a — g ai parameter space the 
IAXO experiment has been proposed. The design and the potential of IAXO will 
be detailed in the following chapter, in which ultra-low background X-ray detectors 
are required. It will push the low background techniques to another stage opening 
unexplored R&D lines. 
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8.1 Introduction 

The CAST experiment is the most sensitive axion helioscope so far, although its 
main components have been "recycled" for axion physics. Beyond CAST, a new he¬ 
lioscope with an improved sensitivity, specifically built for axion and ALPs searches, 
has been proposed: IAXO-the International AXion Observatory. The Letter of In¬ 
tent [135] for IAXO has been submitted to CERN with a positive recommendation 
and the Conceptual Design Review [136] has been already published. IAXO will ex¬ 
ploit the helioscope technique with a dedicated magnet, optics and low background 
detectors, which will be detailed in this chapter. 

IAXO will have sensitivity to the axion-photon coupling of more than one order of 
magnitude beyond CAST, entering a large fraction of unexplored parameter space. 
IAXO could become a generic facility for axion research. The possibility of hosting 
dark matter axion detectors in IAXO is under study. The details of the physics 
potential of IAXO will be also described. 

8.2 The IAXO proposal 

IAXO will enhance the helioscope technique by exploiting all the singularities of 
CAST presented in previous chapters, implemented into a large superconducting 
toroidal magnet, together with X-ray optics and low background detectors attached 
at the end of the magnet bores. In this section these three main elements proposed 
for IAXO will be described. 
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8.2.1 The IAXO superconducting magnet 

The sensitivity of CAST is currently limited by the size of the magnet, a decom¬ 
missioned LHC dipole magnet recycled for axion searches. Indeed, the helioscope 
technique could be enhanced using a longer magnet with a bigger aperture, specif¬ 
ically designed for axion physics. A new toroidal magnet, inspired by the ATLAS 
design, has been proposed for IAXO. Equipped with eight magnet bores, with an 
aperture of 60 cm of diameter each and a length of 21 m. Taking advantage of the 
NbTi superconducting technology, that allows peak magnetic fields up to 5.4 T with 
an stored energy of 500 MJ at an operational current of 12.3 kA. 



Figure 8.1: Left: Mid plane cut of the cryostat of the proposed IAXO magnet, the different 
parts are labeled. Right: Magnetic field simulation of the IAXO magnet, the rectangles are 
the toroid’s coils while the circles represent the position of the optics. Plots taken from 
[135], 

An optimization study [137] shows that the sensitivity of IAXO is affected consid¬ 
erably by the fraction of the aperture of the telescopes exposed to X-rays. So it is 
preferable to use thinner coils, increasing the open aperture in front of the telescopes 
(see figure 8.1 right). However, in this configuration the average magnetic held drops 
to 2.5 T. On the other hand, the design features the decoupling of the magnet sys¬ 
tem from the optical detection systems, which simplifies the system integration (see 
figure 8.1 left). Also, it allows open bores that are centered and aligned in between 
the racetrack coils in accordance with the geometrical study. The inclusion of eight 
warm bores will simplify the use of experimental instrumentation and the periodic 
maintenance of the system. 

The coil windings will be cooled by conduction at a temperature of 4.5 K. The 
conceptual design of the cryogenic system is based on a forced how of sub-cooled 
liquid helium at supercritical pressure. The cold mass operating temperature is 
4.5 K and its mass is approximately 130 tons. It consists of eight coils with two 
double pancakes per coil, which form the toroid geometry and a central cylinder is 
designed to support the magnetic force load. 
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The IAXO detectors will be placed in a light and confined structure, such as a dome 
or a framed tent that will serve as the main site for the experiment. IAXO will need 
to track the Sun for the longest possible period in order to increase the data-taking 
efficiency. Thus, the magnet needs to be rotated both horizontally and vertically 
by the largest possible angles. A vertical inclination of ± 25° is required, while the 
horizontal rotation should be stretched to a full 360° that will allow IAXO to perform 
trackings of about 12 h. The 250 tons magnet system will be supported at the center 
of mass of the whole system at the cryostat central post (see figure 8.2). The vertical 
movement is performed by two semi-circular structures while the rotation of the disk 
is generated by a set of roller drives on a circular rail system. The required magnet 
services providing vacuum, helium supply, current and controls, are placed on top 
of the disk to couple their position to the horizontal rotation of the magnet. 



Cryostat 


Flexible Lines 


Inclination System 


"elescones 


Support Frame 


Rotating Disk 


Services 


Rotation System 


Figure 8.2: Schematic view of the IAXO proposal. The different parts of the magnet, 
magnet support and feeds are labeled. Plot taken from [135]. 


Although the design of the magnet is based on the experience gained on the ATLAS 
toroid, the IAXO magnet will deal with a peak magnetic field of 5.4 T, which is not 
trivial in terms of superconductor development and training behavior of the coil. 
In order to validate the design a single short prototype coil, named TO, has been 
proposed. The assembly of the TO prototype for IAXO is under study, waiting for 
the approval of CERN. 
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8.2.2 X-ray optics for IAXO 

As it was described before, the purpose of the X-ray optics is to focus the expected 
X-ray signal to a small spot in the detector, projecting the signal in a tiny region and 
increasing the signal-to-noise ratio. The performance of X-ray optics is characterized 
mainly by its efficiency, defined as the fraction of the X-rays focused by the optics 
and the size of the spot in the detector. 

In order to archive the smallest spot a the optics should have a short focal length 
/, since the spot area grows quadratically with the focal length. At the same time 
the individual mirrors that compose the optic should have the highest reflectivity to 
X-rays, which increases by decreasing the graze angle a. Since / oc a -1 , the optics 
should have a long enough focal length. Indeed, the complication of the optical 
design is that the efficiency and the size of the spot have a complex dependence on 
the incident energy E and grazing angle a. 

There are different manufacturing technologies of reflective optics, in the case of 
IAXO the segmented and slumped glass optic has been selected, it is a consolidated 
technology used recently for the NuSTAR satellite mission. Also, this approach 
facilitates the deposition of single-layer or multi-layer reflective coatings, being the 
least expensive of the fabrication techniques. Moreover, the requirements for the 
angular resolution of IAXO are gentle, although other optics technologies may have 
better resolution than slumped glass, they would not produce a significantly smaller 
spot of the solar core. 




Figure 8.3: Left: Focal length versus spot size y/a (red squares) and the throughput of 
the optics (blue dots) which is calculated dividing a quantity related with the efficiency of 
the optics (DAF) by yfa. The optimum is found for / = 5 m. Right: Isomorphic side view 
of the telescope in which the "spider-web" structure is shown. Plots taken from [135]. 

After a systematic study of the throughput of the optics (see figure 8.3 right), the 
optimal focal length for IAXO is / = 5 m with a focusing spot of about a — 0.2 cm 2 . 
A detailed calculation of the optimization process can be found in [138]. The IAXO 
optics will be composed by 123 nested layers with a W/B 4 C coating, covering the 
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total magnet aperture of 60 cm of diameter. All the eight magnet bores will be 
equipped by an X-ray optics that will be mounted on a support structure with a 
"spider-web" design, as it is shown in the left part of figure 8.3. 

For instance, the manufacturing process and the technology of the new X-ray optics 
installed at CAST during 2014 in the Sunrise side, presented in section 7.4.2, are 
the same that the one proposed for IAXO. Although it has a small size, it can be 
considered as the first prototype of a telescope for axion research. 

8.2.3 Ultra-low background X-ray detectors for IAXO 

The baseline technology of the low background X-ray detectors for IAXO are the 
Micromegas detectors described before. CAST microbulk detectors have achieved 
record levels in terms of background and they offer the best prospects to meet the 
requirements of IAXO. 
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Figure 8.4: Micromegas background history since 2002. The blue dots correspond to 
the unshielded Micromegas with the classic technology at CAST, the green squares are the 
background levels after the installation of the shielding, the black triangles correspond to 
the shielding upgrade from 2012 to 2014 and the red dots are the levels reached at the LSC. 
The dashed black line marks the goal for IAXO. 


Although the research in low background techniques on Micromegas detectors, pre- 
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sented in chapter 7, has led to a reduction of the background level in a factor ~ 6 
and more than two orders of magnitude since the beginning of CAST (see fig¬ 
ure 8.4). An ultra-low background detector is required for IAXO, with a goal 
of ~ 10~ 7 c cm -2 keV -1 s -1 and down to ~ 10“ 8 if possible. 

The state of art in low background techniques is summarized in the upgrade of the 
Sunrise detector during 2014 at CAST, detailed in section 7.4.2. The new system 
(Micromegas + XRT) represents the first time that both technologies proposed for 
the future IAXO are used in conjunction. It can be considered as a IAXO pathfinder, 
being an important milestone for the technical design phase of IAXO. However, the 
final background is around one order of magnitude above the levels required for 
IAXO. In this way, the design of a DO detector for IAXO has been proposed, it 
will fed by the low background techniques developed until now. Moreover, new 
improvements and research lines have been proposed: 

• Veto coverage: The active muon vetoes installed at CAST are not in its op¬ 
timum configuration due to geometrical constraints in the experimental area. 
These spatial limitations could be removed in IAXO, extending the surface 
area of the muon veto as much as possible. Then, a higher rejection of events 
related with muons is expected. 

• New thin windows: The efficiency of the Micromegas detectors is limited (at 
least at low energies) by the cathode window. Thus the sensitivity could be 
improved by the use of a thinner window with a higher transmission of X-rays. 
Different materials are being investigated for IAXO. 

• New gas mixtures: The background level of the Micromegas measured under¬ 
ground could be limited by the 39 Ar, a natural isotope of the Ar. Thus new 
gas mixtures, like Xe -f- TMA 1 [139], have been proposed. Moreover these gas 
mixtures could improve the quantum efficiency of the detectors. 

• AGET front-end electronics: Although the electronics of the Micromegas de¬ 
tectors were upgraded during 2013 to the new AFTER chip. The novel 
AGET[134] electronics keep its main features, complemented by an auto¬ 
trigger functionality for every single strip. So the low energy threshold of 
the detectors could be increased. It will open the sensibility of IAXO to new 
physics that will be presented in the following section. 

• Resistive Micromegas: Resistive Micromegas may fix the main problem of 
this kind of detectors. Indeed, the intrinsic gain of the regular Micromegas is 
limited by the occurrence of sparks in the detector, which can be a destructive 
process. The use of a Micromegas with a resistive layer between the mesh and 
the anode readout [140] allows to work at higher gains, which can be optimized 
for axion searches. 

1 Trimethylamine 
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These new R&D lines will establish the roadmap on the reduction of the background 
level of the Micromegas detectors. 

Although the features described before complete the baseline of IAXO, the instal¬ 
lation of additional equipment is under study. They offer potential to span the 
detection energy window for solar ALPs or WISP to lower energy ranges: GridPix 
detectors, Transition Edge Sensors (TES) and low-noise Charge Coupled Devices 
(CCD), have been proposed. Moreover, the sensitivity of IAXO to relic dark matter 
axions and ALPs could be performed by the use of microwave cavities or antennas, 
that will be described in section 8.3.2. 

8.3 Physics potential 

In this section the physics potential of IAXO will be described. The sensitivity to 
hadronic solar axions and ALPs, which is the baseline for IAXO, will be calculated. 
Also, the sensitivity of IAXO to non-hadronic solar axions will be evaluated. Finally, 
the potential of IAXO for the detection of relic dark matter axions and ALPs will 
be described. 

8.3.1 Expected sensitivity for solar axions and ALPs 

The expected sensitivity of IAXO to solar axions coming from the axion-photon 
conversion has been computed. For this purpose a complete Monte Carlo simulation 
has been developed, in which the expected background counts in the spot area have 
been taken into account. A coupling limit has been derived, assuming no signal, 
using the unbinned likelihood method presented in section 6.2. Two scenarios have 
been proposed and implemented in the simulations, one "nominal scenario" and 
another "enhanced scenario", their main parameters are shown in table 8.1. 


Parameter 

Units 

IAXO Nominal 

IAXO Enhanced 

Average magnetic field 

T 

2.5 

2.5 

Magnet length 

m 

20 

20 

Magnet bore area 

m 2 

2.3 

2.3 

Background level 

10 -8 c 

keVcm 2 s 

5 

1 

Detector efficiency 


0.7 

0.8 

Optics efficiency 


0.5 

0.7 

Spot area 

9 

cm 

8x0.2 

8x0.15 


Table 8.1: Values of the IAXO parameters taken into account in the simulations. Two 
different scenarios are displayed, the nominal and the enhanced one. 

Regarding exposure, two different phases have been proposed (see table 8.2): IAXO 
Run-I with vacuum in the magnet bores and 3 years of effective data taking (4 years 
of total duration for a 75% assumed duty cycle) will determine the sensitivity of 
IAXO for axion masses below m a < 0.01 eV and IAXO Run-II using 4 He as buffer 
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gas, that will recover the coherence for higher axion masses. The current sensitivity 
curves are calculated assuming that the gas density in Run II is continuously changed 
from 0 to 1 bar of 4 He at room temperature during a total effective data taking time 
of 3 additional years. 


Run: 

Run I 

Run II 

Total duration 

~4 years 

~4 years 

Effective data taking time 

3 years 

3 years 

Effective exposure 

9540 h 

9540 h 

Gas pressure 

0 bar 

0-1 bar 

Axion mass 

m a < 0.01 eV 

0.01 < m a < 0.25 eV 


Table 8.2: Values of the exposure for the different phases of IAXO used to compute the 
sensitivity. 
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Figure 8.5: Expected sensitivity of IAXO for hadronic axions, compared with current 
bounds from CAST and ADMX. Also future prospects of ADMX (dashed brown region) 
and ALPS-II (light blue line) are shown. For the sake of clarity the labels from other bounds 
or regions, which are shown in figure 1.4, have been removed. Plot taken from [135]. 
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The obtained values for the respective scenarios are represented by the couple of lines 
bounding the dashed area in figure 8.5. As shown, IAXO will be more than one order 
of magnitude sensitive than CAST in terms of the axion-photon coupling constant. 
IAXO could be sensitive to coupling constants of about g ai ~ 5 x 10- 12 GeV 1 for 
axion masses up to 10 meV and around g a7 10 11 GeV 1 up to 0.25 eV. 

IAXO will deeply enter into a completely unexplored ALPs and axion parameter 
space excluding a large region of the QCD axion phase space that has yet to be 
explored. Also, IAXO will explore a favored parameter space region for axions 
and ALPs, given by the white dwarf cooling hint and the VHE transparency hint, 
described in section 1.5. 


8.3.1.1 Axion-electron coupling 

IAXO could be also sensitive to non hadronic solar axions, introduced in sec¬ 
tion 1.7.1.2. These axions could be generated in the solar core via axion- 
Bremsstrahlung, Compton and axio-deexcitation processes and the expected flux 
at Earth could be considerably larger than the Primakoff emission. However, the 
differential spectrum is shifted to lower energies ~ 1 keV (see figure 1.10). An he¬ 
lioscope could be sensitive to this kind of axions assuming a non-hadronic axion 
emission in the Sun and the inverse Primakoff conversion inside the magnet. Nev¬ 
ertheless, in this case the expected signal depends on g ae x g a7 , the product of the 
axion electron coupling constant and the axion-photon coupling. 

The plot on the left of figure 8.6 shows the computed sensitivity of IAXO to the 
product gaeda-y , assuming that the Primakoff emission from the Sun is negligible in 
comparison with the non-hadronic. The computation is performed in a similar man¬ 
ner and with the same assumed parameters than in the previous section. However, 
in this case the energy threshold for the detectors is set at 0.5 keV, with background 
levels and efficiencies comparable to the ones in the previous section down to this 
threshold. IAXO could be able to constrain g ae x g ai < 2.5 x 10~ 25 GeV -1 at a 95% 
of C.L. for axion masses m a < 0.01 eV. 

The Primakoff emission could be also included in the calculations, in this case the 
signal depends on three parameters: g ae , g ai and m a . For the vacuum case, in the 
coherence region, for m a < 0.01 eV, the detection is independent of the axion mass 
and a limit can be computed in the g ae — g a7 parameter space, as it is shown in the 
right part of figure 8.6. 

IAXO could directly measure solar axions produced by non-hadronic processes, for 
the first time with sensitivity to relevant g ae values. However, in this case the 
reduction of the low energy threshold and the intrinsic efficiency of the detectors at 
low energies will be crucial. 
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Figure 8.6: Left: IAXO sensitivity line for g ae x g ar/ as a function of m a , assuming a 
non-hadronic solar axion emission. The orange band is related with the DSFZ model. The 
yellow band corresponds to models in which the relation of g ae with m a is more suitable. 
The recent limit of CAST [75] is also shown. Right: IAXO sensitivity on g ae and g ai 
for m a < 0.01 eV. The gray region is excluded by solar neutrino measurements and the 
orange band corresponds to values of g ae ~ 1 — 5 x 10 " l:! . The yellow band represents 
the more plausible models. The recent limit of CAST [75] is also shown. Helium-burning 
stars would be perceptibly affected in the blue band, in which the parameters above are 
excluded. Plots taken from [135]. 


8.3.2 Search for relic CDM axions in IAXO 

As it was presented in chapter 1 axions and ALPs are attractive dark matter can¬ 
didates and could take account of all the amount of the CDM in the Universe. 
The huge magnet required for IAXO offers excellent possibilities to host relic DM 
searches that can be performed using two different techniques: haloscopes and dish 
antennas. 

The haloscope technique was introduced in section 1.6.1.1. Relic axions or ALPs 
could be converted into photons inside strong magnetic fields and detected by the 
use of microwave cavities. The features of the haloscope technique, together with a 
directional effect on CDM axion searches will be described in chapter 9. 

A new concept for axion DM detection consists in a spherical reflecting dish (em¬ 
bedded in a magnetic field) antenna, which reacts to DM axion particles emitting 
radiation focused on its center, where the detector lies [141]. This technique, con¬ 
trary to haloscopes, does not require a tuning of the experiment to the unknown m a . 
The accessible axion mass range is in practice determined by the detector sensitivity. 
The power received in the detector due to axion DM is given by: 
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Pout ~ 9 -^ PDM AB 2 g d ~ 0.5 xC'W (8.1) 

where A is the area of the dish and Q d a geometrical factor with a value of about ~ 1 
and pdm is the total amount of dark matter density in the Universe. The dish search 
does not rely on resonant enhancement available for a cavity search, but it is com¬ 
pensated if a large area for the dish is available. This technique compares favorably 
to resonant cavities for relatively larger axion masses m a > 1 meV, but for these 
large values is not sensitive enough to reach the QCD axion band. Nevertheless, is 
still very interesting for ALPs in certain cosmological scenarios. Dark matter exper¬ 
iments in IAXO would benefit of ultra-low temperatures, reducing the thermal noise 
as much as possible. This suggests to host the experiments in the cold part of the 
magnet. One possibility is to design one of the IAXO bores to remain at liquid-He 
temperature to benefit directly from the cryogenics. 
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9.1 Introduction 

Although helioscopes are one of the most promising techniques for axion discovery, 
axions might be detected using the haloscope technique, in which relic CDM axions 
could be converted into photons in a strong magnetic field via inverse Primakoff 
effect. These photons could be detected in resonant cavities if the frequency of the 
cavity has been properly tuned. The main features of the haloscope technique will 
be described in this chapter. 

On the other hand, a directional effect could be observed in the haloscope technique. 
This effect has been studied in [142], where the dimensions of the magnet and the 
cavities have an important role. These features will be detailed below. 

9.2 The haloscope technique 

Assuming that axions are the dominant component of dark matter, relic axions 
could be directly detected using the haloscope technique. Contrary to solar axions, 
relic axions are non-relativistic, and thus the energies of the resulting photons from 
the conversion are about the corresponding axion mass, in the microwave regime. 
If the conversion happens in a microwave cavity that is resonant to the axion mass, 
the conversion is substantially enhanced and a high sensitivity can be obtained to 
explore realistic QCD axion models. This technique has been already used in a large 
number of experiments, being ADMX the most powerful haloscope until now [36, 37], 
with a sensitivity to QCD axions in the /j>eV range. 
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The haloscope technique is based on a tunable microwave-cavity in a strong magnetic 
field, coupled to an ultra-low-noise microwave sensor. CDM axions and ALPs may 
convert into photons in the microwave regime with an enhanced probability if the 
resonant frequency of the cavity matches the axion energy E = m a + K , where K 
takes account of the velocity distribution of the CDM axions in the halo. Since 
m a is unknown, different axion masses have to be smoothly scanned, by tuning the 
cavity resonances. The power output of the cavity is given by [33]: 

Pout = 9 a 1 VB 2 C — Q (9.1) 

1 m a 

where V is the volume of the cavity, B the magnetic field and Q the quality factor of 
the cavity. p a is the dark matter density of the axions in the halo and C a geometry 
factor involving the precise electric field of relevant resonant modes in the cavity 
E cav {x ) and the magnetic field B(x): 

[f dVE cav (x)B(x)) 2 

“ V\Bl?fdVe(x)El v (x) ' ' 

here e(x) is the dielectric constant of the cavity. The previous equations are valid 
under the basic assumption that the de Broglie wavelength of the relic axions X a 
is longer than the characteristic size of the cavity d\ \ a > d. In this case the 
resonant conversion takes place. CDM axions inside the galactic halo have velocity 
distributions given by the velocity of the Sun referred to the galactic center v ~ 
300 km/s. So approximately, the de Broglie wavelength of the relic axions is given 
by: 


A 


a 


— ~ 12.4m (ITAD f 300km/s\ 
Pa \ m a ) \ Va ) 


(9.3) 


From equation 9.3, the condition A a > d is fulfilled for axion masses below 10 -4 eV, 
for a magnet geometry in the ~m scale. For instance, the ADMX experiment 
employs a cylindrical cavity of 1 nr length and 0.6 m diameter inside a solenoidal 
magnet. The cavity is tunable to axion masses at the few /reV scale, with a sensitivity 
sufficient to exclude the KSVZ model. For higher axion masses the main challenge 
relies on the fact that smaller resonant cavity geometries are needed, which could 
worsen the corresponding sensitivity. Recently the use of long and thin cavities 
(waveguides) inside strong dipole magnets has been proposed as a possibility to 
achieve competitive sensitivity in the 10 -5 — 10 -4 eV range [143, 144], The use of 
small cross-section, but long and powerful magnets like the ones used in particle 
accelerators, can accommodate cavities resonant at these higher frequencies (driven 
by the small dimension of the waveguide) while keeping a large enough volume and 
magnetic held. 

The possibility of using few-meter long cavities for detection of 10 5 eV relic 
axions is closer to the limitation given by the de Broglie wavelength. In this chapter 
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the effect of this limitation on the predicted relic axion signal for different dark 
matter distributions is explored. Particularly for thin and long geometries, the 
orientation of the cavity with respect the main incoming axion direction may affect 
the signal intensity. This effect can be maximized and used as an identification 
signature of the origin of an eventual positive detection and will be detailed in the 
following section. 


9.3 Directional sensitivity of CDM Axions. 

As it was presented before, expression 9.1 is calculated assuming that de Broglie 
wavelength of the incoming axion is larger than the dimensions of the cavity A a > d. 
Thus, the axion is approximated by a spatially constant oscillating held. If this 
condition is relaxed, the spatial variation of the axion held along the cavity volume 
cannot be neglected and thus equation 9.1 has to be modified: 

Pout = ga 1 VB 2 C—Q (9.4) 

1 m a 



Figure 9.1: Left: Sketch of the cavity geometry used for the calculations. Right: \P\ 2 
versus k z L for a single incoming axion momentum k as expressed in equation 9.5. Plots 
taken from [142]. 

The geometry factor C includes now a more complex dependency with the axion 
momentum and the cavity held. In the case of the use long thin rectangular resonant 
cavities along the z direction and inside a dipole magnet with a constant magnetic 
held (see figure 9.1 left), the more relevant cavity mode is the TEioi, the fundamental 
one with the E held parallel to B. This approximation assumes that the typical 
axion wavelengths may be comparable or shorter than the length L of the magnet, 
but otherwise the de Broglie condition is preserved. However, in the case of a 
single incoming axion direction with a momentum k inside a thin cavity, C can be 
computed analytically, by the expression: 
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C = C\F\ 2 


(1 + cos (k z L))ir 4 

2(tt 2 -A;2)2 


T = 


7r 

2L 



(9.5) 


here T is a form factor expressing the loss of coherence due to the axion momentum 
along the length of the cavity and k z is the projection of the axion momentum along 
the z direction of the cavity. 

The form factor \F\ 2 versus k z L is plotted in figure 9.1, in which the anticipated 
behavior is clearly seen. It provides the conventional result C = C for low values 
of k z L < 1, while the signal drops for larger values of k z L. Note that low values 
of k z L are achieved by small L or k, but also for axion directions perpendicular to 
the cavity length. This suggests that full coherence is possible even for long thin 
cavities when they are oriented perpendicularly to the axion direction. However, a 
more realistic treatment using a distribution of the axion incoming directions has 
to be done and will be studied in the following section. 


9.3.1 Isothermal sphere halo model 

The velocity distribution of dark matter axions at the Earth depends on the assump¬ 
tions considered for the halo model. So it has to be consistent with the observed 
rotation curve of the galaxy, keeping this main constraint, a large number of dif¬ 
ferent halo models can be considered. The velocity distribution of the dark matter 
particles at Earth has been studied in the context of WIMP 1 dark matter experi¬ 
ments. For instance, in reference [145] different halo models and their corresponding 
velocity distributions are described. The simplest one is the isothermal sphere halo 
model in which the distribution of velocities of the CDM particles in the halo follows 
a Maxwellian: 


( " 3 “ 2 ^ 

f(v ) = f{y) oc e\ 2v ™ a J (9.6) 

here v = |u| is the normalized velocity and v rms is the root mean square velocity 
v rms = -\/3/2uo, being vq the rotation speed of the galaxy at the solar system radius. 
In addition to the particular shape of f(v), the effect of the movement of the Earth- 
Sun system through the galactic dark matter halo has to be taken into account. 
In the terrestrial frame of reference, the velocity distribution function is given by 
f(w), being w the incoming velocity of the axions at Earth. f(w) is derived from 
f(v), using the relationship v = w + u where u is the velocity of the Earth in the 
galactic rest frame. This translation produces a general anisotropy in the velocity 
distribution at the Earth frame of reference, centered in a certain point in the sky, 
also referred as CYGNUS. The velocity u is of the order 220 km/s with an oscillatory 
component of about 12 km/s due to the rotation of the Earth around the Sun. 

1 Weakly Interacting Massive Particles 
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Figure 9.2: Left: Mollweide projection of the velocity distribution predicted by the 
isothermal sphere halo model. Right: Value of |J-/ 5 | 2 for all possible 2-angle orientation 
of the cavity using the Mollweide projection in galactic coordinates, for m a = 6 x 10 -5 eV 
and L = 15 m. The signal goes to the minimum when the cavity is oriented to the 
CYGNUS point and to the maximum when is perpendicular to it. Plots taken from [142]. 


The distribution of f(w) corresponding to the isothermal sphere model is shown on 
the left part of figure 9.2, in which the Mollweide projection in galactic coordinates 
has been used. As it is shown, the distribution is dominated by the anisotropy 
introduced by the Earth-Sun motion. The constant C has been calculated in this 
distribution, by convoluting the form factor from 9.5 with the momentum distribu¬ 
tion provided by the isothermal sphere model, given by the expression: 

= [ f(k)\T(k)\ 2 dk (9.7) 

Jk 




Figure 9.3: \Fis\ 2 versus magnet orientation with respect the CYGNUS point (solid blue 
line) compared with the idealized situation of a single-direction incoming axion (red dashed 
lines, corresponding from top to bottom to v a = 200, 300 and 400 km/s). The left plot 
corresponds to m a = 4 x 10~ 5 eV and L = 20 m and the right one to m a = 8 x 10~ 5 eV 
and L = 15 m. Plots taken from [142]. 

The factor \J-is\ 2 versus the cavity orientation angle is shown in figure 9.3, in 
which the orientation angle is defined with respect the CYGNUS point in the sky. 
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The curves are compared with idealized curves assuming a monochromatic single¬ 
direction axion beam using expression 9.5. As expected, the momentum dispersion 
of the distribution causes a smoothing of the modulation of the signal, in contrast 
with the idealized case. However, a significant modulation still remains for some 
values of m a and L. In addition, the signal at the maximum decreases for longer 
cavities with respect to full coherence (|J r /s , | 2 ) because the momentum dispersion 
of the distribution prevents to achieve the condition k z L = 0 of the idealized case. 
On the right part of figure 9.2, \J~is\ 2 is plotted for all possible 2-angle orientation 
of the cavity in the particular case of m a = 6 x 10' 5 eV and L = 15 m. 

Moreover the factor |-7 r /s , | 2 has been calculated in a more systematic way, using 
different values of m a and L. On the left part of figure 9.4 the maximum value of 
\J~is \ 2 is plotted, it is the case when the cavity is oriented in perpendicular to the 
CYGNUS point in the sky. As shown, for cavities shorter than ~ 5 X (10~ 4 eV/m ) 
meters the full coherence is preserved. For larger cavities the intensity decreases 
significantly. 



Figure 9.4: Left: Maximum value of | Fis\ 2 versus m a L with the cavity oriented perpen¬ 
dicular to the CYGNUS point. Right: Difference of the maximum and minimum values of 
\F IS \ 2 (with the cavity oriented perpendicular and parallel to the CYGNUS point respec¬ 
tively) versus m a L. Plots taken from [142]. 


On the right part of figure 9.4 the difference of \J 7 is\ 2 at the maximum and the 
minimum of the modulation is plotted, with the cavity oriented perpendicular and 
parallel to the CYGNUS point respectively. As shown, for an adequately chosen 
length of the cavity, a modulation difference as large as 35% of the cavity power is 
expected. This happens for lengths of ~ 10 x (10 _4 eV/m a ) meters, which correspond 
to geometries in which the expected signal is a 60-70% of the full coherence. This re¬ 
gion should be noted as a desirable operating point, providing a strong identificative 
signature of the direction of the CDM axions. 
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9.3.2 Sensitivity to CDM axions 

In this section the experimental issues concerning the axion model parameters which 
can give a detectable signal in the cavities described above will be treated. An 
estimation of the sensitivity on the axion-photon coupling g a ~ has been performed, 
for a range of axion masses from 10~ 5 < m a < 1.5 x 10~ 4 eV. For some reasonable 
experimental parameters and cavity geometries consistent with the requirements 
presented in previous sections. 

The computation has been performed in the usual way for haloscopes [146], but 
using the signal strength given by equation 9.1, with a cavity orientation giving 
the maximum |J r / 5 , | 2 and cavity geometries fixed to the corresponding axion mass. 
Under the assumption that all the CDM is composed by axions with a density of 
p a = 0.3 GeV cm' 3 and using the velocity distribution of the isothermal sphere 
model described above. 

The values of the different parameters concerning the sensitivity calculation are fixed 
after the following considerations: 

• For each axion mass m a , the transversal dimension d y of the cavity is fixed by 
the resonance condition. The remaining dimension d x is fixed at 3 cm for all 
the calculations. 

• The length L is fixed differently in diverse ranges of m a , so it lies approximately 
in the region indicated in the right plot of figure 9.4, having the maximum 
modulation. 

• The specific values taken for L for several values of m a are listed in table 9.1. 


Axion mass (eV) 

Length (m) 

| -pM AX 12 
K/S 1 

| pMAX 12 \-pMlN\2 

K IS 1 rIS 1 

Vol (1) 

2.0 x 10 -5 

20 

0.95 

0.13 

18.59 

4.0 x 10" 5 

20 

0.85 

0.28 

9.30 

6.0 x 10 -5 

15 

0.83 

0.30 

4.65 

8.0 x 10 -5 

15 

0.75 

0.35 

3.49 

1.0 x 10“ 4 

10 

0.80 

0.32 

1.86 

1.2 x 10" 4 

10 

0.75 

0.35 

1.55 

1.4 x 10 -4 

10 

0.70 

0.37 

1.33 


Table 9.1: Values of the length L used for the calculations on the sensitivity for several 
axion masses, as well as the maximum value of the factor [Tis] 2 , its modulation (maxi¬ 
mum minus minimum), and the cavity volume, corresponding to each case. Values taken 
from [142], 


However this procedure is an approximation, in the practice the cavities have to be 
built with a system to tune the resonance frequency by means of movable dielec¬ 
tric pieces inside the cavities. Moreover, the presence of dielectric pieces inside the 
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cavities must be taken into account in the calculation of C, due to the geometrical 
variation of the TE 101 mode and the presence of a dielectric material in the cav¬ 
ity. These issues have been neglected in the estimation of the sensitivity and the 
presented result can be considered within a factor ~ 2 of uncertainty. 

Following the calculations from [143], the quality factor Q of the cavity has been 
assumed to be constant and equal to 1000 with a thermal noise of 3 K. For every mass 
step, an integration time of 30000 s is assumed, this leads to a total of 2.5 effective 
years to scan the presented mass range. All these parameters are considered a priory 
feasible, although a detailed technical study is needed. Finally, a magnetic field of 
10 T has been assumed in the calculations. 



Figure 9.5: Sensitivity regions for the microwave cavities presented, with a geometry 
fixed by the desired directional effect. The conservative option (black region) is calculated 
for Q = 1000 and t s tep = 30000 s, using the parameters from table 9.1. For the red 
region more optimistic parameters have been considered with a quality factor of Q = 
20000, t s t e p = 3000 s and a volume 10 times larger. Some excluded regions by haloscopes 
and the CAST limit are also shown. The yellow band indicates the region favored by the 
QCD axion models, being the green lines two representative benchmark models, the KSVZ 
(upper one) and the DSFZ model. Plot taken from [142]. 


The upper line from figure 9.5 represents the sensitivity calculations on g ai giving a 
signal-to-noise ratio of 5, for the input parameters mentioned above. The expected 









9.3. Directional sensitivity of CDM Axions. 


155 


sensitivity in g ai is already at the level of the realistic QCD axion models. On the 
other hand, a more optimistic scenario has been computed for the lower line, in 
which a quality factor Q of 20000 has been assumed, as well as 3000 s of integration 
time (corresponding to 5 years of effective data taking time). Also a xlO larger 
volume V with respect the conservative case is considered. These improvements 
could be achieved by the use of 10 cavities of the dimensions stated with their 
power combined coherently. This second line deeply enters into the QCD axion 
band and particularly it covers both the KSVZ and DSFZ models in a wide range 
of the parameter space. 

The case studied involving thin cavities inside long magnets is particularly appealing 
because it could be realized in the near future, given that this type of magnets are 
already used by the axion community in experiments looking for solar axions, like 
CAST. Moreover the huge magnet required for IAXO offers excellent possibilities 
and the hosting of microwave cavities for axionic DM searches is under study. 





Chapter 10 


Summary and conclusions 


Axions are well motivated pseudoscalar particles proposed in an extension of the SM 
as a solution to the strong CP problem. Indeed, the theory predicts a CP-symmetry 
violation in the strong interactions that has not been observed experimentally. The 
most compelling solution to the strong CP problem was proposed by Peccei and 
Quinn in 1977, introducing a new global chiral symmetry U(1)pq that is sponta¬ 
neously broken at the energy scale of the symmetry f a . It solves the strong CP 
problem dynamically and a new particle appears as the pseudo Nambu-Goldstone 
boson of the new symmetry, the axiom 

The Peccei Quinn solution fixes some properties of the axions, like it mass and the 
coupling constant, related with the energy scale of the symmetry f a . Axions could 
interact with gluons, photons and fermions. The most interesting case for axion 
searches is its coupling with photons, which is generic to all the models. 

Even though axions are the best motivated particles proposed in the theory, there is 
also the category of Axion Like Particles (ALPs) or more generically WISP (Weakly 
Interacting Slim Particles). They share the same phenomenology of the axion, being 
light particles which couple to two photons. ALPs arise from extensions of the SM 
in which a new symmetry is broken at a high energy scale and also appear in string 
theory as the axion does. In contrast with axions, the coupling constant of ALPs 
is not related with its mass and ALPs might lie in a large region of the parameter 
space. 

Axions and ALPs could have been produced in an early Universe by non thermal 
mechanism like the vacuum realignment and the decay of the topological defects. 
Being electrical neutral particles which interacts weakly with the matter, axions and 
ALPs are attractive Dark Matter candidates, that could explain separately all the 
amount of DM in the Universe. Axions and ALPs properties are constrained by 
astrophysical and cosmological considerations and would play an important role in 
the stellar evolution. On the other hand, different experimental observations could 
be interpreted has a hint of axions and ALPs, such as the excessive transparency of 
the Universe to VHE photons and the anomalous cooling rate of the WD. 

Different techniques have been developed for axion searches: helioscopes looking 
for solar axions; haloscopes that search relic CDM axions and photon regeneration 
experiments, in which axions could be generated and detected in the laboratory. All 
these searching strategies are based on the Primakoff effect in which axions could 
be converted into photons inside strong electromagnetic fields. 
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Figure 10.1: Sketch of the detection of solar axions in the helioscope technique. 


The helioscope technique, which is the main topic of this work, was proposed by 
Sikivie in 1983 and uses the Sun as a powerful axion source. Solar axions could be 
generated in the core of the Sun in the strong electric held of the charged particles 
inside the plasma, via inverse Primakoff effect. Also, further processes like axion- 
Bremsstrahlung, Compton and axio-deexcitation, could be considerably larger than 
the Primakoff emission. These axions could be reconverted into photons inside 
strong magnetic fields via Primakoff effect (see figure 10.1). The related photons, 
which are in the X-ray regime (see figure 10.2 left), could be detected in the X-ray 
detectors placed at the magnet bore ends. 




Figure 10.2: Left: Differential solar axion flux at Earth for hadronic (black line) and 
no-hadronic (red line) axion models. Right: Axion to photon probability conversion in 
function of the axion mass for the vacuum case (black line) and the buffer gas case (red 
line). 

Although the non-hadronic axion emission could be larger in the solar core, only 
the Primakoff processes have been taken into account in this work because is more 
reliable to suppose the same processes involving the generation and detection of 
axions. Also, astrophysical constrains on g ae are more restrictive than the limits 
that could be in principle reached with helioscopes. The solar axion flux is well 
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established due to the detailed description of the Solar Model. The probability of 
the axion to photon conversion inside a strong magnetic held was introduced by van 
Bibber. In the helioscope technique the sensitivity to the axion mass is given by 
the conversion probability, which is constrained to small axion masses if the magnet 
bores are under vacuum (see figure 10.2 right). However, the coherence can be 
restored to higher axion masses using a buffer gas inside the magnet bores. 

The most representative experiment in the helioscope technique is the CAST (CERN 
Solar Axion Telescope) experiment, looking for solar axions since 2003 and being 
the most sensitive helioscope so far. CAST makes use of a decommissioned LHC 
dipole magnet with a length of 9.26 m and a magnetic held up to 9 T. The magnet 
is mounted on a movable platform which allows to tracking the Sun ~ 1.5 h two 
times per day, during sunrise and sunset (see figure 10.3). The magnet is composed 
by two bores with a total of four detectors placed at the magnet bore ends. The 
axion signal would be an excess of X-rays while the magnet is pointing the Sun and 
low background X-ray detectors are mandatory. 



Sunset 

photon detectors 


"Sunrise axions 


Low-background 

shielding 


Figure 10.3: A schematic drawing of the CAST experiment in which the different parts of 
the magnet and the movable platform are labeled. 


The CAST experiment is composed of several additional systems required for the 
data taking. The magnet operates at a nominal temperature of 1.8 K and a cryo¬ 
genic cooling system is required. Also, a vacuum system is installed around the 
magnet bores in order to increase the transparency for the X-rays from the axion 
conversion and to isolate the magnet from the environment. The movement of the 
magnet is controlled by a tracking system program which allows to point the Sun 
automatically. Also, a gas system is installed with the purpose of filling the magnet 
bores with a buffer gas in small steps. Finally, the monitoring of all the systems is 
performed by the slow control. 
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The CAST research program can be divided into two phases, the first one with 
vacuum inside the magnetic bores and a second phase with buffer gas. The Phase I 
data taking period was performed during 2003 and 2004, reaching a limit on the 
coupling constant of g ai < 8.8 x 10 -11 GeV~ 4 at a 95% of C.L., for axion masses 
m a < 0.02 eV. 

During the Phase II, from 2005 to 2006, the magnet was filled with 4 He, providing 
an experimental limit of g ai < 2.17 x 10~ 10 GeV -1 at a 95% of C.L. for axion 
masses 0.02 < m a < 0.39 eV. The 3 He phase started in 2008 and ended in 2011, 
scanning axion masses of 0.39 < m a < 1.17 eV and obtaining an experimental limit 
of g ai < 2.3x 10“ 10 GeV -1 for 0.39 < m a < 0.64 eV and g ai < 3.3xlO~ 10 GeV -1 for 
0.64 < m a < 1.17 eV at a 95% of C.L. CAST crossed first time the KSVZ benchmark 
model, one of the most favored in the theory. Although CAST finished its original 
research program in 2011, the data taking period has been extended. During 2012 
the 4 He phase was revisited, improving the previous limit in a narrow mass range. 
In 2013 CAST started a new data taking campaign revisiting the vacuum phase, 
motivated by the improvement of the background levels of the X-rays detectors and 
an improvement of the sensitivity is expected. 

Three different kind of X-rays detectors have been working at CAST since the 
beginning of the experiment: a Charge Coupled Device (CCD) on the focal plane 
of an X-ray telescope, a Time Projection Chamber (TPC) covering two magnet 
bores on the Sunset side and different types of Micromegas (MICRO MESh GAseous 
Structure) detectors that have been taking data in the Sunrise side and more recently 
in the Sunset side. 

The Micromegas detectors were developed by Giomataris in 1996. Micromegas 
are gaseous ionization detectors with two different regions separated by a mesh: 
the conversion region and the amplification region. In the conversion region the 
interacting particles ionize the gas generating primary charges, in this region an 
electric field is applied and the electrons are drifted to the amplification region. 
Here the avalanche of the primary electrons occurs, due to the strong electric field 
present in the gap (of ~ 50 /mi), generating two readable signals: one in the mesh 
that provides time resolution and another in the anode readout which confers spatial 
resolution to the detectors. A scheme of the working principle of the Micromegas 
detectors is shown in figure 10.4. 

Micromegas is a demanding technology which has experienced a technological evo¬ 
lution in the manufacturing process, from the classical Micromegas to the bulk and 
microbulk technologies. The CAST experiment has been a demanding test bench 
for these different manufacturing techniques. The latest microbulk technology is the 
one which shows the better performance and the lower intrinsic radiopurity. Three 
of the four detectors currently installed at CAST are of this type. 

The design of the microbulk Micromegas detectors installed at CAST during 
2011 has been carefully selected for axion searches. The detector anode is made 
of interconnected square pads which leads to a 2-dimensional strip readout of 
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Figure 10.4: Scheme of the working principle of a Micromegas detector. Two different 
regions are separated by a mesh: the conversion region and the amplification region. 


106 x 106 strips with a pitch of 550 /jrn and an effective area of about 60 x 60 mm 2 . 
The cathode is made of aluminized mylar of about 5 /jm thick, glued to a circular 
aluminum strongback, which is attached to the vacuum side of magnet. The body of 
the detector is made of Plexiglas, the chamber and thus the conversion volume has 
30 mm height. The Micromegas is glued to a Plexiglas base called raquette which 
has a circular shape in the detector area and a thin neck for the strips connections 
to the electronics. The gas in the chamber is an Ar + iC4Hi0 mixture, working at a 
pressure of 1.4 bar. The detectors are placed at the magnet bore ends covering the 
entire aperture of 14.52 cm 2 . Although Sunrise and Sunset detectors have different 
shielding designs, both of them are mainly composed by 5 mm of copper in the 
innermost part, 25 mm of archaeological lead and 2 mm of cadmium foil at the end, 
covered by polyethylene layers. 

The acquisition of the different readouts of the Micromegas detectors at CAST is 
performed by a Labview based program. The mesh pulse is digitized by a Matacq 
board and the anode readout is acquired by the Gassiplex font end electronics. The 
data analysis is done by a dedicated software based on C++ and ROOT. For the 
mesh signal a pulse shape analysis is made, this leads to the definition of different 
observables. For the strips readout a cluster analysis is performed, that allows the 
definition of different parameters related with the shape of the event. 

In a second stage, the discrimination of the X-ray like events is performed. For this 
purpose, the distribution of the background events is compared with the distribution 
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of the 55 Fe calibration events (see figure 10.5 left). The discrimination method 
has been developed computing the log-odds distribution of different observables for 
calibration and background events. It allows to define a cut value in which a certain 
number of events are accepted (see figure 10.5 right), that leads to the definition of 
the software efficiency. 



Figure 10.5: Left: 55 Fe spectra in a Micromegas detector in which the three different 
peaks are fitted. Right: Distribution of the log-likelihood for the 5.9 keV peak (blue), 
3 keV peak (green) and background events (red). The blue filled area corresponds to the 
cut-value at 5.9 keV while the green filled area is the the cut at 3 keV. 

Using this discrimination method, the data of the Micromegas detectors at CAST 
during 2011 have been analyzed. The discrimination capabilities of the detectors 
are strongly dependent on the detector performance. As it was shown in chapter 5 
the Sunrise detector shows an excellent performance in comparison with the ones 
installed in the Sunset side and thus its potential for axion discovery is higher. An 
optimization of the software efficiency and the background level of the detectors has 
been performed, in which the figure of merit of the detectors has been maximized. 
The background and tracking levels of the Micromegas detectors during 2011 are 
presented in table 10.1. As shown, both levels are compatible, in addition the 
occurrence of the events during background and tracking are compatible with the 
expected Poissonian. 


Detector 

Number of 
trackings 

Background level 
c cm 1 s 1 keV 1 

Tracking level 
c cm -1 s -1 keV 1 

Sunrise 

46 

(6.09 ± 0.10)xl0 -6 

(5.71 ± 0.55)xl0~ 6 

Sunset 1 

45 

(5.96 ± 0.10) xl0~ 6 

(6.14 ± 0.57) xlO -6 

Sunset2 

45 

(6.83 ± 0.11)xl0 -6 

(7.58 ± 0.63)xl0~ 6 


Table 10.1: Summary of tracking and background data for all the three Micromegas detec¬ 
tors during the 2011 data taking campaign. Background and tracking levels are computed 
from 2-7 keV and inside the coldbore area (14.52 cm 2 ). 


In order to distinguish the presence of an axion signal in the Micromegas data, an 
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unbinned likelihood method has been developed. In a first stage, for every axion 
mass, the log-likelihood is computed for several values of the coupling constant g or/ , 
in which the minimum and the standard deviation are extracted. If the minimum 
is compatible with the absence of signal (see figure 10.6), an upper limit to the 
coupling constant may be derived. 



Figure 10.6: Left: y 2 as a function of gf n for an axion mass of m„ = 1.061 eV (black 
line). The dashed blue line marks the distance at one unit above the minimum of y 2 and 
the red lines represent the positions of the points at one a. Right: Deviation in a units of 
the minimum of y 2 from zero. 


After discarding a possible axion signal, a limit on the axion to photon coupling is 
extracted. An upper limit to the coupling constant for a given axion mass has been 
calculated by integrating the Bayesian probability at a 95% of C.L. The derived 
limit for the 2011 Micromegas data at the CAST experiment is shown in figure 10.7. 

In order to estimate the systematic uncertainties, the determination of the effective 
coherence length has been considered as the main source of error. Indeed, the 3 He 
dynamics inside the magnet bores has a complex behavior, due to hydrostatic and 
convection effects at a cryogenic temperature. In order to calculate the effective 
length, a detailed study of the density profile of the 3 He has been performed using 
CFD 1 simulations. The effective coherence length has been parameterized for a given 
pressure inside the magnet bores. However, the profile density along the magnet 
cannot be measured and the simulation has some uncertainties. Two extreme cases 
have been taken into account in order to estimate the systematics: one considering all 
the magnet length (9.26 m) as the coherence region and another using a conservative 
scenario for the coherence region. The average values of the coupling limit for the 
nominal analysis and the ones described above for the systematics, for 1 < m a < 
1.17 eV, have been computed: 


1 Computational Fluid Dynamics 
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Figure 10.7: Axion to photon coupling upper limit for the Micromegas detectors during 
2011 data taking campaign (red line). The green line corresponds to the KSVZ model with 

E/N = 0. 


g ai < 4.04 x 10 -10 GeV -1 

Pessimistic 


g ai < 3.90 x 10 _10 GeV _1 

Nominal 

(10.1) 

g ai < 3.65 x 10 -10 GeV -1 

Optimistic 



Using the unbinned likelihood method, the data of the 3 He phase of the CAST 
experiment, from 2008 to 2011, have been computed and a coupling limit has been 
extracted. The results are shown in figure 10.8. 

CAST finished the 3 He phase extending the axion search up to m a = 1.17 eV, 
crossing first time the KSVZ benchmark (E/N—0) model, one of the most favored 
in the theory. Although there was no axion signal, CAST is currently rescanning 
the vacuum phase with an enhanced sensitivity, due to the improved background 
level in the Micromegas detectors. 

CAST microbulk Micromegas exploit different techniques in order to improve its 
background level: the intrinsic radiopurity; the improvements in the manufacturing 
process; the event discrimination and the shielding strategy. These strategies have 
been described in chapter 7. In order to understand the origin of the background 
level of the detectors, different special set-ups have been developed: one in the 
LSC and another in the Zaragoza laboratory. For the underground measurements, 
the muon flux is highly suppressed and its contribution to the background can 
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Figure 10.8: Parameter space m 0 —c/ a7 excluded by the CAST experiment for the vacuum 
and 4 He phase (black line) and the latest 3 He phase (red line). 


be considered negligible, with the proper shielding the background is reduced to a 
~ 10 -7 c cm' 2 keV -1 s _1 level. Moreover, different contributions to the background 
have been measured in this set-up, like the A1 strongback or the effect of the 222 Rn 
(see figure 10.9 left). On the other hand, the measurements performed at surface 
level were crucial in order to measure the contribution of the muons to the final 
background and an important test bench in order to develop the DAQ with the novel 
AFTER front end electronics. These special set-ups together with the simulations 
(see figure 10.9 right) have set the roadmap to the different upgrades developed for 
the Micromegas detectors at CAST. 

Following the low background strategies developed for the Micromegas detectors, the 
Sunset detectors were upgraded during 2012. A new shielding design was proposed 
with an inner copper shielding of 10 mm and a lead shielding of 100 mm (see 
figure 10.10 left), which has been extended along the magnet bore pipes. Moreover, 
the connection to the magnet bores is done by a 10 mm thick copper pipe, which has 
an inner PTFE coating with a thickness of 2.5 mm in order to attenuate the copper 
fluorescence at 8 keV. Additionally, two plastic scintillators were installed around the 
shielding (see figure 10.10 right) in order to discriminate events related with muons. 
Also, during 2013 the novel AFTER front-end electronics for the strip readout were 
installed, which lead to a better discrimination of the background events. 

During the 2014 CAST data taking campaign an X-ray focusing device was installed 
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Date shielding thickness (cm) 

Figure 10.9: Left: Background level evolution at the LSC. An ultra-low background 
level (red anows) was interrupted due to the measurements of the aluminum cathode (blue 
anow) and the 222 Rn intrusion (green an'ow). Right: Comparison between experimen¬ 
tal background levels as a function of the lead shielding thickness in special set-ups and 
simulations. The black line corresponds to a tit of the simulated data (black circles) for en¬ 
vironmental 7’s. The red squares correspond to the measurements performed underground 
and the blue squares are the measurements at surface using an active a muon veto. 



Figure 10.10: Left: Photo of the Sunset shielding partially opened in which the inner 
copper layer and the external lead are visible. Right: Scheme of the Sunset Micromegas 
shielding design in which the plastic scintillators are shown. 


in the Sunrise side, with a Micromegas detector in its focal plane (see figure 10.11). 
The Sunrise Micromegas detector has a novel design, the body and the chamber of 
the detector is made of 20 mm thick radiopure copper and all the gaskets are made 
of PTFE. Following the Sunset design, the shielding has a thickness of 100 mm 
of lead, which is extended along the magnet bore pipe and a plastic scintillator is 
installed on the top for the rejection of the events induced by muons. Furthermore 
a new Micromegas detector has been manufactured for the new line, being the one 
with the better performance working at CAST with a 13% of FWHM in the 5.9 keV 
peak and an excellent spatial resolution and homogeneity of the gain in the active 


area. 
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Figure 10.11: Photo of the new Sunrise Micromegas+XRT system in the CAST experi¬ 
ment. The different parts of the set-up are labeled. 


The upgrades of the Micromegas detectors at CAST have led to a reduction of 
the background level down to ~ 10~ 6 c cm~ 2 keV -1 s -1 . Moreover, the new 
XRT+Micromegas line in the CAST experiment has set a milestone in axion re¬ 
search. CAST will finish the rescanned vacuum phase at the end of 2015 improving 
its previous limit to an expected value of ~ 6 x 10 -10 GeV -1 , due to the reduction 
of the background levels of the detectors and the new XRT. However, the discovery 
potential of CAST is limited by the size of the magnet. In order to scan a wider 
region of the m a — g a -y parameter space the IAXO experiment has been proposed. 

Beyond CAST a new generation helioscope with a improved sensitivity, specifi¬ 
cally built for axion and ALPs searches, has been proposed: IAXO-the Interna¬ 
tional AXion Observatory. It will enhance the helioscope technique by exploiting 
all the singularities of CAST, implemented into a large superconducting toroidal 
magnet, together with X-ray optics and ultra-low background detectors attached 
to the magnet bore ends. IAXO detectors have a goal of background levels of 
~ 10 -7 c cm -2 keV -1 s^ 1 (see figure 10.12) and down to ~ 10~ 8 if possible. 

The new system (Micromegas+XRT) at CAST, can be considered as a IAXO 
pathfinder, being an important milestone for the technical design phase of IAXO. 
However, the final background is around one order of magnitude above the levels 
required for IAXO. In this way, the design of a DO detector for IAXO has been 
proposed, it will feed by the low background techniques developed until now and 
also new research lines have been proposed, such as the improvement of the muon 
veto coverage, new thin cathode windows, new gas mixtures, the novel AGET front 
end electronics and resistive Micromegas. 
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Figure 10.12: Micromegas background history since 2002. The blue dots correspond to 
the unshielded Micromegas with the classic technology at CAST, the green squares are the 
background levels at CAST after the installation of the shielding, the black triangles cor¬ 
respond to the shielding upgrade from 2012 to 2014 and the red dots show the background 
levels reached at the LSC. The dashed black line marks the goal for IAXO. 


Due to a dedicated magnet, optics and ultra-low background detectors, IAXO 
will surpass the sensitivity of CAST in more than one order of magnitude 
(see figure 10.13). IAXO will be sensitive to coupling constants of about 
Sa'y ~ 5 x 1CT 12 GeV -1 , entering into an unexplored parameter space area and by 
first time in a favored region for axions and ALPs. 

IAXO could directly measure the solar flux of axions produced by non-hadronic 
processes, for the first time with sensitivity to relevant g ae values. However, in this 
case the reduction of the low energy threshold in the detectors and the increment of 
the efficiency at low energies will be crucial. On the other hand, the huge magnet 
required for IAXO offers excellent possibilities to host relic dark matter searches 
of axions and ALPs and the use of a dish antenna and resonant cavities are under 
study. 

The haloscope technique consists in a tunable microwave-cavity in a strong magnetic 
held coupled to an ultra-low-noise microwave sensor. CDM axions and ALPs may 
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Figure 10.13: Expected sensitivity of IAXO (black filled areas) for hadronic axions, 
compared with current bounds from CAST and ADMX. Also, future prospects of ADMX 
(dashed brown region) and ALPS-II (light blue line) are shown. 

convert into photons in the microwave regime with an enhanced probability if the 
resonant frequency of the cavity matches the axion energy. A directional effect on 
the detection could be observed by the use of long thin cavities properly tuned. 
This case has been studied in chapter 9 and could provide a strong identihcative 
signature of the direction of the CDM axions. The expected sensitivity of this kind 
of resonant cavities is shown in figure 10.14. 

The use of thin cavities inside long magnets for relic axion detection is particularly 
appealing, because it could be realized in the near future, given that this type of 
magnets are already used by the axion community in experiments looking for solar 
axions, like CAST. Moreover, the huge magnet required for IAXO offers excellent 
possibilities and the hosting of this kind of cavities is under study. 
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Figure 10.14: Sensitivity regions for long thin microwaves cavities, with a geometry fixed 
by the desired directional effect. Two different scenarios has been computed: a conserva¬ 
tive one (black region) and the optimistic one (red region). 
























































































































Chapter 11 


Resumen y conclusiones 


Los axiones son partfculas pseudo-escalares propuestas en una extension del Modelo 
Estandar (SM) como una solucion al problema the CP en las interacciones fuertes. 
La teoria precide una violation de CP en las interacciones fuertes que no ha sido 
observada experimentalmente. La solucion mas convincente a este problema fue 
propuesta por Peccei y Quinn en 1977, en ella se introduce una simetria global 
y quiral U(l)pq que se rompe espontaneamente en la escala de energia f a de la 
simetria. De esta manera se soluciona el problema CP dinamicamente y aparece el 
axion, como en boson pseudoescalar de Nambu-Goldstone de la nueva simetria. 

La solucion de Peccei-Quinn fija algunas propiedades de los axiones como su masa y 
su constante de acoplo, que estan relacionadas con la escala de energia de la nueva 
simetria, f a . Los axiones pueden interacturar con gluones, fotones y fermiones. Sin 
embargo, el acoplo mas interesante es el de axion-foton, que es generico para todos 
los modelos. 

Aun cuando los axiones son las partfculas mejor motivadas en la teoria, tambien 
existe la categoria de partfculas tipo axion (ALPs) o mas genericamente WISP 
(Weakly Interacting Slim Particles). Las cuales comparten la misma fenomenologfa 
que el axion, siendo partfculas ligeras que se acoplan a dos fotones. Las ALPs 
emergen de extensiones del SM donde una nueva simetria se rompe a altas escalas 
de energia, tambien aparecen en la teoria de cuerdas al igual que el axion. Sin 
embargo en este caso la constante de acoplo no esta relacionada con su masa, de esta 
manera las ALPs pueden encontrarse en extensas regiones del espacio de parametros 
Sa'y a- 

Tanto los axiones como las ALPs podrfan haberse producido en un Universo prim¬ 
itive mediante mecanismos no termicos, como el realineamiento del vacio o el de- 
caimiento de los defectos topologicos. Siendo partfculas neutras que interaccionan 
debilmente con la materia, tanto los axiones como las ALPs son atractivos can¬ 
didates a Materia Oscura (DM), que podrfan explicar separadamente toda la can- 
tidad de DM en el Universo. Sin embargo las propiedades de los axiones y de las 
ALPs estan restringidas debido a consideraciones provenientes de la astrofisica y de 
la cosmologia, ya que podrfan tener un papel importante en la evolution estelar. Por 
otra parte diferentes observaciones experiment ales podrfan ser interpretadas como 
un indicio de axiones o ALPs. Como la excesiva transparencia del Universo a fotones 
de alta energia o el enfriamiento anomalo de las Enanas Blancas. 

Diferentes tecnicas han sido desarrolladas en la busqueda de axiones: helioscopios 
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que buscan axiones solares, haloscopios que buscan axiones primigeneos parte de la 
DM y experimentos de regeneration en los que los axiones podrfan ser generados 
y detectados en el laboratorio. Todas ellas estan basadas en el effecto Primakoff, 
donde los axiones podrlan ser convertidos en fotones en presencia de campos elec- 
tromagneticos. 



Figure 11.1: Esquema de deteccion de axiones solares en la tecnica del helioscopio. 




Figure 11.2: Izquierda: Espectro diferencial de los axiones solares en la Tierra para 
modelos hadronicos (llnea negra) y no hadronicos (llnea roja). Derecha: Probabilidad de 
conversion axion-foton en funcion de la masa del axion para el caso de vacfo (lfnea negra) 
y el caso de gas de buffer (lrnea roja). 

La tecnica del helioscopio, la cual es el tema principal de este trabajo, fue propuesta 
por Sikivie en 1983, en ella se utiliza al Sol como una poderosa fuente de axiones. Los 
axiones solares podrlan generarse en el nucleo solar en los fuertes campos electricos 
de plasma solar por medio del efecto Primakoff inverso. Ademas, la emision de 
axiones por procesos adicionales como Bremsstrahlung de axiones, Compton y axio- 
desexcitacion podrla ser considerablemente mas importante que la de los procesos 
Primakoff. Estos axiones podrlan ser reconvertidos en fotones en presencia de fuertes 
campos magneticos mediante el efecto Primakoff (ver figura 11.1). Estos fotones, que 
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estan en el rango de energfa de los rayos-X (ver la parte izquierda de la figura 11.2), 
podrfan ser detectados en detectores de rayos-X colocados en los extremos del iman. 

Aunque la emision de axiones solares mediante procesos no hadronicos podrla ser 
mas importante, en este trabajo solamete se tienen en cuenta los procesos Primakoff. 
Debido a que es mas consistente suponer los mismos procesos en la generacion y 
deteccion de axiones. Ademas, las restricciones provenientes de la astroffsica en g ae 
son mas restrictivas que los lhnites que en principio pueden ser alcanzados por los 
helioscopios. El flujo de axiones solares esta bien establecido debido a la detallada 
descripcion del Modelo Solar. Por otra parte la probabilidad de conversion axion- 
foton dentro de fuertes campos magneticos fue introducida por van Bibber. En 
la tecnica del helioscopio, la sensibilidad a la masa del axion viene dada por la 
probabilidad de conversion, la cual esta restringida a masas pequenas en el caso de 
que el iman este en vacfo (ver la parte derecha de la figura 11.2). Sin embargo, la 
coherencia puede ser restaurada a mayores masas de axiones utilizando gas de buffer 
dentro del iman. 

El mayor exponente dentro de la tecnica de los helioscopios es el experimento CAST 
(CERN Solar Axion Telescope), que busca axiones solares desde 2003, siendo el 
helioscopio mas sensible hasta la fecha. Para ello utiliza un iman obsoleto del LHC 
de tipo dipolo, con una longitud de 9.26 m y con campos magneticos hasta 9 T. El 
iman esta montado en una plataforma movil, la cual permite apuntar al Sol durante 
~ 1.5 h dos veces al dfa, durante la salida y la puesta del Sol (ver figura 11.3). El 
iman esta compuesto por dos cavidades magneticas con un total de cuatro detectores 
de rayos-X colocados en sus extremos. La serial de axiones serf a un exceso de rayos-X 
cuando el iman apunta al Sol y el uso de detectores de bajo fondo es fundamental. 



Sunset 

photon detectors 


"Sunrise axions 


Low-background 

shielding 


Figure 11.3: Esquema del experimento CAST donde las diferentes partes del iman y la 
plataforma movil han sido etiquetadas. 
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El experimento CAST esta compuesto por varios sistemas adicionales que se re- 
quieren para la toma de datos. El iman opera a una temperatura nominal del 1.8 Iv, 
para lo que se requiere un sistema criogenico de enfriamiento. Tambien un sistema 
de vacio ha sido instalado alrededor the las cavidades del iman, con el proposito 
de incrementar la transparencia de los rayos-X provenientes de los axiones y para 
aislar el iman del ambiente. El movimiento del iman se controla mediante un pro- 
grama de seguimiento, el cual permite apuntar automaticamente al Sol. Ademas, 
un sistema de gas ha sido instalado con el proposito de rellenar las cavidades del 
iman con un gas de buffer en pequenas cantidades. Finalmente todos los sistemas 
son monitorizados mediante el slow control. 

El programa de investigation de CAST se puede dividir en dos fases, la primera con 
vacio dentro de las cavidades del iman y una segunda fase utilizando gas de buffer. 
El periodo de toma de datos de la Fase I, se realizo durante 2003 y 2004, obteniendo 
un limite en la constante de acoplo de g ai < 8.8 X 10” 11 GeV” 1 a un 95% de nivel 
de confianza (C.L.), para masas de axiones rn a < 0.02 eV. 

Durante la Fase II, de 2005 a 2006 el iman se relleno de 4 He, obteniendo un 
limite de g ai < 2.17 x 10” 10 GeV -1 con un 95% de CL, para masas de axion 
0.02 < m a < 0.39 eV. Mientras que la fase de 3 He se inicio en 2008 y final- 
izo en 2011, donde se escanearon masas de axiones 0.39 < m a < 1.17 eV, obte¬ 
niendo un limite en la constante de acoplo de g ai < 2.3 x 10” 10 GeV -1 para 
0.39 < m a < 0.64 eV y g ai < 3.3 x 10” 10 GeV” 1 para 0.64 < m a < 1.17 eV 
con un 95% de C.L. Cruzando por primera vez el modelo de referenda KSVZ, uno 
de los mas favorecidos por la teoria. Aunque CAST finalizo su programa de in¬ 
vestigacion en 2011, la toma de datos ha sido ampliada. En 2012 la fase de 4 He 
fue reescaneada, mejorando el limite anterior en una estrecha region. En 2013 una 
nueva fase de vacio se inicio en CAST, motivada por la mejora del nivel de fondo de 
los detectores, donde se espera una mejora en la sensibilidad del experimento. 

Tres tipos diferentes de detectores han operado en CAST desde el comienzo del 
experimento: un dispositivo de carga acoplada (CCD) en el piano focal de un tele- 
scopio de rayos-X, una camara de proyeccion temporal (TPC) que cubria las dos 
cavidades del iman en el lado de la puesta del Sol, que opero hasta 2006, finalmente 
diferentes tipos de detectores Micromegas (MICRO MEsh GAseous Structure) han 
estado tomando datos en el lado de la salida del Sol y mas recientemente en el lado 
de la puesta del Sol. 

Los detectores Micromegas fueron desarrollados por Giomataris en 1996. Son del 
tipo de ionization con dos regiones separadas por una rejilla o mesh: la region de 
conversion y la region de amplification. En la region de conversion las particulas 
interaccionantes ionizan el gas, generando cargas primarias, en esta region esta pre¬ 
sente un campo electrico debido al cual los electrones se derivan hacia la region 
de amplification. En la cual se produce la avalancha de los electrones primarios, 
debido al gran campo electrico presente en este gap (de unos ~ 50 /mi), generando 
dos senales medibles: una en la mesh, que proporciona resolution espacial y otra en 
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Figure 11.4: Esquema del principio de funcionamiento de un detector Micromegas. Dos 
regiones diferentes estan separadas por una rejilla o mesh: la region de conversion y la 
region de amplification. 


el piano de lectura del anodo, la cual confiere resolution espacial a los detectores. 
Un esquema sobre el principio de trabajo de los detectores Micromegas se muestra 
en la figura 11.4. 

Los detectores Micromegas son una tecnologia en alza, la cual ha experimentado una 
evolution tecnologica en el proceso de fabrication, desde las Micromegas clasicas 
hasta las tecnologias bulk y microbulk. El experimento CAST ha sido un deman- 
dante banco de pruebas para estas diferentes tecnicas de fabrication donde la nueva 
tecnologia microbulk es la que muestra un mejor rendimiento y la mayor radiop- 
ureza. Tres de los cuatro detectores instalados en CAST actualmente son de este 
tipo. 

El diseho de los detectores Micromegas del tipo microbulk instalados en CAST 
durante 2011 ha sido cuidadosamente seleccionado para la busqueda de axiones. El 
anodo del detector esta hecho de pistas cuadradas interconectadas. Esto lleva a 
un piano de lectura bidimensional de de 106 x 106 pistas con una distancia entre 
ellas de 550 fi m y un area efectiva de 60 x 60 mm 2 . El catodo esta hecho de mylar 
aluminizado de unos 5 /im de grosor, pegado a una rejilla de aluminio de forma 
circular, la cual se coloca al lado de vatio del iman. El cuerpo del detector esta 
hecho de plexiglass, donde la camara y la region de conversion tienen una altura 
de 30 mm. La Micromegas esta pegada a una base de plexiglass llamada raqueta 
que tiene una forma circular en la region del detector y un cuello estrecho para el 
paso de las conexiones de las pistas a la electronica. El gas en la camara es una 
mezcla de Ar + iC^iO, a una presion de 1.4 bares. Los detectores se colocan en 
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los extremos de las cavidades del iman, cubriendo toda la apertura de 14.52 cm 2 . 
Aunque los detectores del la salida y la puesta de Sol tienen diferentes disenos de 
blindaje, ambos estan compuestos principalmente de 5 mm de cobre en la parte mas 
interna, 25 mm de plomo arqueologico y una lamina de 2 mm de cadmio en la parte 
mas externa, cubierto por diferentes capas de polietileno. 

La adquisicion de los diferentes pianos de lecturas de los detectores Micromegas en 
CAST se realiza mediate un programa basado en Labview. El pulso de la mesh se 
digitaliza mediante una tarjeta Matacq y el piano de lectura del anodo se adquiere 
mediante las tarjetas Gassiplex. El analisis de los datos se realiza mediante un pro¬ 
grama dedicado, basado en C++ y en ROOT. Para la serial de la mesh se realiza un 
analisis de la forma del pulso, lo cual conlleva la definition de diferentes observables. 
Para el piano de lectura del anodo se realiza un analisis por clusters, lo cual permite 
la definition de diferentes parametros relacionados con la forma del evento. 

En una segunda etapa, se procede a la discrimination de los eventos tipo rayos-X. 
Para ello de compara la distribution de los eventos de fondo con las de los eventos de 
la calibration diaria con 55 Fe (ver parte izquierda de la figura 11.5). De esta manera 
se ha desarrollado un metodo de discrimination, computando la distribution de log- 
odds de diferentes observables, tanto para eventos de calibraciones como para los 
de fondo. Esto permite definir un valor de corte para el cual un cierto numero 
de eventos son aceptados (ver parte derecha de la figura 11.5), lo cual conlleva la 
definition de una eficiencia de software. 




Figure 11.5: Izquierda: Espectro de 5 5Fe en un detector Micromegas en el que se realiza 
el ajuste de los tres diferentes picos. Derecha: Distribicion de log-odds para el pico de 
5.9 keV (azul), pico de 3 keV (verde) y eventos de fondo (rojo). El area azul corresponde 
al corte en 5.9 keV, mientras que el area verde corresponde al corte a 3 keV. 


Utilizando este metodo de discrimination, se han analizado los datos de los detec¬ 
tores Micromegas en CAST durante 2011. El poder de discrimination de los detec¬ 
tores es fuertemente dependiente de la calidad del detector. Como se ha descrito en 
el capltulo 5, el detector de la salida del Sol muestra unas cualidades excelentes, en 
contraste con los detectores instalados en el lado de la puesta del Sol y su potential 
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para el descubrimiento de axiones es mayor. Se ha realizado una optimization de 
la eficiencia de software en conjunto con el nivel de fondo de los detectores, de esta 
manera se ha maximizado la figura de merito de los detectores y consecuentemente 
su potencial para el descubrimiento de axiones. Los niveles de fondo y durante 
el seguimiento al Sol para los detectores Micromegas en 2011 se muestran en la 
tabla 11.1. Como se puede observar ambos niveles son compatibles, ademas se ha 
comprobado que la ocurrencia de eventos tanto de fondo como de seguimiento al Sol 
sigue distribution de Poisson. 


Detector 

Number of 
trackings 

Background level 
c cm 1 s 1 keV 1 

Tracking level 
c cm~ 1 s“ 1 keV 1 

Sunrise 

46 

(6.09 ± 0.10) xlO" 6 

(5.71 ± 0.55)xl0~ 6 

Sunset 1 

45 

(5.96 ± 0.10)xl0 -6 

(6.14 ± 0.57)xl0 -6 

Sunset2 

45 

(6.83 ± 0.11)xl0 -6 

(7.58 ± 0.63) xlO -6 


Table 11.1: Niveles de fondo y de seguimiento al Sol para los tres detectores Micromegas 
durante la toma de datos de 2011. Los diferentes niveles se calculan en el rango de energias 
de 2-7 keV y dentro del area de la cavidad del iman (14.52 cm 2 ). 

Con el objetivo de discriminar la presencia de una serial se ha desarrollado un metodo 
unbinned likelihood. En una primera fase, para cada masa de axion, se calcula la 
log-likelihood para diferentes valores de la constante de acoplo g ai , donde se extrae 
el minimo y la desviacion estandar. Si el nhnimo es compatible con la ausencia de 
serial (ver figura 11.6), se puede extraer un lhnite en la constante de acoplo. 



Figure 11.6: Izquierda: y 2 en funcion de para una masa de axion de m a = 1.061 eV 
(linea negra). La linea azul punteada marca la distancia de una unidad por encima del 
minimo de y 2 . Derecha: Desviacion en a del minimo de y 2 desde cero. 

Despues de descartar una posible serial de axiones, se extrae un lhnite en la constante 
axion-foton. Donde se ha calculado un lhnite superior para una masa de axion dada 
mediante la integration de la probabilidad Bayesiana con un 95% de C.L. El lhnite 
extraido para los datos de los detectores Micromegas durante la campana de 2011 
se muestra en la figura 11.7. 
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Figure 11.7: Li'mite superior de la constante de acoplo axion-foton para los detectores 
Micromegas durante la campana de toma de datos de 2011 (lmea roja). La lrnea verde 
corresponde al modelo KSVZ con ( E/N = 0). 


Con el objetivo de estimar el error sistematico en el limite extraido se considera la 
determination de la longitud efectiva de coherencia como la mayor fuente de error. 
La dinamica del 3 He en las cavidades del iman tiene un comportamiento complejo 
debido a efectos hidrostaticos y de convection a una temperatura criogenica. Con el 
objetivo de calendar la longitud efectiva del iman se ha realizado un estudio detallado 
del perfil de densidades dentro del iman utilizando simulaciones CFD 1 . De esta 
manera la longitud de coherencia en el iman ha sido parametrizada para una presion 
dada. Sin embargo, el perfil de densidades a lo largo del iman no puede ser medido 
y la simulation tiene ciertas incertidumbres. Dos casos extremos se han tenido en 
cuenta para estimar el error sistematico en la constante de acoplo: uno de ellos 
considerando como longitud de coherencia toda la longitud del iman (9.26 m) y otro 
utilizando un escenario conservador para la region de coherencia. El valor promedio 
de las constantes de acoplo para el analisis nominal y para los descritos anteriormente 
para el calculo de sistematicos, en el rango de masas 1 < m a < 1.17 eV, se muestran 
a continuation: 


g ai < 4.04 x 10 _10 GeV _1 

Pesimista 


g ai < 3.90 x 10 _1 °GeV _1 

Nominal 

(11.1) 

g ai < 3.65 x 10 _1 °GeV _1 

Optimista 



1 Computacion de Dinamica de Fluidos 
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Utilizando el metodo unbinned likelihood se han analizado los datos de toda la fase 
de de 3 He, de 2008 a 2011, extrayendo un lhnite en la constante de acoplo. Los 
resultados se muestran en la figura 11.8. 



Figure 11.8: Espacio de parametros rn a — g nrf exclufdo por el experimento CAST para 
las fases de vacfo y de 4 He (lfnea negra) y los ultimos resultados de la fase de 3 He (lfnea 
roja). 

CAST finalizo la fase de 3 He extendiendo la busqueda de axiones hasta 
rn a = 1.17 eV, cruzando por primera vez el modelo de referenda KSVZ ( E/N = 0), 
uno de los mas farorecidos en la teorla. Aunque no hubo ninguna serial de axiones, 
CAST actualmente esta reescaneando la fase de vaclo con una sensibilidad mejorada, 
debido a la reduccion del nivel de fondo en los detectores Micromegas. 

Los detectores Micromegas del tipo microbulk explotan diferentes tecnicas para la 
reduccion del nivel de fondo: la radiopureza intrmseca; las mejoras en el proceso de 
fabrication; la discrimination de eventos de fondo y el blindaje. Estas estrategias 
han sido descritas en el capitulo 7. Con el proposito de entender el origen del fondo 
de los detectores se han desarrollado diferentes montajes: uno en el LSC y otro 
en el laboratorio de Zaragoza. Para las medidas bajo tierra, el flujo de muones 
esta altamente suprimido y su contribution se puede considerar despreciable. Con 
el blindaje adecuado el fondo se reduce a un nivel de ~ ICC 7 c cm -2 keV -1 s _1 
(ver parte izquierda de la figura 11.9). Adernas diferentes contribuciones al fondo 
de los detectores han sido medidas, como la del catodo de aluminio o el efecto del 
222 Rn en el montaje. Por otra parte, las medidas realizadas en superficie fueron 










180 


Chapter 11. Resumen y conclusiones 


cruciales para medir la contribution de los muones al fondo de los detectores y 
un importante banco de pruebas para el desarrollo de la adquisicion con la nueva 
electronica AFTER. Estos montajes especiales junto con las simulaciones (ver parte 
derecha de la figura 11.9) marcaron la hoja de ruta para las diferentes mejoras 
desarrolladas para los detectores Micromegas en CAST. 



Date 



shielding thickness (cm) 


Figure 11.9: Izquierda: Evolution del nivel de fondo en el LSC. Donde los niveles de 
ultra-bajo fondo (flechas rojas) fueron interrumpidos debido a las medidas con el catodo 
de aluminio (flecha azul) y la intrusion de 222 Rn (flecha verde). Derecha: Comparacion de 
los niveles de fondo experimentales y simulaciones, en funcion del grosor de plomo. La 
lmea negra corresponde al ajuste de los datos simulados (tirculos rojos) para 7 s del fondo 
natural, mientras que los cuadrados rojos son las medidas bajo tierra y los cfrculos azules 
son las medidas en superficie con un veto de muones. 



Figure 11.10: Izquierda: Foto de los detectores de la puesta del Sol con el blindaje 
parcialmente abierto, en el que se muestran las capas de cobre interno y el plomo extemo. 
Derecha: Esquema del diseno del blindaje de los detectores de la salida del Sol donde se 
puede observar la position de los plasticos centelleadores. 

Siguiendo las estrategias de fondo desarrolladas para los detectores Micromegas, los 
detectores de la puesta de Sol fueron mejorados en 2012 . Con un nuevo diseno del 
blindaje con un capa interna de 10 mm de cobre y una capa externa de 100 mm de 
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plomo (ver parte izquierda de la figura 11.10), la cual se ha extendido a traves de las 
tuberlas hacia las cavidades magneticas. Ademas, la conexion hacia estas cavidades 
se realiza mediante una tuberla de cobre de 10 mm de grosor con un recubrimiento 
de PTFE de 2.5 mm para atenuar la fluorescencia del cobre a 8 keV. Tambien se han 
instalado dos centelleadores plasticos cubriendo en blindaje (ver parte derecha de 
la figura 11.10), con el proposito de discriminar eventos relacionados con muones. 
Ademas, durante 2013 se instalo la nueva electronica AFTER para el piano de 
lectura del anodo, que conlleva a una mejor discrimination de los sucesos de fondo. 

Durante la toma de datos de CAST correspondiente al ano 2014, se ha instalado 
en el lado de la salida del Sol un nuevo telescopio de rayos-X con una Micromegas 
en su piano focal (ver figura 11.11). El detector Micromegas de la salida del Sol 
tiene un nuevo diseno, en el cual la raqueta y la camara estan compuestas por 
cobre radiopuro de un grosor de 20 mm, ademas todas las juntas estan hechas de 
PTFE. Tambien siguiendo el diseno de los detectores de la puesta de Sol, el blindaje 
de plomo tiene un grosor de 100 mm, el cual se extiende a traves de la tuberla 
hacia el iman, tambien se ha instalado un centelleador plastico encima del blindaje 
para discriminar los eventos inducidos por muones. Ademas, un nuevo detector 
Micromegas se ha fabricado para la nueva llnea, siendo uno de los detectores con 
las mejores cualidades instalados en CAST, con un 13% de resolucion FWHM en el 
pico de 5.9 keV y una excelente resolucion espacial y homogeneidad de la ganancia 
en la zona activa. 



Figure 11.11: Fotografla de la nueva llnea del experimento CAST, compuesta por un 
detector Micromegas y un telescopio de rayos-X. Las diferentes partes han sido etiquetadas. 


Las mejoras de los detectores Micromegas en el experimento CAST han conllevado 
una reduction del fondo a unos niveles de ~ 10 -6 c cnL 2 keV -1 s~ 1 . Ademas, la 
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nueva lfnea del XRT+Micromegas de CAST supone un hito en la investigation de 
axiones. CAST acabara la nueva fase de vacfo en 2015 mejorando el anterior lfmite, 
debido a la reduction del fondo de los detectores y el nuevo telescopio de rayos- 
X. Donde se espera una reduction del lnnite en la constante de acoplo a valores 
sobre ~ 6 X 10 -10 GeV -1 . Sin embargo, la sensibilidad de CAST esta limitada 
principalmente por el tarnano del iman. Con el proposito de escanear una region 
mas extensa del espacio de parametros m a — g a 7 , se ha propuesto el experimento 
IAXO. 
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Figure 11.12: Evolucion de los niveles de fondo desde 2002 en los detectores Mi- 
cromegas. Los puntos azules corresponden a Micromegas sin blindaje con la tecnologfa 
clasica, los cuadrados azules representan el nivel de fondo despues de la instalacion del 
blindaje, los triangulos negros corresponden a las mejoras de blindaje de 2012 a 2014, 
mientras que los puntos rojos marcan en nivel de fondo alcanzado en el LSC. La lfnea 
discontinua marca la meta de IAXO. 


Mas alia de CAST, se ha propuesto un helioscopio de nueva generation, con una 
sensibilidad mejorada y espetificamente construido para la busqueda de axiones y 
ALPs: IAXO-el Observatorio International de Axiones. El cual mejorara la tecnica 
del helioscopio explotando todas las peculiaridades de CAST, implementadas en un 
gran iman toroidal superconductor, junto con nuevas opticas de rayos-X y detectores 
de ultra-bajo fondo situados en los extremos del iman. Con el objetivo de alcanzar 
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niveles de fondo de ~ 10 ' c cm 2 keV 1 s 1 (ver figura 11.12) y hasta ~ 10 8 si 
es posible. 
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Figure 11.13: Sensibilidad esperada para IAXO en el caso de axiones hadronicos (re- 
giones rellenas de negro), comparada con los actuales h'mites para CAST y ADMX. Tam- 
bien se muestran las futuras perspectivas para ADMX (regiones marrones) y ALPS-II (lfnea 
azul claro). 

El nuevo sistema (Micromegas+XRT) de CAST, puede ser considerado como pionero 
de IAXO, siendo un hito importante para la fase de diseno tecnico. Sin embargo, 
el nivel de fondo en CAST esta un orden de magnitud por encima de los niveles 
requeridos para IAXO. De este modo se ha propuesto el desarrollo de un detec¬ 
tor DO para IAXO, el cual se proveera de las tecnicas de bajo fondo desarrolladas 
hasta ahora, ademas nuevas mejoras y lineas de investigacion han sido propuestas, 
tales como: mejoras en la cobertura de los vetos de muones, nuevas ventanas para 
el catodo, nuevas mezclas de gases, la novedosa electronica AGET y Micromegas 
resistivas. 

Gracias a un iman dedicado, optica y detectores de ultra-bajo fondo, IAXO superara 
la sensibilidad de CAST en mas de un orden de magnitud (ver figura 11.13). Siendo 
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sensible a constantes de acople entorno a g a -y ~ 5 X 10~ 12 GeV^ 1 . Entrando en 
un area inexplorarada del espacio de parametros y por primera vez en una region 
favorable para axiones y ALPs. 



Figure 11.14: Regiones de sensibilidad para cavidades resonantes estrechas y largas, 
donde la geometria ha sido fijada para conseguir el efecto de direccionalidad deseado. Dos 
diferentes escenarios han sido calculados: uno conservativo (area negra) y otro optimista 
(area roja). 

IAXO podria medir directamente el flujo de axiones solares producidos en procesos 
no hadronicos y ser por primera vez sensible a valores de g ae relevantes. Sin embargo, 
en este caso sera crucial la reduction del umbral de energia de los detectores y 
el aumento de eficiencia a bajas energias. Por otra parte, el inmenso iman que 
se requiere para IAXO ofrece excelentes posibilidades para albergar busquedas de 
axiones y ALPs componentes de DM, de esta forma se estan estudiando el uso de 
una antena parabolica (dish antenna ) y cavidades resonantes. 

La tecnica del haloscopio consiste en una cavidad de microondas ajustable dentro 
de un fuerte campo magnetico, acoplada a sensores de microondas de ultra-bajo 
fondo. Axiones y ALPs CDM podrian convertirse en fotones en el regimen de 
las microondas con una probabilidad mejorada en el caso de que la frecuencia de 
resonancia de la cavidad coincida con la energia del axiom Se podria observar un 
efecto directional mediante el uso de cavidades estrechas y largas apropiadamente 
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ajustadas. Este caso ha sido estudiado en el capltulo 9 y podrla proporcionar una 
fuerte serial de identification de los axiones CDM. La sensibilidad de este tipo de 
cavidades resonantes ha sido estudiada y se muestra en la hgura 11.14. 

El uso de cavidades estrechas dentro de imanes largos para la detection de axiones 
reliquia es un caso particularmente atractivo, ya que se podrla realizar en un futuro 
proximo, dado que este tipo de imanes estan siendo utilizados en la comunidad 
de axiones, en experimentos que buscan axiones solares, conro CAST. Adenras, el 
enorme iman que se requiere para IAXO ofrece excelentes posibilidades y se esta 
estudiando la posibilidad de albergar de cavidades resonantes en el rango de las 
microondas. 
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